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ABSTRACT 
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The  effect  of  cloud  cover  and  surface  albedo  on  the  solar 
insolation  at  the  surface  of  the  Earth  was  investigated.  The 
model  uses  Bird's  (1984)  clear  sky  model  for  directinormal 
and  diffuse  horizontal  irradiance  as  its  basis.  Bird's  model 
was  modified  to  include  inclined  surfaces  and  cloud  cover. 
The  clouds  were  modelled  using  parameter izations  developed 
by  Welch,  Cox  and  Davis  (1980)  and  a  weighting  technique  to 
account  for  cloud  shape  using  plane/parallel  calculations  ■  - 
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suggested  by  Welch  and  Wielicki  (1984).  The  radiative 
transfer  in  the  cloud  was  modelled  by  a  two  stream 

O 

approximation  known  as  the  Shuster-Schwarzchild  method.  The 
model  results  were  compared  with  clear  sky  broadband  and 
spectral  observations,  and  cloudy  broadband  observations. 

The  clear  sky  spectral  observations  were  obtained  from  the 
Solar  Energy  Research  Institute  (SERI)  in  Golden,  Colorado. 
The  broadband  observations  were  obtained  from  the  Solar  and 
Meteorological  Research  Program  at  the  University  of 
Michigari.NThe  effects  of  cloud  amount,  cloud  thickness, 
cloud  type,  and  surface  albedo  are  demonstrated.  Cloud 
amount  and  cloud  thickness  are  the  predominant  cloud  factors 
affecting  insolation.  High  surface  albedos  tend  to  slow  the 
decrease  in  insolation  due  to  increasing  cloud  cover;  while 
low  surface  albedos  cause  greater  variation  in  insolation 
under  the  same  circumstances.^The  model  has  the  capability 
to  simulate  simple  mountain  obstructions.  However,  no  cases 
with  obstructions  are  presented.  The  model  produces 
terrestrial  spectra  between  0.28  and  4.0  micrometers  with  a 
resolution  of  approximately  10  nm.  The  goal  of  this  work  is 
to  provide  researchers  with  the  capability  to  calculate 
spectral  irradiance  for  different  atmospheric  conditions  and 
different  collecting  surface  geometries. 
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CHAPTER  I 


INTRODUCTION 


The  purpose  of  this  research  is  to  develop  a  solar 
spectral  model  that  describes  the  instantaneous  and  daily 
solar  radiation  flux  for  a  given  location  on  the  Earth.  The 
model  was  designed  so  the  results  would  be  valid  at  the 
surface  of  the  Earth  and  the  top  of  the  atmosphere.  The 
author's  desire  was  to  construct  a  solar  spectral  model 
which  would  compute  the  global,  direct,  diffuse,  and 
direct-normal  radiation  incident  on  the  surface.  It  was  also 
desirable  to  obtain  the  spectral  fluxes  of  radiation  leaving 
the  top  of  the  atmosphere  for  those  individuals  interested 
in  remote  sensing  applications.  With  this  purpose  in  mind, 
several  types  of  models  were  studied. 

Many  of  the  models  for  solar  radiation  were  designed 
with  a  specific  purpose  in  mind.  A  model  developed  by  Szeicz 
(1974)  is  specialized  for  use  in  plant  growth.  Other  models 
were  designed  for  particular  locations  (Barbaro,  1979), 
clear  sky  cases  only  (Bird,  1984),  for  slopes  of  different 
orientation  (Temps  and  Coulson,  1977;  Sweat  and  Carroll, 
1983),  and  for  overcast  skies  only  (Atwater  and  Ball,  1981; 
Hasten  and  Dzeplak,  1980;  Liou,  1976).  Some  models  were 
statistical  correlations  between  observed  parameters  and  the 
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solar  radiation  measured  at  the  surface  (Stanhill,  1983). 
Some  models  dealt  only  with  the  ultraviolet  portion  of  the 
solar  spectrum  (Leighton,  1961;  Green,  1974).  There  were 
models  that  used  a  multilayered  atmosphere  (Dave,  1978;  Dave 
and  Brauslau,  1974)  while  others  used  only  a  single  layer 
(Bird,  1984).  In  some  cases,  the  models  were  too  restrictive 
to  be  used  in  a  general  sense  and  others  were  too  difficult 
to  be  used  by  researchers  outside  the  physical  sciences  due 
to  their  complexity. 

There  are  many  parameters  and  factors  which  should  be 
considered  when  developing  a  general  solar  spectral  model. 
Most  of  the  essential  factors  are  listed  in  Figure  I.A.  The 
multiplicity  of  factors  and  their  interaction  make  it 
difficult  to  maintain  a  simple  yet  realistic  model  of  the 
effects  of  the  atmosphere.  However,  an  attempt  will  be  made 
to  include  all  of  the  factors  listed  in  the  figure  as  well 
as  others  when  needed. 

The  intent  behind  developing  this  model  is  to  construct 
a  simple,  solar  spectral  model  which  has  the  basic 
components  so  that  it  may  be  used  in  a  general  fashion  by  a 
large  group  of  users.  The  primary  user  groups  are  expected 
to  be  in  the  fields  of  solar  energy  engineering,  agriculture 
and  plant  ecology,  environmental  photochemistry,  and  remote 
sensing.  The  model  should  distinguish  between  the  global, 
direct,  and  diffuse  components  of  the  solar  energy  incident 
on  a  surface  of  arbitrary  orientation.  Hooper  and  Brunger 
(1980)  showed  that  the  diffuse  fraction  of  the  total  solar 
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Viewing  geometry: 

Solar  zenith  angle,  0 + 

Viewing  angle,  0r 
Azimuthal  angle,  or  <bT 
Relative  azimuthal  angle,  V'i 
where  ^  -  0,  +  180 

Altitude  of  sensor,  H 

Meteorological  parameters: 
Relative  humidity 
Cloud  cover 
Surface  pressure 


Atmospheric  optical  parameters: 
Optical  thickness,  rA 
Aerosol  type  (phase  function) 
Single-scattering  albedo,  u/0 

Target  and  background  parameters: 
Target  size 
Target  reflectance,  pt 
Background  reflectance,  pj, 
Instantaneous  field-of-view,  IFOV 


Figure  1 .A 

Factors  Affecting  Apparent  Reflectance  Determination 
(After  Bowker,  Davis,  Myrick,  Stacy  and  Jones,  1985) 
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radiation  reaching  the  surface  is  very  important,  especially 
for  the  more  northern  (southern)  latitudes  where  the  air 
mass  path  is  large.  They  also  pointed  out  that  local  climate 
(terrain,  proximity  to  bodies  of  water,  local  circulation) 
has  more  effect  than  latitude.  For  tilted  surfaces.  Hooper 
and  Brunger  (1980)  showed  for  11  Canadian  cities  that  about 
40  percent  of  the  annual  incident  short  wave  solar  radiation 
on  such  surfaces  is  diffuse  (sky  plus  ground  diffuse) 
radiation.  The  surface  albedo  is  an  important  factor  in  the 
radiative  interaction  between  the  atmosphere  and  the  ground. 
It  affects  the  amount  of  sky  and  ground  diffuse  radiation 
observed  at  a  particular  location. 

Many  models  account  for  surface  albedo  in  some  simple 
form,  usually  with  an  average  value  near  20  percent  for  the 
entire  spectrum.  However,  in  practice  the  surface  albedo  can 
vary  substantially  with  wavelength.  Figure  1.B  shows  the 
variation  in  reflectance  with  wavelength  for  typical 
vegetation  and  soil.  The  differences  are  rather  pronounced 
in  some  regions  of  the  spectrum.  The  effect  of  different 
spectral  albedos  is  not  usually  included  in  solar  radiation 
models.  The  reflectances  shown  in  Figure  1.B  are  hemispheric 
reflectances  measured  with  a  zenith  angle  relative  to  the 
surface  of  zero  degrees  (normal  to, the  surface).  The  surface 
reflectance  changes  with  the  solar  zenith  angle  and  the 
relative  orientation  of  the  surface.  These  factors  should  be 
included  in  any  general  solar  spectral  model. 


Figure  l.B 

Spectral  reflectance  curves  for  healthy  green  vegetation  and  air- 
dried  soils.  These  curves  represent  averages  of  240  spectra  from 
vegetation  and  154  spectra  from  air-dried  soils.  The  relative 
differences  in  reflectance  in  the  visible  (0.4  to  O.Tfim),  near- 
infrared  (0.7  to  1.3  jim),  and  middle-infrared  (1.3  to  3.0  fim) 
portions  of  the  spectrum  are  clearly  shown  in  this  data. 

(Adapted  from  Lillesand  and  Kiefer,  1979) 


6 


A  desirable  characteristic  of  a  general  solar  spectral 
model  would  be  to  allow  the  user  to  select  spectral  bands 
from  within  the  range  of  0.28  to  4.0  micrometers.  This  would 
enhance  the  usefulness  of  the  model  to  users  who  are  trying 
to  choose  particular  spectral  bands  for  their  application. 
Those  users  interested  in  remote  sensing,  agriculture,  or 
solar  energy  engineering  would  be  particularly  interested  in 
being  able  to  compare  the  radiation  available  in  different 
spectral  bands  on  slopes  of  varying  orientation  at  different 
solar  zenith  angles.  Another  factor  which  directly  affects 
the  amount  of  radiation  reaching  the  surface  is  the  cloud 
cover . 

Clouds  are  not  contained  in  many  of  the  solar  models 
mainly  due  to  the  complexity  of  the  calculations,  and  the 
difficulty  of  obtaining  the  observations  needed  to  verify 
the  model’s  results.  Clouds  have  been  included  so  the  user 
may  obtain  an  idea  of  how  clouds  may  affect  their  project. 
The  effect  of  clouds  is  governed  by  the  thickness  of  the 
cloud,  the  cloud  droplet  size  distribution,  the  size  and 
shape  of  the  clouds,  and  the  cloud  cover.  A  simple  spectral 
cloud  parameterization  has  been  included  in  the  model. 

In  summary,  a  general  spectral  solar  radiation  model 
has  been  constructed  which  has  the  following  features: 

1.  It  covers  the  solar  spectrum  between  0.28  and 
4.0  micrometers  with  126  finite,  non-uniform 
width,  irregularly  spaced  spectral  intervals. 

2.  It  computes  the  global,  direct,  and  diffuse 
radiation  incident  on  a  plane  surface  of 
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arbitrary  orientation  located  at  the  Earth’s 
surface.  The  direct-normal  radiation  at  the 
surface  is  also  computed.  Computations  of  the 
radiation  reaching  a  satellite  observation  point 
at  the  top  of  the  Earth's  atmosphere  are  made. 
This  radiation  is  computed  in  terms  of  path 
radiance,  global,  direct,  and  diffuse.  In  this 
case,  the  direct  radiation  is  the  radiation 
arriving  at  the  satellite  observation  point 
directly  from  the  target.  Diffuse  radiation  is 
radiation  from  the  target  which  is  scattered  to 
the  satellite  by  the  atmosphere.  Global 
radiation  is  the  sum  of  the  direct,  diffuse,  and 
path  radiance  values.  The  path  radiance  term  in 
this  model  xs  the  radiation  scattered  by  the 
atmosphere  to  the  observation  point.  The  path 
radiance  includes  no  radiation  from  the  target. 

3.  The  model  allows  the  user  to  select  any  spectral 
band  between  0.28  and  4.0  micrometers  for  the 
calculations. 

4.  The  model  uses  spectral  surface  albedos. 

5.  The  model  includes  the  dependence  of  albedo  on 
solar  zenith  angle. 

6.  A  simple,  yet  quantitative,  calculation  of  the 
spectral  effects  of  clouds  has  been  included  in 
the  model. 

7.  The  model  determines  the  spectral  solar 
radiation  reaching  an  arbitrarily  oriented 
surface  at  a  given  location  on  the  Earth. 

B.  The  model  adjusts  for  astronomical  factors  such 
as  the  Earth-sun  distance,  declination  of  the 
sun  and  sunrise/sunset  times. 

9.  The  model  accumulates  the  daily  solar  radiation 
on  a  surface  of  fixed  orientation. 

10.  The  model  allows  the  user  to  select  simple 
obstructions  to  sunrise/sunset  so  that  the 
effects  of  solar  radiation  in  mountainous  areas 
might  be  studied. 


The  model  includes  many  of  the  desirable  features 
of  a  general  solar  spectral  radiation  model.  The 
criteria  for  selecting  how  a  particular  component  or 
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feature  would  be  modelled  was  based  on  computer  time, 
how  accurately  the  form  reflected  physical  science,  and 
the  ease  of  inclusion  with  the  other  components  of  the 
model . 

The  unique  features  in  this  work  are: 

1.  the  inclusion  of  a  select  set  of  spectral 
surface  albedos 

2.  the  inclusion  of  spectral  cloud 
parameterizat ions 

3.  the  inclusion  of  inclined  surfaces 

4.  the  inclusion  of  simple  obstructions  to  sunrise 
and  sunset 

5.  the  ability  to  obtain  insolation  values 
representative  of  instantaneous,  hourly  or  daily 
time  intervals. 

The  program  is  written  in  ANSI  Standard  (1977)  FORTRAN.  It 
is  available  from  the  Department  of  Atmospheric  and  Oceanic 
Science,  The  University  of  Michigan. 


CHAPTER  II 
THE  MODEL 

In  this  chapter,  the  basic  components  of  the  model  will 
be  discussed.  In  some  instances,  the  model  allows  the  user 
the  option  of  choosing  which  technique  or  parameters  are  to 
be  used  in  the  model.  These  options  will  be  outlined  and 
discussed  where  necessary.  With  a  large  number  of  composite 
techniques  being  used  together  to  form  a  general  model, 
simplicity  and  accuracy  were  the  primary  criteria  for 
choosing  the  basis  of  the  model.  The  use  of  the  complex  and 
rigorous  SOLTRAN  and  LOWTRAN  radiative  codes  for  clear  skies 
seemed  computationally  unnecessary  as  we  are  striving  for  an 
overall  error  of  less  than  5%.  To  incorporate  clouds  into 
the  more  rigorous  radiative  transfer  schemes  could  lead  to 
computer  computation  times  much  higher  than  would  prove 
economically  feasible  to  many  users.  To  this  end,  the 
foundation  for  this  model  was  chosen  to  be  Bird's  simple 
solar  spectral  model  for  direct-normal  and  diffuse 
horizontal  irradiance  (1984).  Bird  compared  his  model  with 
the  results  of  the  BRITE  and  SOLTRAN  5  models  which  are 
believed  accurate  within  5%  on  the  direct-normal  and  within 
1 5%  on  the  diffuse  horizontal.  Bird's  model  compares  very 
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well  with  these  model  results.  Comparisons  of  this  model 
with  experimental  data  have  not  been  published.  The  Bird 
model  was  chosen  for  use  in  this  study  because  it  is  simple, 
accurate  and  has  good  spectral  resolution  (approximately  10 
nm).  Fodif ications  were  made  to  the  clear  sky  model  to 
account  for  inclined  surfaces. 

2. 1  Extraterrestrial  Solar  Spectrum 

The  Neckel  and  Labs  extraterrestrial  solar  spectrum  was 
used  by  Bird.  According  to  Bird,  Hulstrom,  and  Lewis  (1983), 
who  did  a  comparative  study  between  the  extraterrestrial 
data  sets  available,  the  Neckel  and  Labs  data  set  (1981)  is 
the  best  currently  in  use.  This  data  set  is  designed  to 
reflect  the  extraterrestrial  solar  spectrum  under  normal 
conditions  at  the  mean  Earth-sun  distance.  The  Neckel  and 
Lab's  spectrum  has  a  resolution  of  0.005  micrometers.  This 
is  better  resolution  than  the  spectrum  by  Thekaekara  which 
has  been  used  as  a  standard  by  many  investigators.  For  the 
purposes  of  developing  a  solar  spectral  model,  the  spectrum 
was  divided  into  126  non-uniform  width,  spectral  intervals 
between  0.28  and  4.0  micrometers.  With  the  exception  of  the 
first  four  values,  the  values  of  the  extraterrestrial  solar 
spectrum  are  those  derived  by  Bird  from  the  Labs  and 
Neckel' s  data  set  and  are  shown  in  Table  2.1. A.  The  first 
four  values  of  0.28,  0.285,  0.290,  and  0.295  micrometers 
were  derived  from  the  extraterrestrial  spectrum  of  Frolich 
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Table  2.1. A 


The  Neckel  and  Labs  revised  extraterrestrial 
and  atmospheric  absorption  coefficients  at  126 


spectrum 
wave  lengths . 


(Adapted  from  Bird,  1984) 


and  Werli  (1981)  of  the  World  Radiation  Center  in  Davos, 
Switzerland.  Bird's  original  clear  sky  model  did  not  extend 
to  0.28  micrometers  and  his  model  was  extended  to  include 
the  range  of  0.28  to  4.0  micrometers.  The  Werli  and  Frolich 
data  set  was  largely  based  on  Labs  and  Neckel's  revised 
extraterrestrial  solar  spectrum  (1981)  and  the  two  data  sets 
are  very  compatible.  Also  included  in  Table  2.1. A  are  the 
absorption  coefficients  for  water  vapor  (subscript  w),  ozone 
(subscript  o),  and  the  uniformly  mixed  gases  (subscript  u). 
The  units  of  the  absorption  coefficients  are  square 
centimeters  per  gram,  per  centimeter,  and  per  kilometer 
respectively. 

2.2  Astronomical  Factors 

Since  the  orbit  of  the  Earth  is  not  circular,  one 
should  adjust  the  values  of  the  extraterrestrial  solar 
spectrum  for  the  Earth-sun  distance  on  a  given  day  of  the 
year.  Presently,  the  Earth's  orbit  is  only  slightly 
elliptical  (the  eccentricity  is  small).  However,  the  square 
of  the  Earth-sun  distance  varies  by  +  3.5  percent;  hence, 
the  solar  irradiance  at  the  Earth  varies  by  ±  3.5  percent 
(Iqbal,  1963).  Paltridge  and  Platt  (1976)  give  an  expression 
for  the  square  of  the  ratio  of  the  mean  Earth-sun  distance 
to  the  Earth-sun  distance  with  the  dependent  variable  being 
the  day  of  the  year.  The  accuracy  of  this  expression  is 
quoted  as  being  better  than  0.0001.  The  expression  for  the 
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square  of  the  Earth-sun  ratio  is: 

(Rm  /  R)2  =  1.000110  +  0.034221  cos  (  01  +  0.001280  sin  (  0 

o  o 

*  0.000719  cos  (20}  +  0.000077  sin  (2 O'1 
o  o 

where 

Rm  =  mean  Earth-Sun  distance 

R  =  actual  Earth-Sun  distance 

60  =  an  expression  defined  in  terms  of  the  day  of  the  year 
2  Tdn  /  365 

dn  =  day  of  the  year  (0=Jan  1;  364=Dec  31) 

In  addition  to  the  correction  for  the  Earth-sun 
distance,  other  astronomical  factors  affect  the  solar 
radiation  arriving  at  the  Earth's  surface.  One  of  these 
factors  is  the  declination  of  the  sun.  The  declination  is 
the  angular  distance,  measured  along  an  hour  circle,  between 
the  sun  and  the  celestial  equator.  The  declination  varies  by 
day  of  the  year  and  can  be  expressed  in  the  following  form 
(Spencer,  1971;  Paltridge  and  Platt,  1976): 


6  =  0.006918  -  0.399912  cos 

-  0.006758  cos 

-  0.002697  cos 

where 


(  0) 

+  0.070257 

sin 

( 

0) 

0 

0 

(  2  0) 

+  0.000907 

sin 

( 

201 

o 

0 

(  30) 

+  0.001480 

sin 

( 

30) 

o  o 


6=  declination  in  radians 


The  expression  for  declination  has  a  maximum  error  of  0.0006 
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radians  or  an  error  less  than  3  minutes.  These  expressions 
can  be  found  in  Paltridge  and  Platt  (1976). 

Another  important  astronomical  factor  is  the  oblateness 
of  the  Earth.  Van  Hemelruck  (1983)  studied  the  effects  of 
oblateness  on  the  distribution  of  the  radiation  which 
reaches  the  Earth  and  the  other  planets  in  the  solar  system. 
Van  Hemelruck  found  the  oblateness  effect  on  the  solar 
radiation  is  a  few  tenths  of  a  percent.  The  largest 
differences  being  decreases  in  solar  radiation.  The  maximum 
increase  was  0.1  percent  while  the  largest  decrease  was  1.3 
percent  with  most  of  the  decreases  in  the  range  of  0.1  to 
0.5  percent.  The  study  showed  the  daily  insolation 
accumulated  at  a  given  location  could  be  decreased  by  2 
percent  due  to  oblateness.  The  net  result  of  oblateness  is 
to  decrease  the  amount  of  solar  radiation  reaching  the 
planet.  The  largest  losses  derived  from  this  effect  occur 
between  55  and  85  degrees  latitude  on  the  Earth.  Even  though 
the  effects  of  oblateness  could  have  been  included  in  the 
present  model,  it  was  neglected.  While  the  effects  of  the 
oblateness  can  be  determined  to  have  an  effect  of  a  few 
tenths  of  a  percent  on  the  solar  radiation,  the  effect  was 
not  significant  enough  to  be  included  in  this  study. 

2 . 3  Atmospheric  Transmittance 


The  atmospheric  transmittance  is  calculated  for  the 
regions  of  clear  sky,  i.e.  regions  without  cloud,  using  the 


formulation  by  Bird.  Bird's  model  compares  well  with  the 
BRITE  and  SOLTRAN  5  models  (See  Figures  2.4.F  and  G).  The 
BRITE  and  SOLTRAN  5  models  are  believed  to  be  accurate 
within  ±5%  on  the  direct-normal  radiation  and  within  ±  15% 
on  the  diffuse  radiation  (Bird,  1984).  Bird  uses  the 
multiplication  technique  for  obtaining  the  transmissions. 
The  basic  equation  used  for  the  direct-normal  component  of 
the  solar  radiation  is  : 


I  =  H  T  T  T  T  T 
dA  o  A  rA  «A  oA  uA  wA 

where 

I  =  Direct-normal  component  of  the  solar  radiation 

dA 


H  =  Extraterrestrial  solar  irradiance  at  the  mean  Earth-Sun 

oA 

distance 

T  =  Transmittance  function  for  Rayleigh  scattering 
rA 


T  =  Transmittance  function  for  aerosol  scattering 

•  A 


T  =  Transmittance  function  for  ozone  absorption 

oA 


T  =  Transmittance  function  for  water  vapor  absorption 

wA 


T  =  Transmittance  function  for  the  uniformly  mixed  gases 

uA 

(oxygen  and  carbon  dioxide) 
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Bird  and  Hulstrom  (1981)  studied  the  equations  used  for 
computing  the  direct  irradiance  from  several  different 
models.  Equations  of  similar  structure  have  been  referred  to 
in  the  literature  as  forms  of  the  transport  equation.  Bird 
and  Hulstrom  compared  several  forms  of  the  transport 
equation  developed  by  different  researchers  (Hoyt,  Watt, 
Lacis  and  Hansen,  Atwater  and  Ball,  and  ASHRAE),  they 
concluded  this  form  of  the  transport  equation  provided  the 
closest  agreement  with  the  more  rigorous  SOLTRAN  model. 

Based  on  this  finding,  the  Bird  clear  sky  model,  which 
already  contains  this  formalism,  was  used  as  the  basis  for 
this  solar  spectral  model.  An  inherent  assumption  has  been 
made  when  this  form  of  the  transport  equation  is  used.  The 
assumption  is  that  the  attenuation  by  each  constituent  is 
independent  of  every  other  constituent.  Restated,  this 
implies  that  the  transmittance  measured  through  pure 
materials  can  be  combined  in  the  form  of  the  transport 
equation  to  produce  the  transmittance  through  a  mixture  of 
the  materials.  This  is  not  true.  With  this  methodology,  two 
constituents  which  absorb  in  the  same  spectral  region  may 
cause  over-attenuation  or  under-attentuation  in  the  final 
results.  For  constituents  which  absorb  in  the  same  spectral 
region  in  different  layers  of  the  atmosphere,  this  form  of 
the  transport  equation  usually  results  in  over-attenuation 
in  the  final  results. 
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2.3.1  Scattering  Processes 

The  scattering  processes  discussed  throughout  this 
paper  have  some  underlying  assumptions.  All  the  scattering 
processes  have  been  assumed  to  be  independent  scattering 
which  occurs  when  the  average  separation  between  the 
particles  is  several  times  the  particle  radius  (McCartney, 
1976).  This  insures  that  the  scattering  pattern  from  one 
particle  is  not  affected  by  the  scattering  pattern  from 
neighboring  particles.  The  separation  criterion  is  satisfied 
in  all  typical  meteorological  conditions.  Since  the 
particles  in  the  atmosphere  are  randomly  arranged  and 
randomly  moving,  there  are  no  coherent  phase  relationships 
between  the  separately  scattered  radiation.  This  implies 
that  there  are  no  discernible  interference  patterns.  The 
resultant  simplification  is  that  the  intensities  of  the 
scattered  waves  rather  than  their  amplitudes  are  additive. 
This  type  of  scattering  is  called  incoherent  scattering. 

It  was  assumed  that  the  particles  are  exposed  only  to 
the  light  of  an  incident  or  direct  beam.  This  assumption  is 
referred  to  as  single  scattering.  Hence,  there  will  be 
references  to  single  scattering  ratios  and  cross  sections. 
"No  account  has  been  taken  of  the  fact  that  each  particle  in 
a  scattering  volume  is  exposed  to  and  also  scatters  a  small 
amount  of  the  light  already  scattered  by  the  other 
particles.  Processes  which  account  for  the  rescattered  light 
are  referred  to  as  multiple  or  secondary  scattering 


18 


processes."  (McCartney,  1976).  A  form  of  multiple  scattering 
has  been  included  in  the  model,  however.  The  multiple 
reflections  between  layers  of  atmosphere  is  a  form  of 
multiple  scattering.  The  adding  technique  is  used  to  combine 
the  atmospheric  layers  in  this  model  and  is  a  series 
approximation  to  account  for  the  multiple  reflections 
between  the  layers.  In  this  sense,  some  multiple  scattering 
is  accounted  for  in  this  model.  The  adding  method  was 
developed  without  regard  for  the  thickness  of  the  layers  to 
be  combined.  The  only  requirement  is  that  the  total 
transmittance  and  total  reflection  for  each  layer  be  known. 
Since  the  single  scattering  approximation  is  used  in  the 
model,  each  initial  layer  should  be  thin  enough  that  the 
single  scattering  approximation  for  the  reflection  and 
transmission  be  valid. 

Scattering  processes  are  dependent  upon  the  ratio  of 
the  particle  size  to  the  wavelength  of  the  light.  The  two 
primary  scattering  solutions  are  Rayleigh  and  Mie  scattering 
and  the  model  distinguishes  between  the  two  cases.  The 
effects  of  polarization  have  been  neglected  since  Hansen 
(1971)  has  shown  that  the  values  computed  from  exact  theory 
including  polarization  differ  by  less  than  0.1  percent  from 
the  scalar  approximation  in  which  polarization  is  neglected. 

Hansen  (1969)  demonstrated  that  it  was  not  necessary  to 
use  the  complete  cloud  phase  function  to  determine  the 
results  of  the  scattering  processes.  He  showed  that  the 
truncated  phase  functions,  while  giving  errors  in  the 
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angular  distribution  for  scattering  were  accurate  for  most 
practical  applications.  In  Hansen's  model  for  multiple 
scattering  in  cloudy  and  hazy  atmospheres,  he  performs  his 
calculations  using  only  the  forward  to  backward  scattering 
ratio,  a  procedure  that  is  adopted  here. 

2.3.2  Rayleigh  Scattering 

Bird  used  an  expression  adapted  from  the  LOWTRAN  5 
atmospheric  transmittance  model  to  compute  the  transmittance 
due  to  Rayleigh  scattering.  The  results  agree  well  with 
Pendorf’s  data  on  the  refractive  indices  of  standard  air  and 
the  Rayleigh  scattering  coefficients.  The  expression  for  the 
Rayleigh  scattering  transmittance  is  : 

TrA  =  exp  (  -Hy/[  A  4  (115.6406  -  1.335/ A  2)j) 

where 

H/  =  Pressure  corrected  air  mass 
\  =  Wavelength  in  micrometers 

To  obtain  the  pressure  corrected  air  mass,  one  must 
first  obtain  the  relative  air  mass.  Bird  elected  to  use  the 
formulation  for  the  relative  air  mass  derived  by  Kasten. 
Hasten' s  expression  for  the  relative  air  mass  is: 
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-1  2^3 

M  =  [  cos(Z)  +  0.15  (93.885  -  Z)  X  ] 

where 

M  =  Relative  air  mass  in  atmospheres 
Z  =  Apparent  solar  zenith  angle  in  degrees 


To  obtain  the  pressure  corrected  air  mass,  multiply  the 
relative  air  mass,  M,  by  the  ratio  P  /  P  . 


4. 


M'  = 


P 


M  ( - ) 


P 

o 


where  P  =  surface  pressure  in  millibars  (mb) 

P  =  1013  mb 
o 


2.3.3  Aerosol  Scattering  And  Absorption 

Aerosols  are  solid  or  liquid  particles  that  remain 
suspended  in  the  air  and  tend  to  follow  the  motion  of  the 
air.  The  typical  meteorological  precipitates  (rain,  snow  and 
hail)  are  not  aerosols.  Natural  aerosol  particles  range  in 

-3  2 

size  from  10  to  10  micrometers  (Iqbal,  1983).  The  sources 
of  these  particles  include  smoke,  haze,  clouds,  pollen, 
pollution,  volcanic  eruptions,  sandstorms,  forest  fires, 
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agricultural  burning,  meteoric  dust,  and  marine  based 
particulates  to  include  salt  crystals  and  ocean  spray.  The 
degree  of  turbidity  is  based  on  the  number  of  aerosol 
particles  per  cubic  centimeter  and  the  distribution  of  the 
particle  sizes.  One  of  the  biggest  problems  in  calculating 
solar  irradiance  at  the  surface  is  accounting  for  the 
effects  of  aerosols.  Part  of  the  difficulty  is  specifying 
the  properties  associated  with  the  aerosols.  The  properties 
include  their  size  distribution,  number,  and  optical 
characteristics . 

The  aerosols  produce  Mie  scattering  which  is  dependent 
on  the  ratio  of  the  particle  size  to  the  wavelength  of  the 
incident  radiation.  The  amount  of  scattering  will  also  be 
dependent  on  whether  the  dust  particles  are  wet  or  dry.  The 
optical  properties  of  the  particles  can  change  with  the 
degree  of  'wetness'.  The  more  detailed  treatments  of 
turbidity  consider  whether  the  particles  are  wet  or  dry,  as 
well  as  their  variability  in  form,  size,  distribution  and 
nature  of  the  particles.  Techniques  for  handling  aerosols 
can  be  reviewed  in  van  de  Hulst  (1957)  and  Iqbal  (1983).  A 
simple  technique  for  quantifying  the  effects  of  aerosols  was 
developed  by  Angstrom. 

Angstrom's  technique  is  not  only  simple  in  form,  but 
the  attenuation  effects  (whether  absorption  or  scattering) 
are  included  simultaneously.  In  practice,  it  is  difficult  to 
separate  the  attenuation  effects  of  scattering  and 
absorption  by  dust.  This  technique  includes  the  effects  of 
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both  wet  and  dry  particles;  i.e.  all  aerosols.  Angstrom 
assumed  a  dependence  of  the  attenuation  coefficient  of  the 
form: 


K  = 
a 


0  A  '  a 


where 


K  =  attenuation  coefficient 
a 


/}  =  turbidity  coefficient 
a  =  turbidity  exponents 


The  turbidity  coefficient  is  an  index  representing  the 
amount  of  aerosols  present  while  the  turbidity  exponent  is 
related  to  the  size  distribution  of  the  aerosol  particles.  A 
more  detailed  discussion  of  this  technique  can  be  found  in 
other  sources  (Angstrom  1961,  1964).  The  turbidity 
coefficient  and  exponent  can  be  determined  simultaneously 
with  a  dual  wavelength  sun  photometer.  The  wavelengths 
usually  chosen  for  the  sun  photometer  are  0.38,  0.5,  and 
0.88  micrometers.  The  model  uses  the  value  of  turbidity  at 
0.5  micrometers  as  its  definition  of  turbidity. 

In  the  more  rigorous  atmospheric  transmission  models, 
both  the  scattering  and  extinction  by  aerosols  is  computed 
using  Mie  scattering  theory.  Calculations  of  this  sort 
require  a  knowledge  of  the  aerosol  particle  size 
distribution,  the  complex  index  of  refraction  with 
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wavelength,  and  the  scattering  phase  function.  These 
variables  are  not  readily  available  for  each  individual 
scenario.  As  a  consequence,  a  typical  aerosol  distribution 
is  usually  chosen  for  study  purposes.  Bird  chose  to  use  the 
rural  aerosol  distribution  adopted  by  Shettle  and  Fenn 
(1975)  with  the  particle  size  distribution  shown  in  Figure 
2. 3. 3. A.  Two  aerosol  distributions  are  combined  to  form  the 
rural  aerosol  model.  To  describe  this  distribution,  Leckner 
(1978)  used  an  approximate  expression  of  the  Angstrom 
formalism  to  determine  the  turbidity  or  aerosol  optical 
depth  as  a  function  of  wavelength.  This  technique  assumes 
that  the  plot  of  turbidity  versus  wavelength  on  a  log-log 
scale  will  be  linear.  The  rural  aerosol  distribution  in  this 
model  shows  curvature  when  plotted  on  a  log-log  scale.  To 
compensate  for  this.  Bird  used  the  multiterm  Angstrom 
formalism  as  follows: 


—  a 

n 

T  =  A  \  n=l  for  0.28  -  -1  -  0.5  un 

a.\  n 

n=2  for  0.50  *  -  U.O  um 


where 


t  -  turbidity 

•  A 


fi  =  turbidity  coefficients  directly  proportional 


to  the  turbidity 


=  turbidity  exponent  related  to  the  aerosol  size 


distribut ion 
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distributions  which  combine  to  make  up  the  rural  model 

Figure  2. 3. 3. A 

Aerosol  particle  size  distribution  for  rural  aerosol  model 

(After  Bird,  1984) 


Tf  <0.5  un) 

S,  . 

£2 

0.1 

0.0490 

".0433 

0.27 

0.1324 

0.1170 

0.37 

0.1814 

o 

s£> 

o 

0.51 

0.2501 

0.2210 

Table  2.3.3.B 

Angstrom  turbidity  coefficients  for  several  turbidities 
in  the  rural  aerosol  model  (After  Bird,  1984) . 


25 


The  exponents  are  determined  by  the  slope  of  the  line  on  a 
log-log  plot  of  turbidity  versus  wavelength.  To  approximate 
the  rural  distribution  shown  in  Figure  2. 3. 3. A,  two  straight 
line  segments  were  used.  The  values  of  the  turbidity 
coefficient  and  the  turbidity  exponent  have  been  determined 
based  on  these  two  line  segments.  Table  2.3.3.B  shows  the 
Angstrom  turbidity  coefficients  for  several  turbidities 
using  the  rural  aerosol  model.  The  turbidity  coefficients  in 
the  table  can  be  used  to  generate  coefficients  for  other 
turbidities  by  the  expression: 


f> 


n 


) 


T 


O 


where 


=  new  turbidity  coefficient 

n 


P>o  =  old  turbidity  coefficient 


Tn =  new  turbidity  value 


=  old  turbidity  value 

O  J 

The  turbidity  coefficients  and  turbidity  exponents  are  valid 
for  two  different  regions  of  the  rural  aerosol  distribution, 
each  represented  by  a  straight  line.  The  subscripts  1  and  2 
represent  the  two  line  segments  or  regions  of  the  rural 
aerosol  model.  The  coefficient  and  exponent  with  the 
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subscript  1  are  valid  between  0.28  and  0.5  micrometers.  The 
coefficient  and  exponent  with  the  subscript  2  are  valid 
between  0.5  and  4.0  micrometers.  The  values  of  the  turbidity 
exponents  are: 


a  =  1.0274 

a  =  1.2060 

2 

The  turbidity  at  0.5  micrometers  is  the  input  to  calculate 
the  turbidity  coefficients  and  exponents  by  wavelength.  The 
transmittance  function  for  the  aerosols  is  given  by: 

—  a 

n 

T  =  exp  [  -  P>\  M]  (  for  n  =  1  or  2) 

a  A  n 


where 

M  =  Relative  air  mass  in  atmosphere 

2.3.4  Water  Vapor  Absorption 

There  are  several  atmospheric  constituents  which  affect 
the  absorption  of  the  incoming  solar  radiation.  Figure 
2. 3. 4. A  shows  the  regions  of  the  solar  spectrum  most 
affected  by  absorption.  The  four  primary  absorbing 
constituents  are  water  vapor  (H20),  ozone  (03 ) ,  oxygen  (0 2 ) , 
and  carbon  dioxide  (CC>2).  The  last  two  constituents  are 
often  called  the  uniformly  mixed  gases  because  their  mixing 
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Figure  2. 3. 4. A 

Solar  spectral  irradiance  outside  the  atmosphere  and  at  the 
surface,  for  a  solar  Eenith  angle  of  q° •  Features  due  to  principal 
absorbers  are  identified. 

(From  Bowker,  Davis,  Myrick,  Stacy  and  Jones,  1985) 
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ratio  is  nearly  constant  to  a  height  of  ICO  kilometers.  Bird 
uses  expressions  developed  by  Leckner  to  approximate  the 
influence  of  the  absorbing  gases.  In  each  case,  one  must 
know  the  integrated  amount  of  each  gas  in  a  vertical  column. 
Since  the  calculations  are  based  on  a  vertical  column,  the 
accuracy  of  the  calculated  absorption  will  decrease  with 
increasing  path  distance  through  the  atmosphere  (or 
increasing  zenith  angle).  The  discrepancies  are  usually 
attributed  to  the  variation  in  the  integrated  vertical 
amount  of  each  gas  with  distance  from  the  initial  region 
(Leckner,  1978).  Usually  the  deviations  from  the  calculated 
values  are  small  with  oxygen  and  carbon  dioxide.  Water  vapor 
has  the  most  variable  effect  on  the  absorption  of  the 
incoming  solar  radiation. 

Bird  used  a  modified  expression  of  Leckner 's  in  his 
model  to  compute  the  absorption  by  water  vapor.  The 
transmittance  function  for  water  vapor  used  in  Bird's  model 
was: 

TwA  =  exp  (  -0.3285  awA[w  +  (1.42  -  w)0.5]  M 


where 


0.45 

/  (1.0+  20.07  a  .  M) 

w  A 


) 


a^  ^  =  water  vapor  absorption  coefficient 

w  =  precipitable  water  in  a  vertical  path  in  cm 


2.3.5  Ozone  Absorption 
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The  ozone  transmittance  equation  of  Leckner  was  used  in 
Bird's  model.  The  expression  follows: 


where 


ToA  =  6XP  C 


-a  0- 
oA  3 


M  1 


a  .  =  Ozone  absorption  coefficient 
oA 

(See  Table  2.1 .A) 

0^  =  Ozone  amount  in  a  vertical  path  cm  (NTP) 

M  =  Air  mass  expression  for  ozone 


The  air  mass  expression  for  ozone  is  given  by  Paltridge  and 
Platt  as  : 


M 

o 


35.0  /  [1224.0  cos2  (z) 


->0  •  5 

J 


where 


z  =  zenith  angle  of  the  sun. 


If  the  ozone  amount  or  ozone  data  are  not  available,  the 
model  by  Van  Heuklon  can  be  used  to  determine  the  ozone 
amount  for  a  vertical  column.  The  model  by  T.K.  Van  Heuklon 
is  a  global  model  for  calculating  the  amount  of  ozone  in  a 
vertical  column  at  a  particular  location.  The  information 
required  to  determine  an  ozone  amount  are  latitude, 
longitude,  and  day  of  the  year.  Van  Heuklon  used  the 
observed  distribution  of  ozone  in  a  correlative  type 
expression  which  calculates  the  ozone  content  of  the 
atmosphere  for  any  place  in  the  world.  The  expression  for 
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the  ozone  content  is: 


where 


03  =  J  +  (  A  +  C  sin  [  D(  E+F )  ]  +■ 

G  sin[  H(A+  I)]  )  sin^ 


<p  =  Latitude  in  degrees 
A  =  longitude  in  degrees 


The  other  parameters  (A  thru  I)  are  shown  in  Table 
2. 3. 5. A.  The  values  of  the  parameters  vary  depending  on  the 
hemisphere.  Parameter  E  is  the  day  of  the  year  and  has  the 
same  formulation  regardless  of  hemisphere.  Deviations  in  the 
daily  ozone  amount  as  high  as  thirty  percent  can  occur 
(Hammond  and  Maugh,  1974).  This  model  primarily  predicts  the 
influence  of  the  seasonal  effects  of  the  general  circulation 
on  the  ozone  amount  rather  than  the  effects  of  diurnally 
controlled  forcing  functions.  This  empirical  model  is 
accurate  within  5  percent  for  most  areas  of  the  world.  The 
model  values  typically  range  from  a  1  percent 
underestimation  to  a  5  percent  overestimation.  In  the  higher 
latitudes  (above  80  degrees),  the  model  is  believed  to  be 
fairly  accurate.  However,  the  lack  of  observed  values  for 
comparison  prohibit  a  statement  about  the  accuracy  of  the 
model  for  latitudes  greater  than  60  degrees.  Figure  2.3.5.B 
shows  of  the  deviation  between  the  calculated  and  observed 
annual  averages  for  North  America  (Van  Heuklon,  1979). 
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Parameter 

N.  hemisphere 

Value  of  parameter 
both  hemispheres  S.  hemisphere 

A 

150  0 

_ 

1000 

0 

1.18 

— 

1.5 

c 

40.0 

— 

30.0 

D 

0.9865 

E 

Jan.  1  =  1.0.  2  =  2.0. 

etc. 

F 

-30.0 

— 

152.625 

G 

:o.o 

H 

3.0 

— 

2.0 

1 

:u.u  if  a  =  + 

— 

-75.0 

0.0  if  A  =  - 

J 

235.0 

6 

N  =  + 

S  =  - 

A 

E  =  +,  W  =  - 

Table 

2. 3. 5. A 

Values 

of  parameters  in  the  ozone 

mode  1 

(After  Van 

Heuklon,  1979) 

Comparison  calculated  vs  observed  annual  averages  of  atmospheric 

ozone  for  North  America 


(After  Van  Heuklon,  1979) 


The  solar  spectral  model  developed  by  this  author 
contains  an  option  where  the  user  may  elect  to  use  Van 


Heuklon’s  model  for  estimating  the  ozone  amount  or  enter  the 
ozone  amount  of  his  choice  in  atm-cm. 

2.3.6  Absorption  By  The  Uniformly  Mixed  Gases 

Bird  used  Leckner's  expression  to  compute  the 
transmittance  function  for  the  uniformly  mixed  gases.  The 
uniformly  mixed  gases  in  this  model  are  oxygen  and  carbon 
dioxide.  The  transmittance  function  for  the  uniformly  mixed 
gases  was: 


T  =  exp  [  -1.41  a  M' /  (l  +  118.93  a  ,  m')0-45  ] 
uA  u.v  uA 

where 

a^  =  Absorption  coefficients  for  the  uniformly  mixed 
gases  (See  Table  2.1. A) 

M'  =  Pressure  corrected  air  mass. 

2.4  Diffuse  Horizontal  I rradiance 

The  diffuse  irradiance  is  hard  to  simply  quantify.  The 
distribution  of  the  diffuse  component  is  non-uniform  over 
the  sky  and  changes  with  the  zenith  angle  of  the  sun.  Brine 
and  Iqbal  (1983)  developed  a  method  of  computing  the  diffuse 
radiation  based  on  the  broadband  work  of  Davies  and  Hays. 
Bird  used  the  equations  of  Brine  and  Iqbal  except  that  he 


DIFFUSE  REFLECTED  BACK  TO  SPACE 


Partition  of  radiation  into  direct,  diffuse  propagated  forward, 
and  diffuse  reflected  back  to  space. 

(From  Iqbal,  1983) 
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added  a  correction  factor.  Figure  2. 4. A  from  Iqbal  (1983) 
shows  the  direct  radiation,  total  scattered  radiation,  the 
diffuse  radiation  propagated  forward  and  the  diffuse 
radiation  relected  back  to  space.  The  equations  used  by  Bird 
to  express  the  diffuse  radiation  make  use  of  the  radiation 
components  shown  in  the  figure  to  compute  the  diffuse 
radiation  on  a  horizontal  surface.  The  equations  for  the 
diffuse  components  of  radiation  are: 


I 


S.l 


a  A 


)  C 


A 


I 
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H  cos( Z) 
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where 


I 


sA 


Scattered  components  on  a  horizontal  surface 


g A 


rA 


a  A 


oA 


A 


Ground  /  atmosphere  irradiance  on  a 
horizontal  surface  at  wavelength,  A . 

Rayleigh  scattered  component  on  a  horizontal  surface  at 
wavelength , A . 

Aerosol  scattered  component  on  a  horizontal  surface  at 
wavelength,  A . 

Solar  irradiance  at  the  top  of  the  atmosphere  corrected 
for  Earth-sun  distance 
Single  scattering  albedo 

Forward  to  total  scattering  ratio  of  the  aerosol 
Correction  factor  by  wavelength 
Albedo  of  the  air  (sky) 

Albedo  of  the  ground 


The  primed  transmittance  functions  indicate  the  values  are 
calculated  using  an  air  mass  value  of  1,9.  A  value  of  1.9 
air  masses  produces  the  closest  agreement  for  the  diffuse 
irradiance  when  compared  to  the  BRITE  model.  Bird  used  a 
value  of  0.928  for  the  single  scattering  albedo  and  a  value 
of  0.82  for  the  forward  to  total  scattering  ratio.  The  model 
defaults  to  these  values  for  the  single  scattering  albedo, 


forward  to  total  scattering  ratio,  and  the  air  mass  for  the 
primed  transmittance  quantities  unless  the  user  elects  to 
enter  his  own  values  for  these  quantities  or  use  the 
spectral  values  of  these  parameters  coded  into  the  model. 

A  correction  factor  was  needed  to  fit  the  equations  to 
Bird's  data  for  several  turbidity  and  zenith  angle 
combinations.  Bird  found  it  impossible  to  take  the  equations 
of  this  general  form  and  get  them  to  accurately  represent  a 
variety  of  conditions.  These  simple  equations  exhibit 
several  problems: 

1.  the  cos(Z)  terms  apply  only  to  direct  irradiance 

2.  the  transmittance  terms  apply  to  direct 
irradiance  only  and  are  not  expected  to  apply  to 
scattered  irradiance 

3.  W  and  F  are  wavelength  dependent 

o  a 

Attempts  by  Bird  to  use  a  variable  metric  minimization 
method  to  allow  modified  forms  of  the  equations  were 
unsuccessful.  The  correction  factor  was  computed  for  several 
wavelengths  as  a  function  of  turbidity  and  zenith  angle. 

Bird  assumed  that  linear  interpolation  could  be  used  between 
the  tabulated  wavelengths  to  obtain  the  proper  correction 
factors.  Tables  2.4.B  thru  E  show  the  correction  factor  as  a 
function  of  wavelength,  turbidity,  and  zenith  angle.  Bird 
assumed  a  linear  interpolation  scheme  could  be  used  to 
compute  the  correction  factors.  Figures  2.4.F  and  G  compare 
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Table  2.4.B 

Diffuse  correction  factor  forra(0.5)  =  0.1 
(From  Bird,  1984) 


0.30  0.66  0.93  1.02  1.23  2.00  9.00  6.3 

0.33  1.00  1.07  1.11  l,l»  1. 31  1.97  3. 70 

0.90  1.11  1.13  1.1)  1.29  1.96  1.70  2.61 

0.93  1.09  1.03  1.09  l.ll  1.29  1.39  1.72 

0.30  1.15  1.00  1.00  0.99  1.06  1.07  1.22 

0.33  1.12  0.96  0.96  0.99  0.99  0.96  1.09 

0.71  1.32  1.12  1.07  1.02  1.10  0.90  O.SO 

0.76  1.23  1.09  1.03  1.07  1.00  0.85  0.78 

0.9935  0.98  0.99  0.91  0.83  0.90  0.79  0.79 

2.1  0.69  0.66  0.71  0.69  0.70  0.59  0.32 

9  1  0.69  0.68  0.71  0.69  0.70  0.59  0.52 


Table  2.4.C 

Diffuse  correction  factor  for  ta(0.5)=0.27 
(From  Bird,  1984) 
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Bird's  clear  sky  spectral  model  (SPECTRAL)  against  the  more 
rigorous  BRITE  model.  Because  of  the  close  comparison, 

Bird's  clear  sky  is  expected  to  be  within  the  same  accuracy 
range  as  the  BRITE  model. 

In  the  model  developed  for  this  research,  several 
changes  were  made  which  happen  to  coincide  with  changes  Bird 
made  in  his  most  recent  model  (Bird  and  Riordan,  1984).  The 
single  scattering  albedo  (W^)  and  the  forward  to  total 
scattering  ratio  (F  )  have  been  made  spectrally  dependent 

d 

based  on  the  rural  aerosol  distribution.  Plots  of  the 
spectral  values  of  the  single  scattering  albedo  and  the 
forward  to  total  scattering  ratio  for  a  clear  sky  based  on 
the  rural  aerosol  distribution  used  by  Bird  are  shown  in 
Figure  2.4.1.  The  spectral  values  are  interpolated  from  data 
provided  by  Bird  (1984).  After  reviewing  Bird's  data  and  the 
information  in  Kerker  (1969)  concerning  scattering  ratios, 
it  was  demonstrated  that  the  parameters  varied  slowly  enough 
so  that  linear  interpolation  of  the  single  scattering  albedo 
and  the  forward  to  total  scattering  ratio  could  be  used  over 
the  wavelength  range  of  interest.  If  the  user  desires,  the 
model  gives  the  option  of  choosing  constant  values  for  both 
of  the  parameters.  The  model  also  allows  the  choice  of 
whether  to  use  Bird's  correction  factor.  These  two  options 
give  the  user  the  flexibility  to  compare  results  using 
slightly  different  techniques.  With  the  correction  factor, 
the  equations  represent  Bird’s  clear  sky  model  (1984). 
Without  the  correction  factor,  the  equations  represent  Brine 


Irradiance  (W/mV^im) 


Figure  2.4.F 

Comparison  of  SPECTRAL  and  BRITE  spectra  for  USS  atmosphere 
t  (0.5)=0.1,  ALB=0 .0  and  AMI. 


(From  Bird,  1984) 


Irradiance  (W/mV(jm| 
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Figure  2.4.G 

Comparison  of  SPECTRAL  and  BRITE  spectra  for  USS  atmosphere, 

r  (0.5) =0.1,  ALB=0.0  and  AM  5.6 
a 

(From  Bird,  1984) 
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WZERO  AND  FSUBA 


Figure  2.4.H 

Spectral  single  scattering  albedo  (W^)  and  the  spectral  forward 
to  total  scattering  ratio  (F  )  for  a  clear  sky  based  on  the  rural 
aerosol  distribution. 


and  Iqbal's  clear  sky  formulation  (1983) 


2.5  Model  Geometry 

Most  of  the  geometrical  quantities  needed  for  the  model 
were  calculated  from  equations  developed  by  Enrico  Coffari 
which  were  published  in  Solar  Energy  Engineering  edited  by 
A.A.G.  Sayigh  (1977).  While  the  coordinate  system  used  in 
somewhat  unconventional,  it  has  the  advantage  of  producing 
an  analytic  set  of  equations  which  minimizes  the  number  of 
computational  checks  needed  to  verify  whether  an  angle  is  in 
a  particular  quadrant  or  it  is  in  front  or  behind  a  plane. 
The  equations  are  based  on  the  use  of  vectors  and  the 
direction  cosines  of  the  vectors  combined  with  straight 
forward  principles  of  solid  geometry.  The  equations  allow 
the  computation  of  the  following  quantities: 

(1)  Solar  quantities 

a.  Solar  altitude 

b.  Solar  azimuth 

c.  Sunrise  time 

d.  Sunrise  azimuth 

e.  Sunset  time 

f.  Sunset  azimuth 

(2)  Angle  of  incidence  of  direct  radiation  on  an 
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inclined  plane 

(3)  Zero-irradiation  solar  times  (Planerise  and 
planeset  times;  or  times  when  the  sun  rises  and 
sets  on  an  inclined  plane) 

a.  Solar  altitude  and  azimuth  at  planerise 

b.  Solar  altitude  and  azimuth  at  planeset 

A  list  of  the  equations  used  to  compute  these  quantities  is 
contained  in  Appendix  C.  To  calculate  the  amount  of 
radiation  arriving  on  an  inclined  plane  surface  one  must  be 
able  to  determine  whether: 

a.  Only  diffuse  radiation  is  reaching  the 
inclined  plane 

b.  Both  diffuse  and  direct  radiation  are 
reaching  the  inclined  plane 


A  comparison  of  the  sunrise/sunset  times  with  the 
planerise/planeset  times  allows  the  determination  of  whether 
any  radiation,  diffuse  radiation  only,  or  both  diffuse  and 
direct  radiation  are  reaching  the  inclined  plane  surface.  In 
the  more  northern  and  southern  latitudes  it  is  possible  to 
have  an  inclined  plane  orientation  so  that  the  sun  will 
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begin  and  end  shining  on  the  inclined  plane  twice  in  one 
day.  In  the  polar  regions,  the  sun  may  not  rise  or  set 
during  a  twenty-four  hour  period;  the  sun  either  shines  all 
day  or  not  at  all.  Outside  the  polar  regions,  there  is 
always  one  sunrise/sunset  at  a  particular  location  on  the 
Earth.  The  orientation  of  the  inclined  plane  and  its 
location  on  the  Earth  determine  whether  one  or  two 
planerise/planeset  combinations  will  occur. 

The  angle  of  incidence  of  direct  radiation  on  a 
horizontal  surface  at  the  same  location  as  the  inclined 
plane  is  also  important.  The  radiation  reaching  the  area 
surrounding  the  inclined  plane  is  used  for  computing  the 
diffuse  component  of  the  radiation.  The  solar  zenith  angle 
relative  to  a  horizontal  surface  allows  the  model  to  compute 
the  radiation  from  the  surrounding  area  and  determine  the 
diffuse  radiation  reaching  the  inclined  plane.  The  ground 
diffuse  radiation  from  the  vicinty  is  assumed  to  be  from 
horizontal  or  flat  terrain. 

Both  the  angle  of  incidence  of  direct  radiation  on  a 
horizontal  surface  and  on  an  inclined  plane  affect  the 
albedo  of  their  respective  surfaces.  These  angles  are  used 
to  modify  the  surface  albedo. 

2.6  Surface  Albedo 


The  bidirectional  reflectance  varies  with  the  angle  of 
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incidence  and  the  angle  of  observation. 

In  this  model,  bidirectional  reflectance  is  constant. 

The  ground  (surface)  is  assumed  to  be  a  diffuse  or 
Lambertian  reflector. 

Paltridge  and  Platt  (1976)  developed  an  expression  for 
the  variation  of  albedo  with  the  angle  of  incidence.  Figure 
2. 6. A  shows  the  typical  variation  with  the  angle  of 
incidence  (elevation  angle)  of  the  radiation  source.  The 
figure  demonstates  how  a  change  in  solar  altitude  can  affect 
the  albedo  of  a  surface.  This  variation  has  been  included 
using  the  simple  form  developed  by  Paltridge  and  Platt.  The 
expression  used  by  Paltridge  and  Platt  to  compute  the  change 
in  albedo  with  angle  is: 

P„  (02  }  =  Pa  *  (I  ~  Pa)  expt  ~k(9°  ~  V  J 

-  Hemispheric  reflectance  of  the  surface  or  the  albedo 
of  the  surface  with  the  source  at  an  elevation  angle 
of  90  degrees  with  respect  to  the  surface. 

-  Albedo  or  reflectance  of  the  surface  as  a  function  of 
the  solar  zenith  angle 

-  A  constant  parameter  whose  value  is  of  order  0.1 

-  Solar  zenith  angle  relative  to  an  arbitrarily  oriented 
surface  (See  Figure  A.l  in  Appendix  A) 

A  value  of  k  equal  to  0.1  was  used  in  this  model.  It  is  the 
value  of  k  suggested  for  use  by  Paltridge  and  Platt. 


where 


P  (fl  ) 
g  2 


Figure  2. 6. A 

Angular  variation  of  the  albedo  for  various  ground 
covers  (From  Paltridge  and  Platt,  1976). 


In  addition  to  the  variation  of  the  albedo  with  angle 
of  incidence,  the  albedo  varies  by  wavelength.  This 
variation  is  often  neglected.  Most  models  assume  an  albedo 
which  is  constant  with  wavelength.  A  glance  at  some  of  the 
reflectance  curves  for  typical  materials  shows  the  degree  of 
spectral  variability.  Seven  typical  spectral  surface  albedos 
have  been  installed  in  the  model  and  the  user  selects  the 
surface  albedo  of  his  choice.  The  model  also  gives  the 
option  of  selecting  a  constant  spectral  surface  albedo  or 
the  option  of  entering  a  spectral  albedo  of  the  users 
choice.  In  summary,  the  model  allows  the  user  to  choose  a 
spectral  surface  albedo  from  the  following: 

1 .  Fresh  snow 

2.  Clay  (5%  moisture  by  weight) 

3.  Dry  bare  soil 

4.  Dry  sand 

5.  Wet  sand  (25*  moisture  by  weight) 

6.  Typical  vegetation 

7.  Turbid  water 

8.  A  constant  spectral  albedo  of  the  user's  choice 

9.  A  spectral  albedo  of  the  user’s  choice 

Figure  2.6.B  shows  plots  of  the  first  seven  categories  of 
hemispheric  reflectance  listed  immediately  above.  The  seven 
categories  were  chosen  to  represent  a  variety  of  surface 
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HEMISPHERIC  REFLECTANCE 


HEMISPHERIC  REFLECTANCE 


Legend 


A  ORY_SANO_  _  □  CLAY _ 

X  WET  SAND  B  TURBID  WATER 


Figure  2.6.B 

The  figures  show  the  seven  spectral  reflectance  curves  used 
in  the  model. 
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reflectances  to  give  the  user  the  opportunity  to  make  model 
calculations  and  see  the  effects  different  spectral  albedos 
make. 

A  water  surface  was  included  in  the  list  of  modelled 
albedos  even  though  it  violates  the  assumption  that  the 
surface  is  a  diffuse  (Lambertian)  reflector.  A  water  surface 
was  included  because  most  of  the  surface  of  the  Earth  is 
covered  by  water.  Water  surfaces  have  specular  properties 
which  means  they  tend  to  act  like  mirrors  at  certain  zenith 
angles.  The  albedo  of  water  surfaces  is  also  affected  by 
wind  and  wave  disturbances.  The  wind  and  waves  determine  the 
wave  structure  or  ’choppiness’  of  the  water  surface  which 
has  a  direct  effect  on  the  albedo  of  a  water  surface.  In 
addition,  the  degree  of  pollution  affects  the  observed 
albedo  of  a  water  surface.  For  these  reasons,  a  water 
surface  may  have  an  albedo  very  different  from  the  one 
included  with  this  model.  The  turbid  water  included  in  this 
model  has  a  low  albedo  and  is  assumed  to  be  under 
non-specular  conditions.  The  low  albedo  of  water  is  typical 
of  most  water  surfaces  in  the  natural  environment. 

The  effects  of  surface  albedo  are  complex.  There  are 
models  to  determine  the  albedo  of  a  particular  surface  or 
plant  canopy.  Due  to  the  influence  cf  other  parameters  which 
affect  the  albedo,  more  information  is  required  by  these 
models  than  the  solar  zenith  angle.  Some  of  these  additional 
factors  are: 


51 


(a)  Season 

(b)  Recent  rainfall 

(c)  Leaf  orientation  or  surface  structure 

(d)  Composite  albedo  (i.e.  bare  soil  showing 
through  a  plant  canopy) 

It  should  be  obvious  that  these  factors  are  interrelated. 

For  example,  the  season  will  determine  the  degree  of 
foliation  which  will  affect  the  amount  of  bare  soil  showing 
through  the  plant  canopy.  The  surface  structure  or  leaf 
orientation  may  be  dependent  on  the  wind,  type  of  plant, 
season,  or  rainfall.  The  amount  and  frequency  of  rainfall  is 
dependent  on  season.  The  degree  of  complexity  in  computing 
albedo  can  become  great  quickly. 

For  this  model,  the  surface  is  assumed  to  be  uniformly 
made  of  the  material  chosen  by  the  user.  In  other  words, 
there  are  no  assumptions  concerning  the  amount  of  bare  soil, 
grass,  etc.  that  could  be  showing  through  the  vegetation 
canopy.  The  surface  is  assumed  to  consist  of  a  material 
which  is  uniform  in  texture  and  has  the  spectral  reflectance 
of  the  chosen  category.  The  spectral  reflectance  is  modified 
at  each  wavelength  interval  by  the  effect  of  the  solar 
zenith  angle. 

Although  this  model  does  account  for  changes  in  the 
reflectance  of  the  surface  due  to  the  solar  zenith  angle,  it 
does  not  strictly  account  for  bidirectional  reflectance.  In 
practice,  the  observed  reflectance  depends  on  the  angle  from 
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which  the  material  is  viewed.  Since  every  material  has 
different  characteristics,  it  is  not  feasible  in  a  simple 
model  to  adequately  represent  the  anisotropic  scattering  of 
radiation  by  a  surface.  Most  surfaces  reflect  radiation  in 
an  anisotropic  sense  which  is  the  reason  one  should  account 
for  bidirectional  reflectance.  Since  the  surfaces  in  this 
model  are  assumed  to  be  Lambertian  (perfectly  diffuse 
reflectors),  this  means  that  the  radiation  reflected  off  the 
surface  is  isotropic  and  the  bidirectional  reflectance  is 
constant . 

2.7  Clouds 

Clouds  present  a  number  of  problems.  The  first  and 
foremost  of  these  problems  is  the  complexity  and  detail 
needed  to  do  an  accurate  job  of  computing  their  effects.  It 
would  be  unsuitable  to  include  all  the  aspects  of  clouds 
theorized  in  the  more  current  cloud  models.  Due  to  their 
complexity,  the  amount  of  computation  time  needed  for  some 
of  these  models  would  be  prohibitive.  A  simple,  accurate 
accounting  of  the  effects  of  cloud  cover  on  the  incoming  and 
outgoing  solar  radiation  is  desired.  Several  parameters 
important  to  cloud  modelling  are  listed  below: 

1.  Cloud  geometry 

a.  Shape  of  the  cloud 

b.  The  cloud’s  proximity  to  other  clouds 

c.  The  height  of  the  cloud  base 

d.  The  dimensions  of  the  cloud 
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2.  Cloud  droplet  size  distribution 

3.  Mie  scattering 

4.  Multiple  scattering 

5.  Atmospheric  water  vapor  profile  and  cloud  liquid 
water  content 

6.  Temperature  structure  in  the  cloud 


These  aspects  will  be  discussed  in  relation  to  the  cloud 
model  that  was  developed  for  this  solar  radiation  model. 

2.7.1  Cloud  Characteristics 

Finding  a  simple,  accurate  cloud  model,  proved  to  be  an 
impossible  task.  The  cloud  models  available  were  either  very 
complex  in  terms  of  the  radiative  transfer  scheme,  were 
broadband  versus  spectral  in  nature,  geometrically  complex, 
or  required  data  not  readily  available  to  most  users.  Since 
this  model  is  intended  to  be  spectral,  the  first  priority 
was  to  find  a  cloud  study  or  cloud  model  which  is  either 
spectral  in  nature  or  could  be  expressed  in  spectral  terms. 
The  cloud  characteristics  described  by  Welch,  Cox,  and  Davis 
(1980)  were  spectral.  The  solar  spectrum  was  divided  into 
eight  regions.  The  center  wavelengths  of  the  eight  regions 
were  0.765,  0.950,  1.15,  1.40,  1.85,  2.80,  3.35,  and  6.30 
micrometers.  The  study  used  a  variety  of  cloud  types.  Of  the 
cloud  types  listed  in  the  paper,  four  were  chosen  for  their 
diversity  in  covering  the  range  of  the  most  commonly 


reported  low  clouds.  The  low  cloud  types  chosen  were 
stratocumulus ,  stratus,  nimbostratus ,  and  cumulus. 

The  data  used  by  Welch,  Cox,  and  Davis  includes 
specific  details  about  the  clouds.  One  of  the  primary 
considerations  is  the  cloud  droplet  size  distribution.  The 
drops  are  assumed  to  be  liquid;  no  ice  is  contained  in  these 
clouds.  Ice  clouds  will  have  different  characteristics  from 
water  droplet  clouds.  Using  this  model  for  ice  clouds, 
without  modification,  will  result  in  large  errors.  Welch  and 
Cox  used  the  drop  size  distribution  functions  of  Best  (1951) 
and  Deirmendjian  (1975).  The  drop  size  distribution  varies 
with  the  height  in  the  cloud.  As  a  consequence,  Welch  and 
Cox  list  the  droplet  size  distributions  for  the  base  and  the 
top  of  the  clouds.  The  drop  size  distribution  for  the  four 
cloud  types  used  in  this  model  are  given  in  terms  of  the 
modified  gamma  size  distribution  function  developed  by 
Deirmendjian.  The  modified  gamma  distribution  function  has 
the  following  form; 

n(r)  =  a  ra  exp[-(ct/i)  ( r/r  )“] 

c  J 


where 

n(r)  -  Number  density  of  droplets  with  radius,  r,  expres¬ 
sed  in  number  by  cubic  centimeter  per  micrometer 
r  -  Droplet  radius  in  micrometers 
r 

c 


-  Modal  radius  of  the  distribution  in  micrometers 
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Table  2. 7.1. A  (From  Welch,  Cox,  and  Davis  1980) 


Cloud  Droplet  Size  Distribution  Parameters 


Type  a  a  8  r  W  N 

c  L 

(fxm)  (g/rn  )  (cm-3) 


Stratocumulus 

base 

2.8230 

5 

1.19 

5.33 

0.141 

100 

Stratocumulus 

top 

1.9779 

2 

2.46 

10.19 

0.796 

100 

Stratus 

base 

9.7923 

5 

1 .05 

4.70 

0.114 

100 

Stratus 

top 

3.8180 

3 

1.30 

6.75 

0.379 

100 

N imbostratus 
base 

8.0606 

5 

1.24 

6.41 

0.235 

100 

N imbostratus 
top 

1 .0969 

1 

2.41 

9.67 

1.034 

100 

Cumulus 

0.5481 

4 

1 .00 

6.00 

0.297 

100 

a, a, S  -  Empirically  derived  cloud  droplet  distribution 
constants 

r^_  -  Modal  radius  of  the  cloud  droplet  distribution 

W  -  Cloud  liquid  water  content 
Li 

N  -  Number  of  cloud  droplets  per  cubic  centimeter 
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The  other  parameters  are  empirically  derived  constants. 

Table  2. 7.1. A  gives  the  original  values  of  the  parameters 
used  by  Welch  and  Cox.  Figure  2.7.1.B  shows  the  shape  of  the 
distributions  used  by  Welch  and  Cox.  The  cumulus  cloud 
formulation  was  taken  from  Deirmendjian  and  has  just  one 
droplet  distribution  throughout  the  cloud,  rather  than  a 
different  droplet  distribution  at  the  top  and  base  of  the 
cloud. 

Other  characteristics  of  the  clouds  used  by  Welch  and 
Cox  which  should  be  noted.  The  clouds  are  assumed  to  have 
the  temperature  structure  of  the  standard  atmosphere.  The 
clouds  are  also  assumed  to  be  saturated  at  each  level  within 
the  cloud.  For  this  research,  the  user  is  given  the  option 
of  choosing  the  surface  temperature.  This  option  does  not 
directly  affect  the  modelling  of  the  clouds  or  the 
assumptions  concerning  their  use. 

I  have  defined  the  cloud  model  by  establishing  the 
temperature  of  the  cloud,  the  water  vapor  content,  and  the 
drop  size  distributions.  Next  we  must  consider  the  effects 
of  Mie  scattering  and  geometry. 

2.7.2  Mie  Scattering 

Mie  scattering  defines  a  scatter  process  which  relates 


the  intensity  of  the  scatter  to  the  size  parameter.  The  size 
parameter  is  proportional  to  the  ratio  of  the  drop  radius 
divided  by  the  wavelength  of  light.  When  Welch  and  Cox 
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(1980)  developed  the  cloud  droplet  extinction  and  absorption 
coefficients  for  the  clouds  modelled  in  their  study,  they 
used  a  wavelength  dependent  Mie  efficiency  factor  and  the 
pertinent  cloud  drop  size  distribution.  However,  these 
factors  were  computed  at  only  eight  wavelengths.  The  cloud 
droplet  extinction  coefficients  are  needed  at  126 
non-uniform  width,  irregularly  spaced  intervals  across  the 
region  of  0.28  to  4.0  micrometers.  Deirmendjian  showed  that 
the  droplet  absorption,  scattering,  and  extinction 
parameters  are  relatively  slowly  varying  functions  of 
wavelength.  The  droplet  extinction  and  absorption 
coefficients  for  three  different  types  of  clouds  have  been 
plotted  by  Deirmendjian.  Figure  2. 7. 2. A  shows  the  plot  of 
the  droplet  extinction  and  scattering  coefficients  for 
various  clouds.  Because  the  coefficients  are  slowly  varying, 
it  was  assumed  that  linear  interpolation  could  be  used  to 
obtain  the  values  of  the  droplet  extinction  and  absorption 
coefficients  between  the  wavelengths  where  Welch  and  Cox 
rigorously  computed  their  values.  Linear  interpolation  and 
extrapolation  was  used  to  determine  the  126  droplet 
extinction  and  absorption  coefficients. 

Welch  and  Cox  have  two  droplet  size  distributions  for 
some  of  the  clouds.  They  generated  different  droplet 
extinction  and  absorption  coefficients  at  the  base  and  the 
top  of  the  clouds.  For  simplicity,  it  was  desirable  to 
obtain  only  one  set  of  coefficients  to  represent  each  cloud. 
To  obtain  a  single  droplet  absorption  or  extinction 
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coefficient,  the  values  of  the  droplet  coefficients  at  the 
base  and  top  of  the  cloud  were  averaged.  It  was  assumed  the 
average  droplet  extinction  and  absorption  coefficients  would 
sufficiently  represent  the  entire  cloud.  Table  2.7.2.B  shows 
the  values  computed  by  Welch  and  Cox.  Table  2.7.2.C  shows 
the  averaged  values  used  for  this  model. 

Mie  scattering  affects  the  amount  of  radiation  forward 
scattered  and  backward  scattered.  Welch  and  Cox  (after 
Zdunkowski,  1967)  determined  phase  function  expansion 
coefficients  for  a  series  of  Lengendre  polynomials  which 
represent  the  proportion  of  total  energy  scattered  in  a 
particular  direction.  Assuming  single  scattering  and 
approximating  the  scattering  phase  function  by  a  Dirac  delta 
function  for  the  forward  scattering  peak  (Joseph,  Wiscomb, 
and  Weinman,  1976;  Wiscomb,  1977),  it  can  be  shown  that  the 
fourth  Legendre  coefficient  divided  by  9.0  is  a 
representation  of  the  total  amount  of  energy  forward 
scattered.  The  value,  9.0,  results  from  the  use  of  the 
fourth  term  (n*4)  of  the  approximate  expression  which 
contains  the  factor  (2n+l).  Based  on  the  study  by  Hansen 
(1969)  and  for  simplicity,  only  two  components  of  the 
scattered  radiation,  the  forward  scattered  and  the  backward 
scattered  components  will  be  considered.  Scanning  the  table 
of  expansion  coefficients  of  Welch  and  Cox,  one  can  see  the 
expansion  coefficients  are  slowly  varying  just  like  the 
droplet  extinction  and  absorption  coefficients.  Using  the 
same  reasoning  as  with  the  droplet  extinction  and  absorption 
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Table  2.7.2.B  (After  Welch,  Cox,  and  Davis,  1980) 


Droplet  Extinction  Coefficients  (Per  Kilometer) 


Wavelengths  (pm) 


0.95 

1.15 

1.4 

1 .85 

2.8 

3.35 

6.3 

SC 

base 

29.02 

29.44 

29.80 

30.20 

31.34 

31.69 

34.92 

sc 

top 

90.13 

91.13 

91.59 

92.90 

94.58 

95.83 

105.10 

NS 

base 

40.82 

41.22 

41.58 

42.30 

43.41 

44.05 

50.38 

NS 

top 

101.80 

102.80 

103.20 

104.60 

106.40 

107.70 

115.90 

ST 

base 

24.35 

24.68 

25.02 

25.63 

26.61 

26.81 

28.46 

ST 

top 

53.51 

54.06 

54.49 

55.23 

56.54 

57.32 

64.39 

CU 

45.17 

45.41 

45.95 

46.55 

47.97 

48.61 

54.51 

Droplet  Absorption  Coefficients  (Per  Kilometer) 

Wavelengths  (jim) 


0.95 

1.15 

1.4Q 

1.85 

2.8 

3.35 

6.3 

SC 

base 

0.004 

0.011 

0.204 

0.566 

15.18 

15.25 

11.91 

SC 

top 

0.022 

0.062 

1 .090 

2.960 

46.88 

46.27 

45.68 

NS 

base 

0.008 

0.019 

0.369 

0.940 

21.44 

21.35 

18.06 

NS 

top 

0.028 

0.077 

1 .390 

3.73 

52.46 

51.74 

52.34 

ST 

base 

0.0033 

0.0095 

0.162 

0.45 

12.66 

12.80 

9.64 

ST 

top 

0.0121 

0.037 

0.581 

1.45 

27.94 

27.72 

25.21 

CU 

0.008 

0.028 

0.406 

1.12 

23.61 

23.49 

20.69 

SC  =  Stratocumulus 
ST  =  Stratus 
NS  =  Nimbostratus 
CU  =  Cumulus 
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Table  2.7.2.C 


Average  Droplet  Extinction  Coefficients  (Per  Kilometer) 


Wavelengths  (^m) 


0.95 

1.15 

K40 

1.85 

2.80 

3.35 

6.30 

59.575 

60.285 

60.695 

61.550 

62.96 

63.839 

70.01 

71.31 

72.01 

72.44 

73.45 

74.905 

75.875 

83.14 

38.93 

39.37 

39.755 

40.43 

41.575 

42.065 

46.42 

45.17 

45.41 

45.95 

46.55 

47.97 

48.61 

54.51 

Average  Droplet  Absorption  Coefficients  (Per  Kilometer) 

Wavelengths  (/ urn ) 


0.95 

1.15 

1.40 

1 .85 

2.8 

3.35 

6.30 

SC 

0.013 

0.0365 

0.647 

1.763 

31.03 

30.26 

28.795 

NS 

0.018 

0.048 

0.8795 

2.335 

36.95 

36.545 

35.200 

ST 

0.00765 

0.02325 

0.3715 

0.95 

20.3 

20.26 

17.425 

CU 

0.008 

0.028 

0.406 

1.120 

23.61 

23.49 

20.69 

SC  =  Stratocumulus 
ST  =  Stratus 
NS  =  Nimbostratus 
CU  =  Cumulus 
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Table  2.7.2.D  (After  Welch,  Cox,  and  Davis,  1980) 


Droplet  Phase  Function  Expansion  Coefficients 


Wavelengths  (^m) 


0.765 

0.95  1.15 

1 .40 

1  .85 

2.80 

3.35 

6.3 

SC  base 
SC  top 

5.305 

5.413 

5.250  5.179 

5.395  5.362 

5.133 

5.362 

5.051 

5.372 

7.199 

8.074 

6.926 

7.746 

5.762 

6.670 

NS  base 
NS  top 

5.337 

5.419 

5.294  5.248 

5.399  5.379 

5.205 

5.384 

5.169 

5.420 

7.555 

8.123 

7.258 

7.806 

6.066 

6.892 

ST  base 
ST  top 

5.268 

5.367 

5.206  5.149 

5.339  5.300 

5.096 

5.275 

4.999 

5.253 

7.016 

7.789 

6.713 

7.466 

5.596 

6.310 

CU 

5.360 

5.305  5.269 

5.243 

5.209 

7.661 

7.346 

6.186 

Table  2 

.7.2.E 

Average  Droplet 

Phase  Function 

Expansion 

Coefficients 

Wavelengths  (jim) 

0.765 

0.95  1.15 

1.40 

1.85 

2.80 

3.35 

6,30 

SC 

5.359 

5.3225  5.2705 

5.2475 

5.2115 

7.6365 

7.336 

6.216 

NS 

5.378 

5.3465  5.3135 

5.2945 

5.2945 

7.839 

7.532 

6.479 

ST 

5.3175 

5.2725  5.2245 

5.1855 

5.126 

7.4025 

7.0895 

5.953 

CU 

5.360 

5.305  5.269 

5.243 

5.209 

7.661 

7.346 

6.186 

SC  =  Stratocumulus 

ST  =  Stratus 

NS  =  N imbostratus 

CU  =  Cumulus 

»s 


CLOUD  DROPLET 
EXTINCTION  COEFFICIENTS 


CLOUD  DROPLET 
ABSORPTION  COEFFICIENTS 


CLOUD  DROPLET 
FORWARD  SCATTtRNG 
PERCENTAGE 


Figure  2.7.2.F 

The  spectral  cloud  droplet  extinction,  absorption,  and 
forward  scattering  percentage  as  interpolated  and  extrapolated 
from  Welch,  Cox,  and  Davis  (1980) .  The  wavelength  is  in 
micrometers . 


ST  -  Stratus 

SC  -  Stratocumulus 


NS  -  Nimbostratus 
CU  -  Cumulus 


65 


coefficients,  linear  interpolation  and  extrapolation  was 
used  to  determine  the  expansion  coefficients  at  wavelengths 
not  computed  by  Welch  and  Cox.  An  average  of  the  expansion 
coefficients  at  the  top  and  base  of  the  cloud  was  used  to 
represent  the  expansion  coefficient  for  the  entire  cloud. 
Since  only  the  forward  and  backward  scatter  is  being 
considered,  only  the  fourth  Legendre  phase  function 
expansion  coefficient  was  interpolated.  Table  2.7.2.D  shows 
the  original  values  of  the  fourth  phase  function  coefficient 
used  by  Welch  and  Cox.  Table  2.7.2.E  shows  the  averaged 
values  of  the  fourth  phase  function  coefficient  which  is 
used  to  represent  the  overall  cloud.  The  amount  of  backward 
scattered  radiation  was  determined  by  subtracting  the  amount 
of  forward  scattered  radiation  from  1.0.  The  spectral  cloud 
droplet  absorption  coefficients,  cloud  droplet  extinction 
coefficients,  and  cloud  droplet  forward  scattering 
percentage  as  linearly  interpolated  and  extrapolated  from 
the  data  of  Welch,  Cox,  and  Davis  (1980)  are  plotted  in 
Figure  2.7.2.F. 

2.7.3  Cloud  Geometry 

The  radiative  properties  of  clouds  are  deeply 
intertwined  with  the  geometry  of  the  cloud  field.  A  number 
of  authors  have  studied  the  effects  of  cloud  geometry  and 
shape  on  the  radiative  characteristics  of  the  clouds.  One  of 
the  more  frustrating  aspects  of  dealing  with  this  problem  is 
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the  lack  of  a  convenient  way  of  handling  complicated  cloud 
geometries.  Studies  with  rectangular  clouds  have  been 
performed  (Aida,  1977;  Welch  and  Zdunkowski,  1980)  and  as  a 
result  of  the  relatively  poor  performance  of  rectangular 
type  cloud  geometries,  some  researchers  have  compared  the 
effects  of  clouds  of  different  shapes  (Welch  and  Zdunkowski, 
1981;  Welch  and  Wielicki,  1984). 

A  number  of  interesting  conclusions  have  been  drawn 
from  the  studies.  Aida  (1977)  demonstrated  that  the 
importance  of  cloud  shape  increases  with  increasing  solar 
zenith  angle  and  that  clouds  which  are  separated  by  more 
than  5  times  their  diameter  do  not  alter  the  scattering 
pattern  from  a  single  cloud.  Aida  further  concluded  that  the 
surrounding  clouds  increase  the  fraction  of  reflection  as 
the  distance  between  the  clouds  decreases.  This  effect 
becomes  particularly  evident  at  about  two  cloud  diameters 
separation,  roughly  equivalent  to  a  cloud  amount  of  25 
percent.  Welch  and  Zdunkowski  (1980)  studied  the  cloud  size 
distribution  in  comparison  to  cloud  reflectivity.  They  found 
that  the  cloud  size  which  contributes  most  to  the  cloud 
field  reflectivity  is  near  the  cloud  size  which  contributes 
most  to  the  total  sky  cover.  This  indicates  that  describing 
a  field  of  clouds  as  one  where  all  the  clouds  have  the  same 
dimensions  may  produce  accurate  results  if  the  cloud 
dimensions  are  carefully  chosen.  This  study  did  not  include 
cloud-cloud  radiative  interactions.  However,  the  result 
suggested  cloud-cloud  interactions  may  be  neglected  for 
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cloud  fields  with  30  percent  or  less  cloud  cover.  Welch  and 
Zdunkowski  (1980)  found  that  cloud  droplet  number  density 
has  a  large  influence  on  cloud  reflectivity  and 
transmissivity  and  a  minimal  effect  on  cloud  absorptivity. 
For  very  large  cumulus  clouds,  a  significant  fraction  of  the 
incident  radiation  exits  the  cloud  sides  with  a  majority  in 
the  downward  direction.  At  all  but  the  largest  solar  zenith 
angles,  the  radiation  exiting  the  cloud  sides  is 
anisotropic.  Welch  and  Zdunkowski  (1981)  studied  the  effects 
of  several  cloud  shapes  and  concluded  that  rectangular 
geometries  produce  results  which  are  marginally  accurate  for 
cylindr ically  shaped  clouds.  They  also  showed  that  the  exact 
geometry  should  be  taken  into  account  to  obtain  accurate 
radiative  characteristics. 

Wendling  (1977)  studied  the  effects  of  striated  cloud 
fields  and  his  conclusions  indicated  the  effects  were 
primarily  geometry  dependent.  The  radiance  emerging  from 
striated  clouds  can  be  greater  than  that  reflected  from  a 
plane-parallel  cloud  with  a  thickness  equal  to  the  thickness 
of  the  striated  cloud.  This  is  due  to  the  backscatter ing 
from  the  vertical  walls  of  the  cloud  columns  occurring  at 
intermediate  sun  elevation  angles.  By  the  same  token,  the 
radiance  reflected  from  the  striated  clouds  is  reduced 
compared  to  that  of  a  plane-parallel  cloud  when  shadowing 
occurs.  Wendling  pointed  out  that  plane-parallel  cloud 
albedos  are  a  good  approximation  for  striated  cloud  albedos 
only  when  the  striations  are  not  tco  deep  and  the  sun 
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illuminates  the  clouds  at  low  angles. 

Welch  and  Wielicki  (1984)  compared  the  effect  of  cloud 
shape  on  the  reflected  flux  differences  between  broken  and 
plane-parallel  cloudiness.  The  results  indicated 
plane-parallel  cloud  calculations  are  not  accurate  at  most 
values  of  cloud  cover,  particularly  those  with  cloud  cover 
between  10  and  90  percent.  They  also  found  the  choice  of 
cloud  shape  led  to  large  differences  in  the  reflected 
fluxes.  The  differences  were  attributed  to  the  anisotropic 
intensity  pattern  out  of  the  cloud  sides,  the  shape  of  the 
holes  (spaces)  between  the  clouds,  and  to  variations  in 
cloud  area  as  viewed  from  the  solar  direction.  Based  on  this 
finding,  Welch  and  Wielicki  were  able  to  develop  an 
expression  for  the  effective  cloud  cover  for  a  zenith  angle 
of  60  degrees.  Using  this,  they  computed  the  broken  field 
reflected  fluxes  from  plane-parallel  calculations.  They  were 
able  to  make  these  calculations  regardless  of  the  assumed 
cloud  shapes.  This  information  gives  credence  to  the 
techniques  which  weight  the  areas  of  the  clouds  by  the 
relative  proportions  which  have  been  illuminated  by  the  sun. 

Based  on  the  study  by  Welch  and  Wielicki  (1984),  the 
decision  was  made  to  use  a  simple  cloud  geometry  and  adjust 
the  effective  cloud  cover  based  on  solar  zenith  angle  in  an 
attempt  to  obtain  cloud  fluxes  which  represent  broken  cloud 
fields  using  plane-parallel  calculations.  Based  on  the  study 
by  Welch  and  Zdunkowski  (1980),  cloud  fields  containing  a 
single  cloud  size  will  be  used. 
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2.7.4.  The  Cloud  Model 

Based  on  the  paper  by  Welch  and  Wielicki  (1984),  the 
simplest  cloud  geometry  possible  is  used.  The  clouds  in  this 
model  are  rectangular  in  shape.  I  believe  that  adjusting  or 
weighting  the  areas  of  the  top  of  the  cloud,  the  side  of  the 
cloud  and  the  clear  sky  with  the  amounts  of  radiation  from 
plane-parallel  calculations  it  is  possible  to  obtain 
accurate  answers  regardless  of  the  cloud  shape  and  geometry. 
A  description  of  the  cloud  model  I  developed  follows. 

The  clouds  are  basically  rectangular.  For  rectangular 
clouds,  the  base  of  the  clouds  are  assumed  to  have  the  same 
dimensions  in  both  length  and  width  (a  square  base).  The 
height  of  the  cloud  base  and  the  separation  between  the 
sidewalls  of  the  clouds  are  determined  by  the  user.  The 
cloud  separation  is  assumed  to  be  uniform  in  the  directions 
perpendicular  to  the  sidewalls  of  the  clouds.  When  the 
length  of  the  cloud  base,  the  height  of  the  cloud,  the 
separation  between  the  clouds,  and  the  height  of  the  cloud 
base  are  specified,  the  geometry  associated  with  an  array  of 
rectangular  clouds  is  determined.  In  addition  to  an  array  of 
rectangular  clouds  (which  has  the  flavor  of  a  typical 
cumulus  field),  the  model  allows  for  a  striated  cloud  field. 
In  this  geometry,  the  clouds  are  formed  in  infinite  bands 
whose  width  and  separation  between  the  bands  is  defined  by 
the  user.  This  model  is  capable  of  handling  four  distinctly 


70 


different  geometries: 


1. 

Clear  sky  (no  clouds) 

2. 

Rectangular  array 

of  clouds 

3. 

Striated  array  of 

Clouds 

4. 

Overcast  sky  ( 100 
cover ) 

percent  cloud 

For  geometric  simplicity,  the  sun's  rays  are  assumed  to  be 
parallel.  Only  one  face  (side)  of  the  rectangular  cloud  is 
illuminated.  It  is  assumed  that  this  cloud  orientation  is 
maintained  throughout  the  day  when  the  model  computes  daily 
values.  The  top  of  the  cloud  is  always  illuminated  between 
sunrise  and  sunset.  The  portion  of  the  side  of  the  cloud 
which  is  directly  illuminated  by  the  sun  is  determined  by 
the  solar  zenith  angle.  Shadowing  of  the  sides  of  the  clouds 
is  allowed. 

2.7.4. 1  Cloud  Model  Geometry 

One  of  the  critical  factors  which  must  be  determined  is 
the  path  distance  that  radiation  must  travel  through  the 
cloud.  The  cloud  has  two  distinct  components  to  its 
radiation  field.  One  of  the  components  consists  of  radiation 
which  first  passes  through  the  top  of  the  cloud  and 
eventually  emerges  at  the  level  of  the  cloud  base  through 
either  the  side  or  bottom  of  the  same  (or  another)  cloud. 
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The  second  component  consists  of  radiation  which  first 
passes  through  the  side  of  the  cloud  and  eventually  emerges 
at  the  level  of  the  cloud  base.  The  geometric  light  ray, 
which  would  travel  the  paths  described,  may  cut  through  any 
number  of  individual  clouds  and  spaces  between  the  clouds 
depending  on  the  solar  zenith  angle  and  the  dimensions  of 
the  clouds.  The  average  equivalent  path  distances  through 
the  cloud  only  are  computed  numerically. 

The  definition  of  average  equivalent  path  distance  is 
the  distance  through  cloud  that  a  light  ray  would  travel 
beginning  at  the  level  of  the  top  of  the  cloud  and  exiting 
at  the  level  cloud  base.  The  path  distances  will  vary 
depending  on  whether  the  top  of  the  cloud  is  entered  first 
or  whether  the  side  of  the  cloud  is  entered  first.  It  will 
also  be  dependent  on  the  solar  zenith  angle.  The  solar 
zenith  angle  and  the  dimensions  of  the  clouds  will  determine 
how  much  shadowing  of  the  cloud  sides  will  occur,  how  many 
complete  clouds  are  intersected  by  the  path,  and  whether  the 
path  through  the  top  of  the  cloud  is  longer  or  shorter  than 
the  path  through  the  side  of  the  cloud.  Obviously,  there  are 
many  possible  paths  through  the  clouds  in  the  manner 
described.  The  average  path  distances  were  computed  by 
calculating  100  distinct  paths  (50  through  the  cloud  top, 
and  50  through  the  cloud  side)  and  averaging  the  path 
distances  for  the  respective  parts  of  the  cloud.  Figures 
2.7.4. 1. A  and  B  show  the  two  possible  average  equivalent 
path  distances  computed  by  the  model.  The  transmission  and 
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Figure  2. 7. 4.1 

The  equivalent  path  distances  for  rays  entering  the  side  and  the 
top  of  the  cloud  are  indicated  by  the  bold  lines  inside  the 
cloud  borders.  An  average  equivalent  path  distance  is  used  in  the 
model . 


reflection  from  the  clear  sky  area  between  the  clouds  will 
be  determined  by  the  clear  sky  portion  of  the  model. 


2. 7. 4. 2  Cloud  Model  Radiative  Transfer 

Multiple  scattering  contributes  to  both  the  upward  and 
downward  radiation  fluxes.  Neglecting  multiple  scattering 
leads  to  an  underestimate  of  the  radiation  leaving  the  cloud 
base  and  ultimately  an  underestimate  of  the  radiation 
reaching  the  surface  of  the  Earth.  The  same  is  true  of  the 
cloud  reflectivity. 

The  necessity  of  including  multiple  scattering  led  to 
the  review  of  solutions  to  the  basic  radiative  transfer 
equation.  A  simple,  analytic  solution  to  the  equation  of 
radiative  would  conform  with  the  original  intent  of 
developing  this  model.  After  reviewing  the  solutions  in  the 
literature,  two  techniques  produce  simple,  analytic 
solutions,  the  two  stream  technique  and  the  Eddington 
technique.  Both  solutions  handle  optically  thin  cases  and 
large  absorption  poorly.  I  chose  the  solution  developed  by 
Chu  and  Churchill  published  by  the  Radiation  Commission 
(1975)  in  Standard  Procedures  To  Compute  Atmospheric 
Radiative  Transfer  In  a  Scattering  Atmosphere  ;  the  solution 
was  recommended  by  Paltridge  and  Platt  (1976).  This  solution 
uses  the  two  stream  technique  of  decomposing  the  radiation 
field  into  components  (one  upward;  one  downward)  to  arrive 
at  an  approximate  solution  for  radiation  passing  through  a 
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plane-parallel  slab  of  a  given  thickness.  The  solutions 
allow  for  the  computation  of  the  reflectivity  or  albedo  at 
the  top  of  the  plane-parallel  layer  and  the  transmission 
through  the  plane-parallel  layer.  The  solutions  for  the 
albedo  and  transmission  of  the  layer  are: 


Albedo  = 


-2Kt 

G  (1.0  -  e _ )_ 

2  -2Kt 
1.0  -  G  e 


Transmission  = 


(1.0  -  G  ) 


1  .0  -  G 


where 


G  =  ( r  -  s )  /  ( r  +  s) 


r  =  1.0  -  W  f  +  W  b 
o  o 


K  =  s  / 


Mo 


-Kt 

e 


-2Kt 


2  2 


1/2 


s  =  [  (1.0  -  W  f  )  -  W  b  1 

o  o 


f  -  percentage  of  radiation  forward  scattered 
b  -  percentage  of  radiation  backscattered 
wq  -  single  scattering  albedo 
-  t  cos(solar  zenith  angle) 
t  -  thickness  of  the  layer 


The  use  of  this  solution  was  suggested  by  Paltridge  and 
Platt  (1976).  The  two  stream  approximation  is  generally  used 
in  applications  where  the  scatter  is  not  highly  anisotropic. 
However,  the  solutions  are  very  valuable  because  of  the 
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simplicity  of  the  analytic  solution.  In  addition,  this 
approach  can  be  justified  within  the  accuracy  of  the  model. 

According  to  Paltridge  and  Platt  (1976),  "a  method  of 

parameterizing  cloud-radiation  interaction  on  the  basis  that 

the  prediction  or  observation  of  cloud  amount  and  cloud 

character  is  inherently  inaccurate  to  much  worse  than  10 

percent.  Generally  they  (the  solutions)  give  answers  to 

better  than  about  5  percent  for  quite  anisotropic  scattering  * 

when  absorption  is  negligible  and  when  the  layer  optical 

depth  is  greater  than  about  5.  This  last  is  usually  the  case 

for  water  clouds  in  the  low  atmosphere,  but  it  is  certainly 

not  true  for  aerosol  layers  and  thin  cirrus."  The  solutions 

assume  a  non-reflecting  surface  or  layer  immediately  below 

the  atmospheric  layer  in  question.  To  account  for  a 

non-black  surface  below  the  layer  requires  an  additional 

technique.  The  complete  radiative  transfer  problem  and  its 

solution  for  the  plane-parallel  case  is  reproduced  from  the 

Radiation  Commission  (1975)  in  Appendix  D. 

In  addition  to  multiple  scattering  within  the  cloud 
itself,  multiple  scattering  (reflections)  between  the  cloud, 
atmosphere,  and  ground  are  important.  The  multiple 
reflections  play  a  key  role  in  determining  the  intensity  of 
the  radiation  which  reaches  the  surface  of  the  Earth. 

Neglecting  multiple  reflections  between  the  clouds  and  the 
Earth's  surface  will  lead  to  an  underestimate  of  the 
downward  flux  of  solar  radiation  reaching  the  Earth's 
surface  (Schneider  and  Dickinson,  1976).  According  to 
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Schneider  and  Dickinson,  the  underestimate  would  be  most 
pronounced  in  regions  of  persistent  cloud  cover  and  high 
surface  albedos.  To  include  the  effect  of  the  multiple 
reflections  between  the  surface,  the  clouds  and  the 
atmospheric  layers,  the  adding  technique  was  used  (Appendix 
B).  The  adding  technique  (Liou,  1980)  combines  two  layers  of 
known  reflectivity  and  transmissivity  through  a  series 
approximation  of  the  multiple  reflections  which  occur 
between  the  layers. 

2. 7. 4. 3  The  Atmospheric  Layers 

The  atmosphere  will  be  divided  into  three  layers  plus  the 
surface  (the  fourth  layer).  The  layer  above  the  clouds  will 
be  treated  as  clear  sky.  The  second  layer  will  contain  the 
clouds.  The  layer  below  the  clouds  and  treated  as  clear  sky. 
The  bottom  layer  will  be  the  ground  or  surface.  The 
transmissivity  and  reflectivity  of  the  second  or  cloudy 
layer  will  be  computed  as  described  in  the  cloud  radiative 
transfer  section.  The  transmissivity  and  reflectivity  of  the 
clear  sky  layers  above  and  below  the  clouds  must  be  computed 
separately. 

The  amount  of  atmosphere  in  the  layers  above  and  below 
the  clouds  affects  the  transmissivity  and  reflectivity  of 
the  individual  layers.  The  temperature,  the  pressure  at  the 
top  of  the  cloud,  and  the  pressure  at  the  base  of  the  cloud 
affect  the  amount  of  atmosphere  in  each  layer.  The  user 
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determines  the  surface  temperature,  surface  pressure,  height 
of  the  cloud  base,  and  the  thickness  of  the  cloud  for  the 
model.  From  these  values,  the  atmospheric  pressure  at  the 
base  and  top  of  the  cloud  layer  is  computed.  Based  on  these 
reference  pressures,  the  transmissivity  and  reflectivity  of 
the  atmospheric  layers  are  determined.  The  calculation  of 
the  reference  pressures  is  performed  via  the  formula  for 
constant  lapse  rate  atmospheres  (Hess,  1959;  Iribarne  and 
Godson,  1973). 

P  =  P  (T  /  T  )(§/R~)  =  p  [  (T  -«*)/T 
oo  o  o  o 


where 

g  -  Gravitational  constant 

P  -  Pressure  at  height  z  above  the  surface 

P  -  Pressure  at  the  surface 
o 

T  -  Temperature  at  height  z  above  the  surface 
Tq  -  Temperature  at  the  surface 
R  -  Universal  Gas  Constant 
z  -  Height  above  the  surface 
3  -  Constant  atmospheric  lapse  rate 

Some  assumptions  are  implicit  in  the  model,  when  using  this 
formula.  First,  the  U.S.  Standard  Atmosphere  tropospheric 
lapse  rate  (-6.5  degrees  Kelvin  /  Km)  is  used.  This  agrees 
with  the  assumed  lapse  rate  in  the  cloud  parameter izations 
of  Welch,  Cox,  and  Davis  (1980).  The  atmosphere  is  assumed 
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to  have  a  constant  lapse  rate  to  the  tropopause.  Second,  it 
is  assumed  that  the  clouds  exist  below  the  tropopause.  This 
model  is  not  designed  for  clouds  which  break  through  the 
tropopause  as  can  be  the  case  in  some  thunderstorm 
situations  or  with  noctilucent  clouds.  Both  of  these 
assumptions  are  very  good  for  most  of  the  usual  atmospheric 
conditions  encountered. 

Since  the  transmission  for  an  atmospheric  layer  is 
computed  using  an  exponential  function,  the  exponent  of  each 
transmission  function  was  weighted  by  a  ratio  of  the 
pressures  to  give  the  correct  exponent.  For  example,  the 
total  transmission  through  three  atmospheric  layers  might  be 
written  as  an  exponential  function  whose  exponent  includes 
the  product  of  an  extinction  coefficient  (k)  and  an 
atmospheric  path  distance  (M) .  The  total  transmission  (T) 
for  an  atmospheric  layer  could  be  written  as  : 


This  single  layer  could  be  broken  down  into  component  layers 
if  we  assume  that  the  product  of  the  transmittances  gives 
the  total  transmittance.  While  this  technique  could  be  very 
poor  for  molecular  absorption  and  scattering,  it  is  often 
used.  Iqbal  summarized  the  reasoning  very  succinctly.  "In 
order  to  compute  the  direct  spectral  irradiance  on  the 
Earth,  we  need  the  values  of  the  monochromatic 
transmittances  due  to  the  various  molecular  absorbers. 


79 


Molecular  scattering  depends  primarily  on  the  number  density 
of  molecules  in  the  path,  whereas  the  molecular  absorption 
is  a  function  of  local  pressure  and  temperature  as  well.  The 
wavelength  dependence  of  molecular  (Rayleigh)  scattering  is 
very  nearly  proportional  to  On  the  other  hand,  the 
spectral  variation  of  the  molecular  absorption  coefficient 
is  a  highly  oscillatory  function  of  wavelength.  For 
molecular  absorption,  it  is  necessary  to  know  the 
frequencies,  intensities,  shapes,  and  widths  of  spectral 
lines  in  the  region  of  interest  in  order  to  be  able  to 
evaluate  the  spectral  transmittance.  The  simple  Bouguer's 
law  does  not  apply  in  this  case.  Nevertheless,  since 
elaborate  procedures  require  a  great  deal  of  computer  time, 
it  is  quite  common  to  apply  this  law  or  a  variation  of  it  to 
compute  molecular  absorption."  (Iqbal,  1983).  Bird  made  this 
assumption  in  choosing  the  form  of  the  transport  equation 
for  his  clear  sky  model.  In  a  similar  fashion,  each  layer 
will  be  formulated  using  this  assumption.  If  each  layer  is 
represented  by  Tx  then  the  total  transmittance  for  the 
atmosphere  could  be  represented  by: 
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T  T 
1  2 
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Each  layer  can  be  represented  by: 
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T_  =  e 


•[(P  k-P„J/PJKM 

CD  Ct  O 


-t<Pn~P^K)/P^]KM 
0  CD  O 


where 


T  -  Transmittance  of  each  layer 
x 

P  ,  -  Pressure  at  the  cloud  base 
cb 

P  -  Pressure  at  the  cloud  top 
ct 

P  -  Pressure  at  the  surface 


(x  =  1,  2,  3) 


Since  the  exponents  in  a  product  are  additive,  the  product 
of  the  resultant  layer  transmissivities  yields  the  original 
total  transmissivity  for  the  entire  layer. 
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T  T  T 
1  2  3 


-KM 
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This  is  the  procedure  used  to  compute  the  atmospheric 
transmissivities  for  the  individual  layers  of  the 
atmosphere.  This  technique  is  applied  to  each  atmospheric 
absorber,  scattering  process,  or  constituent  in  turn.  For 
example,  the  transmission  for  aerosols  was  computed  by 
layer,  the  transmission  by  the  uniformly  mixed  gases  was 
computed  by  layer,  etc.  The  transmission  for  water  vapor  was 
accomplished  in  an  identical  fashion  except  that  only  two 
atmospheric  layers  were  used.  It  was  assumed  that  the  water 
vapor  in  the  atmosphere  exists  below  the  tropopause  or  the 
top  of  the  cloud  (whichever  is  lower).  This  is  a  good 
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assumption  under  most  of  the  normal  meteorological 
conditions.  Making  this  assumption  means  that  the  highest  of 
the  three  component  atmospheric  layers  used  in  this  model 
has  its  base  at  the  top  of  the  tropopause  or  at  the  top  of 
the  cloud.  Another  assumption  is  that  all  the  absorption  by 
ozone  takes  place  above  the  tropopause  or  the  top  of  the 
cloud.  This  means  the  terms  which  account  for  the  effects  of 
ozone  will  be  in  the  top  layer  and  the  terms  which  account 
for  the  effects  of  water  vapor  will  only  be  in  the  lower  two 
layers.  Once  the  transmission  for  each  of  the  components  has 
been  broken  down  by  layer,  the  resultant  layer  transmission 
was  computed  by  multiplying  the  individual  layer  components 
together  in  the  formulation  devised  by  Iqbal  (See  Figure 
2. 4. A).  Iqbal's  formulation  was  modified  slightly  by  layer 
depending  on  the  assumptions  mentioned  above.  The  resultant 
transmission  calculations  for  the  clear  sky  layers  are: 


Layer  Above  The  Clouds: 

T ,  =  T  T  T  T  ,,  +  T  T  T  v, 

1  rAl  aAl  oAl  uAl  aAl  oAl  uAl 


U.°  -  TrXl)  0 


+  T 


oAl  TrAl  Tu\l  (1-°  "  ^Al*  Wo 


Layer  With  The  Clouds  (Clear  Air  Portion  Only) : 

T2  =  TrA2  TaA2  TwA2  TuA2  +  TaA2  TuA2  TwA2 


+  TrA2  TuA2  TwA2  {1*°  -  W  "o 


'  TrA2 
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Layer  Below  The  Clouds: 


T3  TrX3  TaX3  TwA3  Tu\3  +  TaA3  TuA3  TwA3  (1'°  TrA3)  °'5 


+  TrX3  TaX3  TwA3  (1‘°  ~  W  Wo  Fa 


All  of  the  quantities  in  the  above  equations  have  been 
previously  defined  with  the  exception  of  the  additional 
subscript  to  indicate  the  level  of  the  atmosphere  and  will 
not  be  redefined  here. 

The  reflectivities  for  each  layer  were  computed 
similarly  and  the  equations  from  Iqbal  were  modified  by 
layer  to  include  or  exclude  water  vapor  and  ozone  as 
appropriate.  One  might  wonder  why  go  to  such  trouble  for  a 
simple  model.  The  answer  lies  in  the  combining  of  the 
atmospheric  layers.  The  explanation  is  easiest  if  one 
studies  the  equations  of  the  adding  technique  for  combining 
atmospheric  layers.  One  can  see  that  the  order  of  combining 
the  atmospheric  layers  makes  absolutely  no  difference  to  the 
total  transmissivity  calculated.  However,  the  order  will 
directly  affect  the  atmospheric  reflectivity  for  the 
combination  of  the  layers.  For  this  reason,  it  was  very 
important  to  insure  that  the  regions  of  dominant  influence 
of  water  vapor  and  ozone  were  constructed  in  the  proper 
layer.  The  adding  technique  of  combining  layers  requires 
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both  the  total  (direct  and  diffuse  combined)  transmittance 
and  the  total  reflectance  of  each  layer.  At  the  surface, 
only  the  reflectivity  is  known. 

To  account  for  the  multiple  scattering  effects  of  water 
vapor,  the  uniformly  mixed  gases,  etc.  in  the  layer 
containing  the  cloud,  the  transmission  of  the  layer  was 
calculated  as  if  no  clouds  existed.  The  layer  transmission 
was  equated  to  an  exponential  and  solved  to  obtain  an 
extinction  coefficient  for  the  'clear  sky'  in  this  layer. 
This  extinction  coefficient  was  added  to  the  cloud  droplet 
extinction  coefficient  to  obtain  the  extinction  coefficient 
for  the  entire  layer.  The  layer  extinction  coefficient  was 
used  in  the  two  stream  solution  which  accounts  for  multiple 
scattering.  In  this  way,  the  effects  of  multiple  scattering 
were  applied  to  all  the  constituents  in  the  layer  containing 
the  clouds. 

2. 7.4.4  Combining  The  Atmospheric  Layers  And  The  Surface 
Reflectivity 

Once  the  reflectivity  and  transmissivity  of  each  layer 
(including  the  surface)  have  been  determined,  then  the 
layers  are  combined  using  the  adding  technique,  two  layers 
at  a  time.  The  adding  technique  produces  four  equations  for 
combining  the  layers.  The  equations  reperesent: 

1.  The  radiation  reflected  from  the  combination  of 
the  two  layers  (R). 
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2.  The  radiation  transmitted  by  the  combination  of 
the  two  layers  (T). 

3.  The  radiation  reflected  upward  onto  the  bottom 
of  the  top  layer  (U). 

4.  The  radiation  transmitted  downward  onto  the  top 
of  the  bottom  layer  (D). 


Appendix  B  develops  the  equations  for  each  situation.  For 
the  incoming  solar  radiation,  the  layers  are  combined  from 
the  top  down,  two  layers  at  a  time.  By  choosing  the  desired 
equation  from  the  adding  technique  at  each  step,  the 
radiation  reflected  from  the  atmosphere-surface  combination 
as  well  as  the  radiation  incident  on  the  surface  can  be 
obtained. 

These  steps  are  accomplished  for  two  different  cases, 
from  the  different  equivalent  path  lengths  through  the  side 
and  the  top  of  the  cloud.  The  solution  to  the  radiative 
transfer  equation  assumes  a  plane-parallel  layer.  The 
radiation  passing  through  the  side  of  the  cloud  is  treated 
as  if  it  was  a  plane-parallel  layer  with  a  thickness  equal 
to  the  equivalent  path  distance  through  the  side  of  the 

cloud.  The  radiation  passing  through  the  top  of  the  cloud  is 

/ 

treated  similarly.  The  model  produces  three  sets  of 
reflectances  and  transmittances  indicative  of  the  entire 
atmosphere.  The  three  sets  of  reflectances  and 
transmittances  represent: 


1.  the  clear  sky  atmosphere 
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2.  the  combination  of  the  clear  layer  above  the 
cloud,  the  radiation  passing  through  that  top  of 
the  cloud,  the  clear  layer  below  the  cloud,  and 
the  surface  reflectivity 

3.  the  combination  of  the  clear  layer  above  the 
cloud,  the  radiation  passing  through  the  side  of 
the  cloud,  and  the  clear  layer  below  the  cloud, 
and  the  surface  reflectivity. 


Using  the  irradiance  falling  on  the  top  of  the  atmosphere, 
the  radiation  reaching  the  surface  of  the  Earth  can  be 
computed  for  the  three  cases.  The  radiation  reaching  the 
surface  of  the  Earth  from  the  three  cases  is  weighted  by  the 
relative  area  each  represents  on  a  horizontal  surface.  As 
suggested  by  Welch  and  Wielicki,  performing  this  calculation 
enables  one  to  compute  the  effects  of  broken  cloud  cover 
regardless  of  the  cloud  shape.  For  clear  sky  cases,  the  Bird 
model  values  will  be  obtained.  For  overcast  sky  cases,  the 
cloud  layer  will  be  treated  as  plane-parallel  with  the 
radiation  entering  the  top  of  the  clouds. 

2. 7. 4. 5  Satellite  Observation  Of  The  Surface 

For  a  satellite,  the  process  described  is  repeated  for 
a  different  set  of  equivalent  path  distances  through  the  top 
and  side  of  the  cloud  as  well  as  a  different  air  mass  for 
the  clear  sky.  These  differences  are  the  result  of  viewing 
geometry.  All  the  layer  transmissivities  and  reflectivities 
are  computed  again  using  the  new  atmospheric  path  distances. 
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Because  the  radiation  is  coming  from  the  surface,  the 
atmospheric  layers  are  combined  from  the  bottom  (surface)  to 
the  top  of  the  atmosphere.  The  exact  opposite  of  the 
incoming  radiation.  Once  the  transmissions  for  the 
atmosphere  have  been  recomputed,  the  radiation  leaving  the 
surface  of  the  Earth  is  taken  as  the  source  irradiance  and 
is  used  to  compute  the  irradiance  arriving  at  the  top  of  the 
atmosphere.  The  radiation  arriving  at  the  top  of  the 
atmosphere  is  weighted  according  to  the  relative  areas  of 
the  clear  sky,  top,  and  side  of  the  cloud  as  viewed  from  the 
satellite  to  obtain  the  radiation  reaching  .he  satellite. 

2. 7. 4. 6  Direct  And  Diffuse  Radiation  Beneath  The  Cloud 

The  adding  technique  uses  the  total  transmittance 
(direct  plus  diffuse)  of  the  layers  and  computes  a  total 
transmittance  for  the  combination  of  layers.  At  the  surface, 
the  values  of  direct  and  diffuse  radiation  are  desired.  To 
determine  the  direct  radiation  portion  of  the  total  combined 
transmittance  including  the  cloud,  a  simple  Beer's  law 
approximation  is  used.  The  Beer’s  law  solution  is  for  the 
direct  radiation  (plus  a  small  amount  of  forward  scattered 
radiation).  The  Beer's  law  solution  has  the  form: 


Transmittance  =  e 


where 


t  -  optical  depth 
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By  subtracting  the  Beer's  law  solution  from  the  total 
solution  of  the  equation  of  radiative  transfer,  the 
proportion  of  direct  and  diffuse  components  can  be 
approximated.  To  summarize,  the  direct  and  diffuse  radiation 
components  beneath  the  cloud  are  computed  as  follows: 


Total 


(1.0 


G2) 


-Kt 

e 


1 .0 


-Kt 

e 


Direct  =  e 

Diffuse  =  Total  -  Direct 

where 

K  -  Absorption  coefficient 
t  -  Thickness  of  che  cloud 

G  =  (r-s)  /  (r+s)  (See  Section  2. 7. 4. 2) 

This  simple  technique  separates  the  direct  and  diffuse 
components  of  the  radiation  underneath  the  cloud.  This  same 
technique  determines  the  amount  of  radiation  from  the  target 
at  the  surface  which  reaches  the  satellite  observation  point 
directly. 

2.8  Radiation  Reaching  The  Surface 

The  radiation  from  the  clear  s Jcy  and  cloud  must  be 
weighted  by  their  respective  areas  to  get  accurate  results 
using  plane  parallel  solutions. 
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2.8.1  Weighting  Each  Component  Of  Radiation 

My  technique  for  weighting  the  three  components  of 
radiation  is  based  on  the  projected  area  of  each  component 
onto  the  ground  (See  Figure  2.8.1. A).  The  projected  area  of 
clear  sky  between  the  clouds  could  vanish  if  the  solar 
zenith  angle  and  the  dimensions  of  the  clouds  are  such  that 
shadowing  occurs.  However,  there  will  always  be  lanes  of 
clear  sky  in  the  array  of  rectangular  clouds.  There  may  not 
be  projected  areas  of  clear  sky  associated  with  a  striated 
cloud  field.  For  demonstration  purposes,  the  following 
discussion  will  consider  a  rectangular  array  of  clouds  with 
radiation  reaching  the  ground  in  the  clear  space  between  the 
clouds.  The  first  step  is  to  calculate  the  projected 
dimensions  of  the  areas  on  the  ground.  (Refer  to  Figure 
2.8. 1.B).  To  determine  the  weights  of  each  component  of 
radiation,  one  only  needs  to  compute  the  areas  on  the  ground 
relative  to  the  total  area.  For  this  particular  case,  the 
areas  turn  out  to  be: 


FADR 


FRACTIONAL  AREA  D  (W  +D)  +  W  ( D  -  H  tan  $  ) 


(  W  +  D)2 


TOTAL  AREA 
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H  TAN  8  D-  H  TAN  8 


Figure  2.8.1.B 

Fractional,  area  weights  for  the  components  of  radiation  through 
the  clear  sky  (FADR),  the  side  of  the  cloud  (FASC),  and  the  top  of 
the  cloud  (FATC) 


FADR 


D(W+D)  +  W(D-  H  tan  6 ) 
(W+D) 2 


W  ( H  tan  0  ) 

FASC  =  fw7 op 

W2 

FATC  =  ( W*  D )  ? 


A 
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W  (H  tan  9) 

FASC  =  _ 

(  W  +  D  )2 

w2 

FATC  =  _ 

(  W  *  D)2 

FADR  -  Fractional  area  of  direct  (clear  sky)  radiation 
FASC  -  Fractional  area  of  radiation  through  the  side  of 
the  cloud 

FATC  -  Fractional  area  of  radiation  through  the  cop  of  the 
c  loud 

W  -  Width  of  the  cloud 
D  -  Separation  between  the  clouds 
H  -  Height  of  the  clouds 

&  -  Solar  zenith  angle  relative  to  a  flat  surface 

Once  the  weighting  for  each  component  of  radiation  has  been 
determined,  the  radiation  reaching  the  surface  can  be 
computed. 

2.8.2  Calculating  The  Radiation  Reaching  The  Surface 

Each  of  the  components  of  the  radiation  {clear  sky, 
cloud  side,  and  cloud  top)  can  be  broken  into  direct-normal 
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and  diffuse  horizontal  radiation.  This  is  accomplished  for 
the  areas  underneath  the  clouds  by  the  Beer’s  law  technique 
described  in  Section  2. 7. 4. 6.  The  average  direct  normal  and 
diffuse  horizontal  radiation  can  be  computed  in  the 
following  way: 


where 


DN  =  (FADR  *  DNC)  +  ( FASC  *  DNSC)  +  ( FATC  *  DNTG ) 

DF  =  (FADR  *  DFC)  +  (FASC  *  DFSC)  +  (FATC  *  DFTC ) 

DF  -  Total  diffuse  radiation 

DFC  -  Diffuse  radiation  from  the  clear  sky 

DFSC  -  Diffuse  radiation  from  the  side  of  the  cloud 

DFTC  -  Diffuse  radiation  from  the  top  of  the  cloud 

DN  -  Total  direct-normal  radiation 

DNC  -  Direct-normal  radiation  from  the  clear  sky 

DNSC  -  Direct-normal  radiation  from  the  side  of  the  cloud 

DNTC  -  Direct-normal  radiation  from  the  top  fo  the  cloud 


Once  the  total  direct-normal  (DN)  and  the  total  diffuse  (DF) 
radiation  reaching  the  surface  have  been  computed,  then  the 
amount  of  radiation  reaching  an  inclined  plane  surface  can 
be  determined. 


2.8.3  Radiation  Arriving  At  An  Inclined  Plane 
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The  diffuse  radiation  is  assumed  to  be  isotropically 
distributed  throughout  the  sky.  This  is  not  true.  Making 
this  assumption  simplifies  the  solution.  Hooper  and  Brunger 
(1980)  compared  the  results  of  several  modelling 
assumptions.  Among  them  were  the  isotropic  assumption  and 
the  heliocentric  assumption  for  the  diffuse  radiation. 
According  to  Hooper  and  Brunger,  making  the  isotropic 
assumption  will  give  conservative  predictions 
(underestimates)  for  partially  cloudy  conditions,  but  it 
gives  good  approximations  for  overcast  conditions.  They 
indicate  that  the  isotropic  diffuse  assumption  will  result 
in  underestimates  of  the  diffuse  irradiance  on  sun-facing 
slopes,  and  overestimates  of  diffuse  irradiance  on  slopes 
not  facing  the  sun.  Under  varying  sky  conditions  and  a 
sloping  surface,  the  error  between  calculated  diffuse  using 
the  isotropic  assumption  and  the  measured  diffuse  can  be  as 
large  as  56  percent  (Hooper  and  Brunger,  1980). 

To  compute  the  radiation  arriving  at  an  inclined  plane 
surface,  the  percentage  of  the  sky  visible  to  the  plane  must 
be  computed.  To  compute  the  percentage  of  the  sky  visible  to 
the  plane,  the  ratio  of  the  surface  area  of  the  lune-shaped 
portion  of  the  sky  visible  to  the  inclined  plane  divided  by 
the  surface  area  of  a  hemisphere  was  calculated.  To  compute 
the  percentage  of  sky  and  ground  visible  to  the  inclined 
plane,  the  surface  areas  of  the  two  regions  shown  in  Figure 
2. 8. 3. A  need  to  be  determined. 
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Figure  2. 8. 3. A 

Regions  of  sky  and  ground  visible  to  inclined  plane 


95 


St  =  2  R 2  /}  S2  =  2*R2  -  2  R2  /J 

where 

-  Surface  area  of  the  sky  hemisphere  not  visible  to 
the  inclined  plane 

S2  -  Surface  area  of  the  sky  hemisphere  visible  to  the 
inclined  plane 

ji  -  Zenith  angle  of  the  inclined  plane  measured  relative 
to  the  normal  of  a  flat  surface 

At  this  point,  the  two  solid  angles  of  interest  can  be 
computed.  The  solid  angles  representing  the  portions  of  the 
inclined  plane's  hemisphere  visible  to  the  sky  and  visible 
to  the  ground  may  be  computed  as  follows: 


SASKY  = 


2  2 
2  R  -  2R  /} 


R 


2 


2(  t-H) 


SAGND 


S 


1 


R 


2 


2(  /J  ) 


The  direct  radiation  must  be  modified  by  the  solar 
zenith  angle  relative  to  the  inclined  plane  and  the  diffuse 
radiation  must  be  modified  by  the  percentage  of  the  sky 
visible  to  the  plane.  The  equations  for  the  direct  and 


diffuse  radiation  reaching  the  surface  are  : 

DR  =  DN  *  cos(  ZPLANE) 

DFS  =  OF  *  PCTSKY 

where 

DFS  -  Diffuse  radiation  reaching  the  inclined  plane  from 
the  sky 

DR  -  Direct  radiation  reaching  the  inclined  plane 
ZPLANE  -  Solar  zenith  angle  relative  to  the  inclined  plane 

Because  the  plane  is  inclined,  diffuse  radiation  is  received 
from  both  the  sky  and  the  ground.  The  ground  diffuse 
radiation  is  computed  in  the  following  manner: 

GNDDIF  =  [  DN  *  cos  (ZFLAT )  *  SFCRFL 
DF  *  GNDRFL] 

where 

GNDDIF  -  Diffuse  radiation  reflected  by  the  ground 
GNDRFL  -  Hemispheric  or  diffuse  reflectance,  not  angle 
dependent 

SFCRFL  -  Direct  radiation  surface  reflectance,  dependent  on 
the  solar  zenith  angle  relative  to  a  flat  surface 

The  'irst  term  is  the  reflected  portion  of  the  radiation 
directly  incident  on  the  horizontal  ground.  The  second  term 
is  the  reflected  portion  of  the  diffuse  radiation  incident 
on  the  horizontal  ground.  To  obtain  the  ground  diffuse 
contribution  to  the  radiation  incident  on  the  inclined 
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plane,  the  ground  diffuse  radiation  must  be  divided  by  the 
solid  angle  subtended  by  a  horizontal  surface  (Lambertian 
surface  assumed)  and  multiplied  by  the  solid  angle  of  the 
ground  as  viewed  from  the  inclined  plane.  The  global 
radiation  arriving  at  an  inclined  plane  is: 

GBL  =  DN  *  cos(ZPLANE) 

+  DF  *  PCTSKY 

+  [  DN  *  cos(ZFLAT)  *  SFCRFL 

+  DF  *  GNDRFL  ]  *  SAGND  / t 

where 

GBL  -  Global  radiation  arriving  at  the  inclined  plane 

which  includes  direct,  sky  diffuse,  and  ground 
•  diffuse  radiation 

SAGND  -  Solid  angle  subtended  by  the  ground  as  viewed 
from  the  inclined  plane 

ZPLANE  -  Solar  zenith  angle  relative  to  the  normal  vector 
of  the  inclined  plane 

PCTSKY  -  Percentage  of  sky  visible  to  inclined  plane 

To  compute  the  radiation  reaching  the  satellite  observation 
point,  the  radiation  leaving  the  surface  of  the  inclined 
plane  must  be  calculated.  It  is  evaluated  using  the 
following  equation: 
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SHZERO  =  (  DN  *  cos(ZPLANE)  *  PLNRFL 
+  DF  *  PCTSKY  *  GNDRFL 
+  [  DN  *  cos(ZFLAT)  *  SFCRFL  *  PLNRFL 
where  +  DF  *  GNDRFL  *  GNDRFL  ]  *  SAGND  /  tt  )  /t 

SHZERO  -  Radiation  leaving  the  inclined  plane  located  at 
the  surface  of  the  Earth 

PLNRFL  -  Direct  radiation  surface  reflectance  of  the 

inclined  plane,  dependent  on  the  solar  zenith 
angle  relative  to  the  normal  vector  of  the 
inclined  plane 

2.9  Satellite  Observation  Of  The  Surface 

The  solar  spectral  model  computes  the  irradiance 
reaching  a  satellite  in  orbit.  To  define  the  satellite 
viewing  geometry,  the  satellite  zenith  angle  relative  to  the 
horizontal  surface  at  the  inclined  plane’s  location  must  be 
given.  In  order  to  determine  the  solid  angles  needed  for  the 
computation  of  the  radiation  reaching  the  satellite,  the 
horizontal  extent  of  the  plane  surface  must  be  given.  The 
model  uses  a  fixed  viewing  geometry  of  the  satellite.  In 
other  words  the  viewing  geometry  does  not  change  relative  to 
the  inclined  plane  or  target  on  an  hourly  basis.  However, 
the  model  will  sum  the  total  radiation  reaching  the 
satellite  position  with  the  fixed  viewing  geometry 
throughout  the  day.  While  this  is  not  representative  of 
practical  situations  where  the  geometry  of  the  satellite 
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relative  to  a  fixed  ground  location  changes  with  each  orbit, 
it  will  give  an  estimate  of  the  radiation  which  might  be 
accumulated  at  satellite  altitude  through  the  course  of  a 
day.  The  user  may  get  around  this  by  computing  several 
instantaneous  or  hourly  cases  with  different  viewing  angles 
and  summing  them. 

Since  the  geometry  for  the  satellite  will  probably  be 
different  from  the  geometry  of  the  sun  relative  to  the 
inclined  plane,  all  the  calculations  concerning  the 
fractional  areas  of  clear  and  cloudy  skies,  the  equivalent 
path  distances  throughout  the  clouds,  and  the  zenith  angle 
relative  to  the  inclined  plane  must  be  calculated  for  the 
path  between  the  Earth  and  the  satellite.  Since  the  path 
distance  through  the  atmosphere  will  be  different,  the 
atmospheric  parameters  affected  by  the  atmospheric  path 
distance  must  be  computed  again.  Exactly  the  same  procedures 
for  computing  the  radiation  reaching  the  surface  from  the 
sun  have  been  followed  in  determining  the  radiation  reaching 
the  satellite.  The  only  difference  is  that  the  source  of 
radiation  is  the  radiation  leaving  the  inclined  plane’s 
surface  (SHZERO)  instead  of  the  extraterrestrial  solar 
spectrum  (HZERO) .  The  radiation  leaving  the  inclined  plane’s 
surface  is  assumed  to  be  a  direct  source  of  radiation.  The 
direct-normal  radiation  reaching  the  satellite  from  the 
target  plane  (DNO)  and  the  diffuse  radiation  reaching  the 
satellite  from  the  target  plane  (DFO)  constitute  the 
radiation  from  the  target  itself.  Scattered  radiation  from 
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the  atmosphere  also  reaches  the  satellite.  Depending  on  the 
position  of  the  sun  relative  to  the  satellite,  this 
radiation  may  be  mostly  forward  scattered  or  backscattered 
radiation  (See  Figures  A. 3  and  A. 4).  The  radiance  of  the 
atmosphere  is  calculated  based  on  whether  the  geometry 
indicates  forward  or  backward  scattering  of  radiation.  The 
radiance  of  the  atmosphere  includes  any  radiation  scattered 
by  the  atmosphere  before  reaching  the  ground  on  the  incoming 
path  from  the  sun  to  the  surface.  Radiation  scattered  by 
clouds  is  included  in  the  computation  of  the  radiance  of  the 
atmosphere.  Since  the  radiance  of  the  atmosphere  is  computed 
from  the  adding  technique  as  if  it  included  all  the 
radiation  scattered  over  an  entire  hemisphere,  it  was 
necessary  to  multiply  the  atmospheric  radiance  by  the  solid 
angle  of  the  plane.  It  is  assumed  the  radiation  reaching  the 
satellite  is  determined  by  the  instantaneous  f ield-of-view 
of  the  satellite. 

The  satellite  sensor  is  assumed  to  be  pointing  at  the 
target  so  there  is  no  equivalent  term  for  the  direct 
radiation  (DR)  for  the  satellite  as  with  the  surface 
radiation.  One  significant  difference  in  the  calculations  of 
the  incoming  and  outgoing  radiation  is  the  combination  of 
the  layers  of  the  atmosphere.  In  one  case,  the  layers  are 
combined  from  the  top  of  the  atmosphere  downward  and  in  the 
other  case  the  layers  are  combined  from  the  bottom  of  the 
atmosphere  upward.  A  second  significant  difference  is  the 
inclusion  of  a  term  to  account  for  the  radiance  of  the 
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atmosphere.  The  equations  for  the  radiation  reaching  the 
satellite  are: 


DNO  =  (FADRO  *  DNCO)  +■  (FASCO  *  DNSCO)  +  ( FATCO  *  DNTCO) 
+  LATM  *  SATGT 

DFO  =  (FADRO  *  DFCO)  +  (FASCO  *  DFSCO)  +  (FATCO  *  DFTCO) 


where 

DFCO 

DFSCO 

DFTCO 

DFO 

DNCO 

DNTCO 

DNO 

FADRO 


-  Diffuse  radiation  observed  by  the  satellite  in  the  clear 
sky 

-  Diffuse  radiation  observed  by  the  satellite  through  the 
side  of  the  cloud 

-  Diffuse  radiation  observed  by  the  satellite  through  the  top 
of  the  cloud 

-  Diffuse  radiation  observed  by  the  satellite  from  the  target 

-  Direct-normal  radiation  observed  by  the  satellite  in  the 
clear  sky 

-  Direct-normal  radiation  observed  by  the  satellite  through 
the  top  of  the  cloud 

-  Direct-normal  radiation  observed  by  the  satellite  from  the 
target 

-  Fractional  area  of  the  clear  sky  radiation  as  observed  from 


the  satellite 


I 
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FASCO  -  Fractional  area  of  the  radiation  through  the  side  of  the 
cloud  as  observed  from  the  satellite 


FATCO  -  Fractional  area  of  the  radiation  through  the  top  of  the 
cloud  as  observed  from  the  satellite 
LATM  -  Radiance  of  the  atmosphere 
SATGT  -  Solid  angle  of  the  target 


2. 10  Obstructions  And  Mountain  Modelling 

Obstructions  have  been  included  in  the  model.  These 
obstructions  allow  the  computation  of  daily  radiation  by 
modifying  the  beginning  and  ending  times  that  direct 
radiation  may  fall  on  the  inclined  plane.  The  obstructions 
do  not  modify  the  length  or  amount  of  diffuse  radiation 
time.  From  sunrise  to  sunset,  diffuse  radiation  is  received 
even  though  direct  radiation  from  the  sun  is  blocked.  These 
obstructions  were  designed  to  estimate  the  radiation 
reaching  the  surface  in  a  mountainous  area  or  an  area  with 
buildings.  These  calculations  will  likely  be  underestimates 
in  every  case.  The  model  does  not  account  for  the  effects  of 
multiple  reflections  due  to  nearby  inclined  surfaces  which 
contribute  to  the  total  radiation  received  at  a  particular 
location.  The  following  assumptions  are  made  when  including 
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the  obstructions: 


1.  The  obstructions  are  assumed  to  be  far  enough 
from  the  inclined  plane  that  the  radiation 
reflected  from  the  obstructions  can  be 
considered  negligible. 

2.  The  obstructions  exactly  block  sunrise  and 
sunset  until  the  desired  elevation  angle  is 
reached. 

3.  The  obstruction  is  assumed  to  block  an 
insignificant  portion  of  the  diffuse  sky 
radiation  regardless  of  the  elevation  angle  to 
the  top  of  the  obstruction. 

4.  The  inclined  plane  is  resting  on  a  horizontal 
surface  between  either  one  or  two  obstructions 
and  does  not  receive  radiation  reflected  from 
any  other  nearby  inclined  object. 


The  obstructions  are  defined  by  the  elevation  angles  to 
the  top  of  the  obstruction  as  measured  from  a  horizontal 
surface  at  the  location  of  the  inclined  plane.  While  the 
usefulness  of  this  type  of  obstruction  is  very  limited,  it 
provides  a  starting  point  for  some  preliminary  studies  of 
how  obstructions  affect  the  radiation  reaching  an  inclined 
plane  surface. 


CHAPTER  III 

ANALYSIS  OF  MODEL  RESULTS 

t. 

In  this  chapter,  the  discussion  will  focus  on  the  model 
results.  The  analysis  of  the  model  will  contain  comparisons 
with  data  when  available.  The  specific  goals  are  to: 

1.  Compare  the  model  results  with  observations 

2.  Assess  the  effects  of  spectral  surface  albedos 

3.  Assess  the  effects  of  cloud  type,  thickness,  and 
amount 

4.  Assess  the  effects  of  solar  zenith  angle  (which 
is  directly  affected  by  latitude,  and  hour 
angle ) 


The  effects  of  turbidity,  water  vapor  amount,  air  mass  path, 
solar  zenith  angle,  and  ozone  amount  for  clear  skies  have 
been  published  by  others  (Bird  and  Riordan  1984;  Hatfield, 
Giorgis,  Flocchini,  1981).  Since  the  intent  of  this  work  is 
to  concentrate  on  the  effects  of  surface  albedo,  cloud  type 
and  amount,  cloud  thickness,  solar  zenith  angle,  and  cloud 
albedo,  the  effect  of  the  other  parameters  will  be  discussed 
only  when  directly  pertinent  to  interpreting  the  results. 
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3 . 1  Comparison  With  Observation 

The  first  step  is  to  compare  the  model  against  data. 

The  Bird  simple  clear  sky  model  (Bird,  1984)  is  the  basis 
for  the  overall  model  and  was  modified  to  contain  spectral 
aerosol  single  scattering  albedos,  spectral  aerosol  forward 
to  total  scattering  ratios,  spectral  ground  albedos  and 
tilted  surfaces.  To  confirm  that  these  changes  have 
maintained  the  integrity  of  the  original  model,  comparisons 
with  observed  clear  sky  data  were  completed. 

3.1.1  Clear  Sky  Data 

There  are  two  sources  of  data  for  the  clear  sky 
comparisons.  Broadband  data  was  obtained  from  the  Solar  and 
Meteorological  Research  Program  at  the  University  of 
Michigan.  Spectral  data  was  obtained  from  the  Solar  Energy 
Research  Institute  (SERI)  in  Golden,  Colorado. 

The  broadband  data  from  the  University  of  Michigan 
consists  of  data  which  covers  several  spectral  bands.  The 
instrumentation  used  to  collect  this  data  was  the  Eppley 
Precision  Spectral  Pyranometer  (PSP),  the  Eppley 
Pyheliometer ,  and  the,  Eppley  Ultraviolet  (UV)  radiometer. 

The  global  and  diffuse  radiation  on  a  horizontal  surface  was 
measured  by  the  Eppley  PSP's.  The  Eppley  PSP  measuring  the 
diffuse  radiation  was  shaded  by  an  occulting  disk  mounted  on 
a  solar  tracker.  This  method  was  preferred  to  the  shadow 
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band  technique  because  it  was  unnecessary  to  apply  shadow 
band  corrections  and  introduce  uncertainties  associated  with 
those  corrections.  The  alignment  of  the  occulting  disk  is 
very  important  in  correctly  measuring  the  diffuse  radiation 
and  only  those  days  with  proper  alignment  were  used  for 
comparison.  The  occulting  disk  is  5  cm  in  diameter  and  50  cm 
from  the  pyranometer.  This  produces  a  5.7  degree 
f ield-of-view  which  is  blocked  by  the  occulting  disk  and  was 
chosen  so  that  the  blocked  f ield-of-view  matched  the  angular 
f ield-of-view  of  the  Eppley  Pyrheliometer  which  measured  the 
direct-normal  radiation.  Another  Eppley  PSP  was  mounted  on  a 
south-facing  surface  tilted  at  42.3  degrees  (the  latitude  of 
Ann  Arbor,  Michigan).  This  pyranometer  measured  the  total 
solar  radiation  arriving  at  the  tilted  surface  in  the 
spectral  band  0.28  to  2.8  micrometers.  (Department  of 
Atmospheric  and  Oceanic  Science,  University  of  Michigan, 
1979) . 


Summary  Of  Data  From  The  University  Of  Michigan 

Global  radiation  0.28  to  2.8  micrometers 
Diffuse  radiation  0.28  to  2.8  micrometers 
Direct-normal  0.28  to  2.8  micrometers 
Global  Ultraviolet  0.290  to  0.385  micrometers 
Global  with  red  filter  (RG-2)  0.63  to  2.8  micrometers 
Total  radiation  with  42.3  Deg  Tilt  0.28  to  2.8  micrometers 


The  spectral  data  was  obtained  from  the  Solar  Energy 
Research  Institute  (SERI).  "Spectral  irradiance  measurements 
were  taken  with  SERI-modif ied  LI -COR  model  LI-1800  portable 
spectroradiometers .  These  units  were  modified  by 
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incorporating  integrating  spheres  (used  in  measurements  of 
global  irradiance),  direct  beam  modules  (used  to  measure 
direct  normal  irradiance  in  a  5  degree  f ield-of-view,  and 
temperature  controllers.  The  instruments  are  calibrated  and 
operated  at  40  degrees  Celsius  and  are  shielded  (except 
aperture)  from  direct  sunlight  to  prevent  heating  above  the 
control  temperature.  The  optical  system  of  the  LI-1800 
consists  of  a  filter  wheel,  holographic  grating 
monochromator,  and  a  silicon  detector.  The  wavelength  range 
is  300-1100  nm  with  a  bandwidth  of  6  nm  and  a  step  size  of 
1,  2,  5,  or  10  nm,  which  is  selected  by  the  user.  The  total 
scan  time  for  the  300  to  1100  nm  at  2  nm  steps  is  about  27 
seconds.  A  portable  terminal  and  an  internal  microprocessor 
are  used  to  start  the  scans  and  to  perform  data  storage, 
reduction,  and  retrieval  functions."  (Riordan,  1986). 

Summary  Of  Data  From  SERI 

Global  spectral  300  to  1100  nanometers 

(0.3  to  1.1  micrometers) 

Direct-normal  spectral  300  to  1100  nanometers 

(0.3  to  1.1  micrometers) 

Global  spectral  with  37.0  deg  tilt  300  to  1100  nanometers 

(0.3  to  1.1  micrometers) 

3.1.2  Comparison  Wi th  Spectral  Clear  Sky  Observations 

The  data  obtained  from  Dr.  Riordan  at  SERI  consisted  of 
observations  made  on  August  7,  1985  (Day  219)  at  Golden, 
Colorado  (latitude  39.75  N,  longitude  105.16  W).  Spectral 
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observations  of  the  global  total  radiation  on  a  surface 
tilted  at  37  degrees  zenith  angle  were  made.  The  tilted 
surface  was  aligned  so  that  it  was  sun-facing  at  the  time 
each  measurement  was  taken.  At  the  same  time,  the 
corresponding  direct-normal  radiation  was  measured.  Spectral 
observations  of  the  global  radiation  on  a  horizontal  surface 
and  the  corresponding  direct-normal  radiation  were  also 
made.  From  the  10  sets  of  observations  obtained  from  Dr. 
Riordan,  only  one  was  for  a  horizontal  surface.  Three  sets 
of  observations  will  be  compared.  They  were  chosen  to  span 
the  range  of  solar  zenith  angle  covered  by  the  observations. 
Table  3.1.2  lists  some  of  the  pertinent  atmospheric 
characteristics  at  the  time  of  the  observations.  The  water 
vapor  value  for  this  day  was  obtained  from  the  National 
Weather  Service  (NWS)  and  the  atmospheric  turbidity  (aerosol 
optical  depth  at  500  nm)  was  calculated  using  a  sun 
photometer  (Riordan,  1986).  The  model  requires  a  surface 
albedo  and  the  surface  albedo  thought  to  be  most 
representative  of  Colorado  was  the  bare  soil  surface. 
Colorado  lacks  heavy  broadleaf  foliage  and  has  predominantly 
pine,  aspen,  and  evergreen  trees.  The  vegetation  is  also 
relatively  sparse;  even  in  the  months  of  high  growth.  These 
conclusions  were  drawn  from  the  personal  experience  of  the 
author  who  spent  four  years  in  Colorado. 

The  spectral  observation  data  for  the  horizontal 
surface  are  plotted  against  the  modelled  spectral  values  and 
are  shown  on  Figure  3. 1.2. A.  The  solid  lines  are  the 
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observed  values  and  the  dotted  lines  are  the  modelled 
values.  The  majority  of  the  disagreement  of  the  direct 
normal  radiation  with  observations  is  that  the  modelled 
values  are  almost  always  lower  than  measurement.  This  is 
expected  since  the  model  computes  the  direct  radiation  while 
the  measurement  contains  both  direct  and  circumsolar  diffuse 
radiation.  This  also  is  in  agreement  with  the  assumptions 
concerning  the  formulation  of  the  transport  equation  used 
for  this  model.  The  form  of  the  transport  equation  selected 
for  use  in  this  model  tends  to  over-attenuate  the  radiation 
reaching  the  surface  of  the  Earth.  After  scanning  the 
modelled  and  observed  values,  the  largest  disagreement  is 
approximately  80  watts  /  sq  m  /Vm  out  of  a  total  observed  of 
nearly  1100  watts  /  sq  m  /fxm  at  the  peak.  This  corresponds 
to  a  maximum  error  of  7.5%.  The  majority  of  the  spectral 
comparisons  are  much  smaller  than  this  maximum  and  fall 
below  the  5%  error  level.  Based  on  the  comparison  plots  of 
the  direct-normal  and  the  global  radiation  with 
measurements,  it  is  reasonable  to  infer  that  the  overall 
accuracy  is  within  5  percent.  While  it's  desirable  to  obtain 
spectral  diffuse  observations  to  compare  with  the  model 
diffuse  values,  spectral  diffuse  observations  are  not 
available  at  this  time. 

The  spectral  observations  for  a  tilted  surface  are 
plotted  against  the  model  values  in  Figures  3.1.2.B  and 
3.I.2.C.  The  close  agreement  between  the  modelled  and 
observed  values  indicate  the  model  properly  accounts  for 
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tilted  surfaces.  Bird's  original  simple  clear  sky  model 
(Bird,  1984)  did  nqt  account  for  tilted  surfaces  or  spectral 
albedos.  Bird's  latest  model  (Bird  and  Riordan,  1984) 
handles  tilted  surfaces  by  using  the  technique  of  Justus  and 
Paris  (1984).  Figure  3.1.2.B  compares  the  data  for  a 
sun-facing  37  degree  tilted  surface  where  the  solar  zenith 
angle  relative  to  the  tilted  surface  is  24.36  degrees. 

Figure  3.1.2.C  compares  the  data  when  the  solar  incidence 
angle  is  12.02  degrees.  In  both  instances,  the  comparison  is 
good.  Based  on  the  comparisons,  the  overall  accuracy  is 
likely  within  the  desired  5%  error  for  the  majority  of 
tilted  surface  orientations.  Note  that  the  global  radiation 
reaching  the  tilted  surface  is  greater  than  the 
direct-normal  radiation  arriving  at  the  same  location.  This 
is  attributable  to  the  sky  and  ground  diffuse  radiation.  The 
ground  diffuse  radiation  is  defined  as  the  radiation 
reflected  from  the  ground  which  reaches  the  collector 
surface.  The  maximum  direct  radiation  a  surface  could 
receive  is  when  the  surface  is  direct-normal  to  the  incoming 
solar  beam.  Since  the  tilted  surface  is  not  oriented  normal 
to  the  incoming  sunlight,  the  direct  radiation  must  be  less 
than  the  direct-normal  radiation.  This  means  the  addition  of 
the  sky  and  ground  diffuse  radiation  to  the  direct  radiation 
causes  the  global  total  radiation  to  be  greater  than  the 
direct-normal  radiation.  This  contrasts  with  a  horizontal 
surface  and  clear  sky  conditions  where  the  global  radiation 
(direct  plus  sky  diffuse  only)  is  usually  less  than  the 
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tilted  surfaces.  Bird's  original  simple  clear  sky  model 
(Bird,  1984)  did  nqt  account  for  tilted  surfaces  or  spectral 
albedos.  Bird's  latest  model  (Bird  and  Riordan,  1984) 
handles  tilted  surfaces  by  using  the  technique  of  Justus  and 
Paris  (1984).  Figure  3.1.2.B  compares  the  data  for  a 
sun-facing  37  degree  tilted  surface  where  the  solar  zenith 
angle  relative  to  the  tilted  surface  is  24.36  degrees. 

Figure  3.1.2.C  compares  the  data  when  the  solar  incidence 
angle  is  12.02  degrees.  In  both  instances,  the  comparison  is 
good.  Based  on  the  comparisons,  the  overall  accuracy  is 
likely  within  the  desired  5X  error  for  the  majority  of 
tilted  surface  orientations.  Note  that  the  global  radiation 
reaching  the  tilted  surface  is  greater  than  the 
direct-normal  radiation  arriving  at  the  same  location.  This 
is  attributable  to  the  sky  and  ground  diffuse  radiation.  The 
ground  diffuse  radiation  is  defined  as  the  radiation 
reflected  from  the  ground  which  reaches  the  collector 
surface.  The  maximum  direct  radiation  a  surface  could 
receive  is  when  the  surface  is  direct-normal  to  the  incoming 
solar  beam.  Since  the  tilted  surface  is  not  oriented  normal 
to  the  incoming  sunlight,  the  direct  radiation  must  be  less 
than  the  direct-normal  radiation.  This  means  the  addition  of 
the  sky  and  ground  diffuse  radiation  to  the  direct  radiation 
causes  the  global  total  radiation  to  be  greater  than  the 
direct-normal  radiation.  This  contrasts  with  a  horizontal 
surface  and  clear  sky  conditions  where  the  global  radiation 
(direct  plus  sky  diffuse  only)  is  usually  less  than  the 


direct-normal  radiation. 

3.1.3  Comparison  With  Broadband  Clear  Sky  Data 

The  model  can  determine  the  radiation  in  any  spectral 
band  between  0.28  and  4.0  micrometers.  In  addition,  the 
model  can  accumulate  the  amount  of  radiation  over  a  time 
interval  of  less  than  one  hour  or  over  the  entire  day.  Daily 
values  of  broadband  radiation  were  tabulated  from  broadband 
observations  made  at  the  University  of  Michigan  for  the 
Solar  and  Meteorological  Research  Program.  To  test  the  model 
against  observed  daily  values,  it  was  necessary  to  choose 
days  where  the  sky  was  clear  and  the  variation  in  the 
turbidity  and  water  vapor  was  reasonably  minimized.  Ed 
Ryznar,  a  researcher  at  the  University  of  Michigan,  had  for 
his  own  studies  selected  a  set  of  clear  days  from  the 
observations  (Ryznar,  1981).  A  clear  day  is  defined  as  a 
cloudless  day.  From  this  data  set,  one  day  each  month  was 
selected  as  representative  of  a  clear  sky  day.  Days  were 
selected  where  all  (or  all  but  one)  of  the  instruments  were 
aligned,  calibrated,  and  working.  The  values  of  turbidity 
(at  500  nm)  and  water  vapor  amount  were  made  using  a  sun 
photometer.  The  sun  photometer  measurements  were  made 
irregularly  during  the  day.  Some  days  had  more  measurements 
than  others.  Based  on  the  measurements  available  for  each 
day,  a  typical  value  of  turbidity  and  a  typical  value  of  the 
water  vapor  amount  was  chosen  as  representative  of  that  day. 
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The  largest  variation  in  the  water  vapor  amount  over  a  day 
was  0.27  cm  and  the  largest  variation  in  the  turbidity  (at 
500  nm)  over  a  day  was  0.045.  Table  3. 1.3. A  contains  the 
atmospheric  parameters  observed  for  clear  sky  days.  The 
model  calculated  the  daily  amount  of  radiation  based  on 
these  typical  atmospheric  values  which  were  assumed  constant 
throughout  the  day. 

The  spectral  surface  albedo  representing  bare  soil  was 
used  in  the  model  computations.  The  bare  soil  albedo  was 
chosen  for  comparison  purposes  because  it  was  thought  to 
best  represent  Ann  Arbor  throughout  the  year.  Using  the  same 
spectral  albedo  throughout  the  year  provides  a  baseline  for 
comparison  purposes.  It  may  be  possible  to  obtain  closer 
agreement  between  the  observed  and  modelled  values  by 
changing  the  surface  albedo  on  a  monthly  or  seasonal  basis. 
However,  performing  such  an  analysis  would  perhaps  create 
the  impression  of  having  selectively  chosen  the  surface 
albedo  to  obtain  the  best  results.  From  the  author's  point 
of  view,  this  is  not  desirable. 

In  late  1982,  a  volcano  (El  Chichon)  erupted  in  Mexico 
and  spewed  ash  and  other  constituents  into  the  atmosphere. 
The  volcanic  cloud  associated  with  this  eruption  reached  Ann 
Arbor,  Michigan  during  October  1982.  For  comparison  with  the 
clear  sky  days  chosen  previously,  a  set  of  clear  sky  days 
during  the  year  following  the  volcanic  eruption  was  also 
selected.  These  days  were  chosen  using  the  same  criteria  as 
the  previous  set  of  clear  sky  days.  While  there  is  no 
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absolute  guarantee,  it  is  hoped  that  the  primary  difference 
between  the  two  data  sets  will  be  the  aerosol  distribution. 
As  before,  one  clear  sky  day  each  month  was  chosen.  Table 

3.1.3. B  contains  the  atmospheric  parameters  used  for  the  El 
Chichon  clear  sky  days.  The  modelled  daily  values  were  based 
on  these  typical  atmospheric  values  which  were  assumed 
constant  throughout  the  day. 

The  observed  versus  the  modelled  values  of  the  global, 
direct-normal,  direct,  diffuse,  UV,  RG2 ,  and  tilted  surface 
(42.3  degrees)  have  been  plotted  for  the  normal  clear  sky 
days  and  the  El  Chichon  clear  sky  days.  Figures  3.1.3.C  and 
D  show  the  normal  clear  sky  comparison  plots.  Figures 

3.1.3. E  and  F  show  the  El  Chichon  clear  sky  comparison 
plots.  Note  that  the  scales  have  been  varied  from  plot  to 
plot  to  highlight  the  differences  between  the  observed  and 
modelled  values.  Tables  3.1.3.G  and  3.1.3.H  contain  a  simple 
error  analysis  of  the  normal  clear  sky  and  the  El  Chichon 
clear  sky  cases  respectively. 

For  the  typical  clear  sky  cases,  the  broadband  global 
and  total  measurements  (global,  RG2 ,  and  tilted)  agree  with 
the  modelled  daily  values.  The  other  measurements 
(direct-normal,  direct,  diffuse,  and  UV)  do  not  agree  as 
well.  The  direct-normal  and  direct  modelled  values  appear  to 
consistently  underestimate  the  observed  values;  while  the 
modelled  diffuse  values  appear  to  consistently  overestimate 
the  observed  diffuse.  In  the  worst  case,  this  overestimate 
is  nearly  100  percent  of  the  observed  daily  values.  However, 
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MONTH 

YEAR 

DAY 

OF  THE 
YEAR 

DAY 

OF  THE 
MONTH 

TURBIDITY 
( 500  nm) 

WATER 
VAPOR 
( cm) 

JAN 

81 

18 

18 

0.1 

1  .0 

FEB 

82 

34 

4 

0.16 

0.19 

MAR 

82 

72 

13 

0.2 

0.3 

APR 

81 

105 

15 

0.13 

0.38 

MAY 

81 

140 

20 

0.305 

0.80 

JUN 

81 

168 

17 

0.25 

1.0 

JUL 

80 

188 

6 

0.3 

2.20 

AUG 

81 

231 

19 

0.2 

1 .0 

SEP 

80 

262 

18 

0.2 

2.20 

OCT 

80 

284 

10 

0.14 

1.3 

NOV 

81 

315 

11 

0.14 

0.79 

DEC 

81 

344 

10 

0.12 

0.3 

Table  3. 1.3. A 

Atmospheric  Parameters  For  The  Clear  Sky  Days 
At  Ann  Arbor,  Michigan 
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MONTH 

YEAR 

DAY 

OF  THE 
YEAR 

DAY 

OF  THE 
MONTH 

TURBIDITY 
( 500  nm) 

WATER 

VAPOR 

(cm) 

OCT 

82 

299 

26 

0.31 

1.14 

NOV 

82 

319 

15 

0.17 

0.55 

JAN 

83 

1 

1 

0.28 

0.51 

JAN 

83 

18 

18 

0.22 

0.14 

FEB 

83 

44 

14 

0.41 

0.97 

MAR 

83 

71 

12 

0.18 

0.56 

APR 

83 

110 

20 

0.30 

0.50 

MAY 

83 

129 

9 

0.23 

0.40 

JUN 

83 

173 

22 

0.263 

2.31 

JUL 

83 

194 

13 

0.168 

1.38 

AUG 

83 

236 

24 

0.241 

2.73 

OCT 

83 

282 

9 

0.125 

0.88 

Table  3.1.3.B 


Atmospheric  Parameters  For  The  Clear  Sky  Days  With 
The  El  Chicon  Volcanic  Cloud  At  Ann  Arbor,  Michigan 
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Figure  3.1.3.C 
Daily  Irradiance  Comparison 

The  daily  observed  irradiance  is  compared  with  the  daily  modelled 
irradiance.  The  spectral  band  is  0.28  to  2.8  micrometers. 
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Figure  3.1.3.D 
Daily  Irradiance  Comparison 

The  daily  observed  irradiance  is  compared  with  the  daily  modelled 
irradiance.  The  spectral  band  is  0.28  to  2.8  micormeters . 
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Figure  3.1.3.E 

Daily  Irradiance  Comparison 
With  The  El  Chichon  Volcanic  Cloud 

The  daily  observed  irradiance  is  compared  with  the  daily  modelled 
irradiance.  The  spectral  band  is  0.28  to  2.8  micrometers. 
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Figure  3.1.3.F 

Daily  Irradiance  Comparison 
With  The  El  Chichon  Volcanic  Cloud 

The  daily  observed  irradiance  is  compared  with  the  daily  modelled 
irradiance.  The  spectral  band  is  0.28  to  2.8  micrometers. 
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RADIATION 

TYPE 

SMALLEST 
DIFFERENCE 
(MJ  /  DAY) 

ERROR 

(7.) 

LARGEST 
DIFFERENCE 
(MJ  /  DAY) 

ERROR 

(7.) 

Global 

0.044 

0.17 

1.287 

5.10 

Direct-Normal 

1.37 

4.86 

6.253 

15.30 

Direct 

0.38 

5.7 

4.0 

15.5 

Diffuse 

0.5 

29.9 

3.55 

93.2 

Ultraviolet 

0.011 

1.04 

0.1 

33.22 

RG2  (Red  Filter) 

0.0 

0.0 

0.97 

15.67 

Tilted 

0.4 

1.58 

2.845 

13.29 

Table  3.1.3  .G 


Model  Errors  Associated  With  The  Clear  Sky  Days 
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RADIATION 

TYPE 

SMALLEST 
DIFFERENCE 
(MJ  /  DAY) 

ERROR 

(%) 

LARGEST 
DIFFERENCE 
(MJ  /  DAY) 

ERROR 

(%) 

Global 

0.24 

0.8 

1.608 

5.61 

Direct-Normal 

0.03 

0.21 

2.193 

7.022 

Direct 

0.022 

0.097 

0.886 

3.709 

Diffuse 

0.042 

0.67 

1.465 

34.641 

Ultraviolet 

0.003 

0.22 

0.255 

20.383 

RG2  {Red  Filter) 

0.051 

0.32 

1.188 

7.048 

Tilted 

0.476 

1.7 

4.012 

18.945 

Table  3.1.3.H 


Model  Errors  Associated  With  T5ie  El  Chichon  Clear  Sky  Days 
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RADIATION 

TYPE 

ESTIMATED 

PERCENTAGE  ERROR 
(RURAL  AEROSOL) 
DISTRIBUTION 

ESTIMATED 
PERCENTAGE  ERROR 
(UNKNOWN  AEROSOL) 
DISTRIBUTION 

Global 

0  - 

7 

+ 

15 

Direct-Normal 

♦ 

5  - 

10 

+ 

25 

Direct 

♦ 

5  - 

10 

+ 

25 

Diffuse 

+ 

15  - 

20 

+ 

100 

Ultraviolet 

+ 

15  - 

20 

+ 

50 

RG2  (  Red  Filter) 

+ 

5  - 

10 

50 

Tilted 

+ 

10  - 

20 

+ 

50 

Table  3. 1.3. I 

Qualitative  Estimates  of  Model  Error 
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the  overestimate  in  the  modelled  diffuse  values  and  the 
underestimate  in  the  direct  and  direct-normal  radiation 
values  nearly  compensate  for  each  other.  The  result  is 
consistently  more  accurate  broadband  global  and  total 
radiation  values.  This  is  expected,  due  to  the  formulation 
of  the  transport  equation,  where  the  direct  radiation  is 
dependent  on  the  aerosol  and  Rayleigh  transmittances  and  the 
diffuse  radiation  is  dependent  on  the  complement  of  those 
transmittances. 

The  most  influential  factor  in  the  clear  sky  model  is 
the  aerosol  distribution.  In  this  model,  the  aerosol 
distribution  was  assumed  to  resemble  the  rural  aerosol 
distribution.  The  clear  sky  days  probably  have  an  aerosol 
distribution  different  from  the  rural  aerosol  distribution. 
Since  both  the  direct-normal,  direct,  and  diffuse  values 
depend  on  the  aerosol  transmittance,  an  inappropriate 
aerosol  distribution  will  affect  the  three  quantities.  The 
clear  days  probably  have  an  aerosol  distribution  which  has 
fewer  particles  at  every  radius  than  the  rural  aerosol 
distribution.  This  will  mean  a  higher  transmittance  through 
the  atmosphere  and  less  scattering.  This  is  observed  in  the 
plots;  reduced  direct-normal  and  direct  values  and  increased 
diffuse  values.  Fortunately,  the  amount  of  underestimate  in 
the  direct  case  and  the  overestimate  in  the  diffuse  case 
tend  to  offset  each  other  so  that  the  global  and  total 
radiation  values  show  consistently  less  error  when  compared 
with  observation.  The  aerosol  distribution  should  be 
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significantly  different  for  the  El  Chichon  clear  sky  cases 
and  one  should  see  a  definite  change  in  the  comparison  of 
observed  and  modelled  values. 

Reviewing  the  plots  of  the  modelled  versus  observed 
values  for  the  El  Chichon  clear  sky  cases  and  the  typical 
clear  sky  cases,  one  can  see  the  El  Chichon  comparisons 
indicate  closer  agreement  especially  in  the  direct-normal, 
direct,  and  diffuse  quantities.  This  indicates  the  effects 
of  the  volcanic  cloud  over  Ann  Arbor  tends  to  more  closely 
resemble  the  model  calculations  using  the  rural  aerosol 
distribution  than  the  distribution  from  the  clear  sky  days 
without  the  volcanic  cloud.  One  obvious  conclusion  that  can 
be  drawn  is  that  the  aerosol  distribution  is  very  important. 
This  model  should  only  be  used  with  those  cases  which 
closely  approximate  the  rural  aerosol  distribution.  One 
turbidity  value  at  500  nm  is  not  sufficient  to  determine 
whether  the  rural  aerosol  distribution  should  be  used.  Three 
to  five  turbidity  measurements  should  be  obtained  at 
different  wavelengths  and  compared  with  the  turbidities 
generated  using  the  rural  aerosol  model  to  be  more  certain 
whether  the  aerosol  distribution  resembles  the  modelled 
rural  aerosol  distribution.  Since  the  model  has  only  one 
aerosol  distribution,  it  is  limited  in  application  to  those 
situations  which  closely  approximate  a  rural  aerosol 
distribution  (particularly  if  numeric  accuracy  is  required). 
The  model's  usefulness  could  be  expanded  by  having  more  than 
one  aerosol  distribution  available. 
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Because  there  is  not  sufficient  data  available  to 
determine  how  close  either  the  clear  sky  or  the  El  Chichon 
clear  sky  aerosol  distributions  are  to  the  rural  aerosol 
distribution  used  in  the  model,  it  is  impossible  to 
quantitatively  assign  error  to  the  model.  However,  Tables 
3.1.3.G  and  3.1.3.H  show  the  maximum  and  minimum  absolute 
error  for  the  days  used  in  the  plots.  The  smallest  and 
largest  differences  in  megajoules  per  day  were  divided  by 
the  observed  daily  values  to  obtain  the  percentage  error. 

The  tables  indicate  the  errors  between  the  observed  values 
and  the  modelled  values.  Realizing  a  difference  in  the 
aerosol  distributions  affected  the  errors,  one  can  infer 
what  the  range  of  error  for  the  model  might  be  based  on 
whether  the  aerosol  distribution  resembles  the  rural  aerosol 
distribution.  Table  3. 1.3.1  consists  of  inferred  estimates 
of  how  good  the  model  might  typically  be  for  both  a  rural 
aerosol  distribution  and  an  unknown  aerosol  distribution. 

The  error  estimates  are  subjective  but  reflect  the  error  one 
could  expect  when  using  the  model.  The  observed  versus 
modelled  comparisons  only  cover  a  few  of  the  possible 
situations.  The  inferred  error  estimates  are  an  attempt  to 
give  a  typical  error  for  the  model  and  not  an  upper  bound 
for  the  model  error. 

3.1.4  Cloudy  Sky  Data 
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There  are  no  good  cloudy  sky  data  against  which  to  I 

compare  this  model  against.  A  good  data  set  will  be  defined  j 

as  one  where  enough  information  is  available  to  completely  , 

describe  a  particular  set  of  circumstances.  There  are  many 
problems  inherent  in  comparing  this  model  with  observations. 

To  begin  with,  there  are  no  rectangular  clouds  in  the  real  j 

atmosphere.  Based  on  the  study  by  Welch  and  Zdunkowski  ! 

(1980),  suitable  dimensions  for  rectangular  clouds  could  be  | 

used  to  obtain  realistic  calculations  provided  data 
regarding  the  size  and  separation  distances  between  the 
clouds  were  available.  The  model  assumes  a  particular 

geometry  for  the  clouds  (rectangular  array  of  clouds,  i 

striated  array  of  clouds  or  an  overcast).  Very  few 
observations  of  cloud  cover  describe  the  cloud  pattern  in 
the  sky.  The  model  calculates  the  spatial  average  of  the 

radiation  reaching  the  ground  in  the  clear  sky  and  under  the  ‘ 

cloud.  The  instrumentation  measuring  the  radiation  is 

basically  a  point  collector.  The  instrumentation  is  either 

in  the  sunlight  or  in  the  shadow  of  a  cloud  (not  both  at  the 

same  time).  This  presents  a  problem  in  comparing  the 

observations  to  the  modelled  values.  If  one  assumes  that 

over  an  hour  the  instruments  will  pass  in  and  out  of  the 

cloud  shadows,  the  integration  of  the  instrument 

measurements  into  hourly  totals  may  provide  a  data  set  for 

comparison.  However,  the  variability  in  such  data  would 

probably  be  large  since  the  accuracy  of  such  measurements 

would  depend  on  the  size  of  the  clouds  and  their  velocity. 
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Overcast  skies  present  fewer  problems  than  partially  cloudy 
skies.  In  overcast  situations,  the  instrumentation  is  always 
in  the  ’shadow'  of  the  clouds  and  the  spatial  average  under 
the  clouds  is  the  same  value  the  instrumentation  will 
measure.  The  problem  with  the  overcast  situation  is 
determining  the  cloud  thickness.  While  rawinsondes  can  give 
information  about  the  cloud  bases  and  tops  it  is  not  highly 
accurate.  A  rawinsonde  passing  out  of  the  top  of  a  cloud  at 
saturation  will  continue  to  show  high  dewpoints  until  the 
sensor  dries  sufficiently.  By  the  time  this  occurs,  the 
rawinsonde  could  have  travelled  several  tens  of  meters  or 
more.  The  cloud  amounts  reported  in  observations  are  highly 
subjective  and  not  necessarily  numerically  accurate 
assessments  of  the  amount  of  cloud  cover.  The  combination  of 
all  these  problems  casts  doubt  on  whether  the  model  can  be 
checked  for  accuracy.  The  accuracy  of  the  model  will  not  be 
decided  by  the  data  currently  available.  More  complete 
observations  must  be  taken  in  order  to  test  this  and  future 
models  of  this  sort.  However,  in  the  interests  of  verifying 
trends  shown  by  the  model,  several  days  of  cloudy  sky  data 
will  be  evaluated.  There  are  many  unknown  quantities  which 
must  be  estimated  and  given  to  the  model  in  lieu  of 
observation.  The  ability  of  any  model  to  reproduce 
observation  under  these  less  than  ideal  conditions  is  very 
low. 

The  data  was  obtained  from  the  Solar  and  Meteorological 
Research  Program  at  the  University  of  Michigan.  The 


instrumentation  and  measurements  were  exactly  the  same  as 
those  described  in  Section  3.1.1  for  clear  sky  data.  The 
nearest  complete  National  Weather  Service  surface 
observation  site  (Detroit  Metropolitan  Airport)  is 
approximately  15  miles  east.  The  nearest  rawinsonde  station 
is  located  in  Flint,  Michigan  approximately  50  miles  north. 
Some  observations  are  kept  at  the  University  of  Michigan, 
but  the  only  cloud  parameters  kept  are  the  cloud  type  and 
amount.  The  cloud  amount  is  recorded  in  terms  of  scattered, 
broken  and  overcast.  There  is  no  estimate  of  cloud  tops.  As 
with  most  observations,  no  information  is  kept  on  the  size 
or  structure  of  the  observed  cloud  formations.  There  is  also 
no  information  on  where  the  clouds  may  be  located  in  the  sky 
relative  to  the  observation  site.  The  estimates  of  cloud 
amount  and  cloud  thickness  are  very  important  in  getting 
model  results.  The  variation  between  0.1  to  0.4  for 
scattered  clouds  and  0.5  to  0.9  for  broken  clouds  can  cause 
large  differences  in  the  model  results. 

The  cloudy  observations  were  developed  from  composite 
information  from  limited  observations  taken  in  Ann  Arbor 
(NWS),  observations  from  Detroit  Metro  (NWS),  and  rawinsonde 
data  from  Flint,  Michigan  (NWS).  All  the  cloudy  observations 
were  during  the  hour  of  1400  to  1500  (2  p.m.  to  3  p.m.) 
local  time.  The  cloud  amounts  were  matched  between  the 
Detroit  (Metro)  and  Ann  Arbor  observations  in  terms  of  the 
categories  scattered,  broken,  and  overcast.  When  they 
matched,  the  Detroit  observation  was  used  to  obtain  the 


cloud  amount  in  tenths.  The  Detroit  observation  was  also 
used  to  determine  the  cloud  type  and  cloud  bases.  The  Flint 
rawinsonde  was  used  to  determine  the  cloud  bases,  cloud 
tops,  cloud  thickness  and  to  determine  if  more  than  one 
cloud  layer  existed.  Only  those  cases  where  the  Detroit 
observation  and  rawinsonde  data  gave  the  same  cloud  base  was 
used.  Since  this  model  is  for  a  single  layer  of  cloud,  those 
days  on  which  the  rawinsonde  indicated  more  than  one 
possible  layer  of  cloud  were  discarded.  Air  Weather  Service 
Manual  105-124  on  the  use  of  the  Skew  T,  Log  P  thermodynamic 
diagram  in  analysis  and  forecasting  indicates  that  clouds 
are  usually  present  when  the  dewpoint  depression  is  less 
than  6  degrees  Celsius.  While  there  are  several  methods  for 
determining  cloud  tops,  bases,  and  thicknesses,  this  method 
is  simple  and  just  as  reliable  as  the  other  techniques.  The 
cloud  tops  and  thicknesses  were  determined  from  the  Flint 
rawinsonde  using  this  technique.  No  information  is  available 
about  the  dimensions  or  the  location  of  the  clouds.  A  cloud 
size  was  chosen  that  represents  a  typical  cloud  (500  m).  The 
spacing  between  the  clouds  was  varied  to  obtain  the  correct 
cloud  amounts.  The  precipitable  water  was  obtained  from  the 
National  Climatic  Center  for  Detroit,  Michigan  for  the 
desired  days.  The  surface  pressure  was  obtained  from  the 
Detroit  observation.  Table  3. 1.4. A  lists  the  observed 
parameters  for  the  cloudy  observations. 

Turbidity  measurements  were  taken  on  the  days  where  the 
mornings  were  clear  and  the  aftermoons  were  cloudy  or  vice 
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versa.  Choosing  an  appropriate  turbidity  value  for  the  model 
was  accomplished  by  picking  days  where  the  clear  sky 
turbidity  measurements  were  made  in  the  morning  and  the  sky 
became  cloudy  in  the  afternoon.  A  water  vapor  amount  was 
inferred  for  the  afternoon  from  the  results  of  the  morning 
sun  photometer  measurements.  These  water  vapor  amounts 
agreed  with  the  precipitable  water  amounts  obtained  from  the 
National  Climatic  Center.  The  largest  difference  between  the 
National  Climatic  Center  values  and  the  sun  photometer  was 
0.2  cm.  For  consistency,  the  National  Climatic  Center  values 
of  precipitable  water  were  used.  On  days  where  the  turbidity 
measurements  could  not  be  obtained  (overcast  days),  a 
turbidity  value  of  0.1  at  500  nm  was  used.  These  values  are 
enclosed  in  parentheses  in  Table  3. 1.4. A. 

3.1.5  Comparison  With  Broadband  Cloudy  Sky  Observations 

The  observed  versus  modelled  results  are  listed  in 
Table  3. 1.5. A.  The  agreement  between  observation  and  model 
is  poor.  This  was  expected.  Conventional  observations  do  not 
provide  the  detail  necessary  to  accurately  compare  the  model 
with  observations.  The  model  was  designed  for  very  specific 
geometric  cases  (an  array  of  clouds,  striated  clouds  or 
overcast  clouds).  In  the  scattered  and  broken  cases,  there 
is  no  information  about  the  geometry  of  the  cloud  formation. 
There  is  also  no  information  about  the  length  of  time  the 
sensor  was  in  the  cloud  shadow  or  in  the  direct  sunlight. 
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This  information  is  very  important  when  comparing  with 
observation.  The  time  the  sensor  is  in  the  direct  sunlight 
versus  the  cloud  shadow  has  a  large  impact  on  the  measured 
irradiance.  A  significant  portion  of  the  difference  between 
the  observations  and  the  model  calculations  is  probably 
attributable  to  this  factor.  Coupled  with  the  fact  that  the 
cloud  geometries  were  probably  not  the  same  as  those  used  in 
the  model,  there  is  reason  to  expect  substantial  differences 
between  the  observations  and  the  modelled  values. 

To  accurately  compare  a  model  with  measurements,  very 
specialized  and  detailed  observations  must  be  kept. 
Information  on  both  the  time  and  spatial  characteristics  of 
the  cloud  pattern  as  well  as  the  location  of  the  clouds,  the 
cloud  geometry,  and  the  cloud  dimensions  will  be  needed. 

Even  though  the  model  comparison  with  observation  is  poor, 
there  is  sufficient  reason  to  expect  it  to  be  and  this 
indicates  that  present  observations  do  not  sufficiently 
characterize  each  situation  to  allow  accurate  comparison. 

One  piece  of  information  can  be  derived  from  the 
comparison  with  observation,  however.  While  all  of  the 
modelled  values  are  too  low  in  comparison  to  actual 
observation,  the  modelled  values  for  the  overcast  cases  are 
much  too  low.  The  observed  cases  are  for  very  thick  clouds 
(699  and  2584  meters).  This  indicates  a  problem  in  the 
modelling  of  the  clouds.  The  simple  solution  to  the 
radiative  transfer  equation  used  to  represent  the  cloud  is 
not  indicative  of  actual  cloud  characteristics.  Based  on  the 
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comparison,  the  results  indicate  that  the  extinction  of 
radiation  through  the  cloud  is  probably  too  large  as  the 
cloud  becomes  thick.  Since  both  the  cloud  reflection  and 
transmission  properties  are  determined  from  the  solution  to 
the  same  radiative  transfer  equation,  the  solution  to  the 
cloud  albedo  should  exhibit  differences  from  measurement. 

3.1.6  Cloud  Reflectance  Data 

A  limited  amount  of  cloud  albedo  data  is  available  for 
comparison;  Kondratev  presents  cloud  albedo  measurements 
made  from  aircraft  by  Chel'tsov  (Kondratev,  1969, 
translation,  1973).  The  data  is  based  on  measureme  s  taken 
above  an  overcast.  Kondratev  indicates  cloud  albedo  varies 
with  cloud  thickness,  cloud  type,  water  content  of  the 
cloud,  solar  altitude,  and  the  underlying  surface  albedo. 
The  underlying  surface  albedo  can  produce  a  variation  of  30 
to  40  percent  in  the  cloud  albedo  (Kondratev, 1969) . 
Kondratev  indicates  the  higher  the  surface  albedo,  the 
higher  the  cloud  albedo;  and  the  lower  the  surface  albedo, 
the  lower  the  cloud  albedo.  The  model  confirms  the  basic 
trends  indicated  by  Kondratev.  This  model  computes  the 
reflectance  of  the  cloud-atmosphere-ground  combination  as 
well  as  the  cloud  albedo  (based  on  the  characteristics  of 
the  cloud  alone).  The  modelled  cloud  albedo  will  be  smaller 
than  the  modelled  cloud-atmosphere-ground  albedo.  The 
thickness  of  the  cloud  will  determine  the  degree  of 
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influence  of  the  surface  albedo.  This  suggests  a  strong 
connection  between  the  underlying  surface  albedo  and  the 
cloud-atmosphere-ground  albedo  (total  albedo).  Since  no 
specific  data  sets  against  which  to  test  the  model  in  this 
respect  are  available,  it  will  only  be  noted  here  that  the 
trends  produced  by  the  model  are  in  agreement  with  observed 
trends.  Table  3. 1.6. A  from  Kondratev  (1969)  indicates 
changes  in  cloud  albedo  with  different  underlying  surfaces. 
For  comparison,  the  cloud  albedo  as  a  function  of  cloud 
thickness  was  computed  for  three  different  underlying 
surfaces  for  Archangel,  Soviet  Union.  The  model  produced 
cloud-atmosphere-ground  albedos  for  Archangel,  Soviet  Union 
with  the  following  parameters  set: 


Latitude 
Longitude 
Hour  Angle 
Cloud  Type 
Cloud  Amount 
Cloud  Height 
Turbidity  (500  nm) 
Water  Vapor 


64.66  N 
41.00  E 
0.00 

Stratocumulus 
Overcast 
1000.00  m 
0.1 

2.0  cm 


Both  Kondratev' s  data  in  Table  3.  1.6. A  and  the  modelled 
albedos  in  Table  3.1.6.B  exhibit  similar  qualities.  The 
largest  differences  associated  with  cloud  albedos  occurred 
when  the  underlying  surface  had  a  low  albedo.  When  the 
underlying  surface  has  a  high  albedo,  the  distinction 
between  the  cloud  albedos  virtually  disappears  regardless  of 
the  thickness  of  the  cloud.  The  thickness  of  the  cloud  has  a 
larger  effect  on  the  cloud  albedo  for  lower  surface 
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Albedo  (%)  of  clouds  of  lower  strstnm  with  different 
density  for  versons  types  of  underlying  surface 


Underlying  surface 

Breaking 
of  snow, 
tenths 

dense 

Clouds 

inter¬ 

mediate 

translucent 

Clear  sky 

Snow,  9-10  tenths 

0 

86 

83 

79 

79 

1 

82 

76 

70 

70 

1-2 

73 

67 

58 

57 

2-3 

71 

61 

51 

45 

Tundra 

— 

61 

46 

30 

17 

Water 

— 

58 

45 

25 

8 

Table  3 . 1 .6. A 

Albedo  of  Clouds  With  Various  Underlying  Surfaces 


(From  Kondratev,  1969) 
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CLOUD 

THICKNESS 

(M) 

KONDRATEV 

TURBID 

WATER 

BARE 

SOIL 

FRESH 

SNOW 

50 

0.29 

0.53 

0.55 

0.63 

100 

0.37 

0.62 

0.63 

0.67 

200 

0.50 

0.68 

0.68 

0.70 

300 

0.62 

0.70 

0.70 

0.71 

400 

0.71 

0.71 

0.71 

0.72 

500 

0.79 

0.72 

0.72 

0.72 

Table  3.1.6.B 

Comparison  Of  Changes  In  Cloud  Albedo  tflch  Cloud  Thickness 
For  Different  Underlying  Surfaces 
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reflectivities.  It  is  the  combination  of  the  cloud  thickness 
and  the  surface  albedo  which  regulates  the  cloud  albedo 
until  the  cloud  becomes  very  thick.  When  the  cloud  becomes 
very  thick  the  effect  of  the  underlying  surface  on  the  cloud 
albedo  is  very  small.  Figure  3.1.6.C  shows  the  cloud  albedo 
computed  by  the  model  based  on  the  cloud  characteristics 
alone.  Since  the  cloud  is  absorptive,  the  albedo  approaches 
some  limit  but  does  not  reach  100%  regardless  of  how  thick 
the  cloud  becomes.  Of  the  four  low  cloud  types  modelled,  the 
clouds  shown  in  the  figure  represent  the  least  absorbing 
cloud  type  (stratus)  and  the  most  absorbing  cloud  type 
(nimbostratus)  as  determined  from  the  absorption 
coefficients.  Figures  3.1.6.D  gives  an  indication  of  the 
changes  cloud  type  can  make  in  affecting  the 
cloud-atmosphere-ground  albedo  when  the  underlying  surface 
is  changed,  while  the  changes  are  on  the  order  of  a  few 
percent,  a  higher  cloud-atmosphere-ground  albedo  is  obtained 
for  the  most  reflective  cloud  type  when  compared  to  a  less 
reflective  cloud  type  over  the  same  surface  albedo. 

Chel’tsov's  measurements  of  cloud  albedo  and  cloud 
thickness  provide  a  benchmark  for  comparing  the  model 
produced  albedos.  Chel'tsov’s  measurements  were  taken  over 
Archangel  (latitude  64.66  N,  longitude  41.99  E),  Soviet 
Union  for  several  types  of  clouds.  Of  the  categories  of 
clouds  used  by  Chel'tsov,  the  overcast  stratocumulus 
category  was  the  only  one  which  contained  a  single  cloud 
type  directly  analogous  to  one  used  in  this  model.  For  this 
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reason,  only  the  stratocumulus  category  will  be  compared. 
Chel'tsov's  values  of  cloud  albedo  versus  cloud  thickness 
are  shown  in  Figure  3.1.6.E  as  reproduced  from  Kondratev 
(1969).  The  data  indicates  quite  a  bit  of  scatter.  Since 
Chel'tsov's  measurements  were  probably  made  on  different 
days,  the  effects  of  solar  zenith  angle  (solar  altitude)  are 
more  than  likely  exhibited  in  the  data.  Kondratev  does  not 
state  when  Chel'tsov’s  measurements  were  taken.  Differences 
in  cloud  conditions  on  distinct  days  also  affect  the 
measurements.  These  two  facts  may  explain  the  scatter  in  the 
data.  Because  of  Archangel's  latitude,  the  solar  zenith 
angle  can  get  quite  large  for  several  months  of  the  year. 

The  model  calculated  the  cloud-atmosphere-ground 
albedos  at  Archangel,  Soviet  Union  with  the  following 
parameters  set: 


Latitude : 

Longitude: 

Hour  Angle: 

Cloud  Type: 

Cloud  Amount: 

Cloud  Height: 
Turbidity  (500  nm): 
Water  Vapor: 

Surface  Albedo  Type: 


64.66  N 
41.00  E 
0.00 

Stratocumulus 
Overcast 
1000.0  m 
0.1 

2.0  cm 
Bare  Soil 


The  cloud  albedo  is  dependent  on  the  solar  zenith  angle 
since  it  determines  the  path  distance  a  light  ray  must 
travel  through  the  cloud.  Day  166  was  chosen  as  a 
representative  summer  day  and  is  plotted  against  Chel’tsov's 
measurements  in  Figure  3.1.6.F,  The  curve  plotted  and 
labelled  Kondratev  in  Figure  3.1.6.F  is  the  curve  Kondratev 


CLOUD  ALBEDO 


ALBEDO 


Legend 

A  STRATUS 

X  NIMBOSTRATUS 


Figure  3.1.6.C 

Model  Calculation  Of  Cloud  Albedo  Based  On  Cloud  Thickness 
And  Cloud  Type 

The  cloud  albedo  is  based  on  cloud  characteristics  alone 
and  the  cloud  albedo  goes  to  zero  as  the  cloud  thickness 
approaches  zero. 


CLOUD  THICKNESS  CLOUD  THICKNESS 
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CUXJD-ATMOSPHERE-GROUW)  albedo 
NkBOSTRATUS  CLOUD 


CLOUD- ATVKTSPHERF -GROUND  ALBEDO 
STRATUS  CLOUD 


Legend 

A  mow  lurfoci 
X  wit  *and  lurfaci 


Figure  3.1.6.D 

Combined  Cloud-Atmosphere-Ground  Albedo  As  A  Function 
Of  Cloud  Thickness 

Day  174  was  used  for  these  calculations . 
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Figure  3.1.6.E 

Cloud  Thickness  And  Cloud  Albedo 

This  figure  represents  data  taken  from  aircraft  above  an  overcast 
deck  of  clouds  over  Archangel,  Russia  (  lat  64.66  N ,  Ion  41.00  E). 
The  curves  fitted  through  the  data  highlight  the  change  in  albedo 
with  thickness.  (Kondratev,  1969). 


CLOUD  THICKNESS 
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CLOUD  REFLECTANCE 


0.25  0.45  0.65  0.85 


REFLECTANCE 


Figure  3.1.6.F 

Comparison  Of  The  Observed  And  Modelled  Relationship  Between 
Cloud  Thickness  And  Cloud  Albedo 

This  figure  contains  the  Kondratev  curve  for  stratocumulus  type 
cloud  and  the  modelled  curve  for  stratocumulus  type  cloud.  The 
distinction  between  the  curves  indicates  how  the  simplified  solu¬ 
tion  for  an  abbreviated  equation  of  radiative  transfer  differs 
from  observation. 


CLOUD  THICKNESS 


CLOUD  REFLECTANCE 


Legend 

A  DAY  166  H  DAY  319 

X  DAY_258_  H  DAY  349 

□  DAY  288 


Figure  3.1.6.G 


The  Effects  Of  Solar  Zenith  Angle  And  Cloud  Thickness 
On  Cloud  Albedo 

Five  curves  are  plotted  for  cloud  thickness  versus  cloud  albedo 
for  different  days  of  the  year.  Each  day  represents  a  different 
solar  zenith  angle.  The  curves  are  representative  of  noon  local 
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fitted  to  Chel'tsov's  data  for  stratocumulus  cloud.  The 
difference  in  the  shape  of  the  modelled  and  observed  cloud 
albedo  curves  indicates  that  the  calculated  albedos  can  be 
improved. 

The  effects  of  solar  zenith  angle  on  cloud  albedo  can 
be  demonstrated.  The  exact  same  scenario  was  modelled  except 
the  day  of  the  year  was  varied.  During  approximately  9 
months  of  the  year  at  Archangel,  the  effect  of  the  solar 
zenith  angle  is  small.  Figure  3.1.6.G  demonstrates  the 
effects  of  solar  zenith  angle  between  days  166  and  349. 
However,  the  solar  zenith  angle  changes  most  rapidly  during 
the  3  winter  months  and  the  results  of  its  effect  on  cloud 
albedo  is  shown. 

Several  conclusions  can  be  drawn  from  the  above 
comparison.  The  solution  used  by  the  model  to  compute  the 
cloud  transmittances  and  reflectances  exhibits  more 
exponential  quality  than  the  observations.  Even  so,  the 
basic  effects  caused  by  variations  in  cloud  thickness,  solar 
zenith  angle,  cloud  type,  and  surface  albedo  seem  to 
preserved. 

3 . 2  The  Effects  Of  Surface  Albedo  With  A  Clear  Sky 

The  effects  of  surface  albedo  on  the  radiation  received 
by  a  plane  at  the  Earth's  surface  can  be  evaluated  with  this 
model.  The  effects  of  changing  the  spectral  surface  albedo 
will  vary  in  each  part  of  the  spectrum  and  it  becomes  nearly 
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impossible  to  quantify  the  error.  For  analysis  purposes,  the 
two  spectral  albedos  available  in  the  model  which  tend  to 
represent  the  extremes  were  compared.  For  this  and  the 
following  sections,  the  degree  of  change  will  be  stated  in 
terms  of  the  largest  differences  noted. 

3.2.1  Daily  Effects  Of  Surface  Albedo  With  A  Clear  Sky 

The  effects  of  surface  albedo  are  expected  to  be  the 
smallest  for  a  horizontal  surface  and  larger  for  tilted 
surfaces.  The  changes  should  also  be  the  smallest  for  a 
clear  sky.  In  a  clear  sky,  there  is  no  multiple  scattering 
between  the  clouds  and  the  ground,  while  there  is  some 
multiple  scattering  between  the  ground  and  the  atmosphere, 
the  magnitude  of  it  does  not  compare  with  the  changes  in 
multiple  scattering  between  cloud  and  ground.  To  verify 
these  common  sense  deductions,  the  set  of  clear  sky  days  for 
Ann  Arbor  without  the  volcanic  cloud  was  modelled  using  the 
exact  same  parameters  except  for  the  surface  albedo.  Instead 
of  using  the  bare  soil  surface  albedo,  the  snow  and  the  wet 
sand  albedos  were  used.  These  two  surfaces  were  chosen 
because  they  represent  the  largest  difference  between  the 
seven  surface  albedo  types  in  the  model.  The  turbid  water 
surface  is  the  lowest  surface  albedo  of  all  the  albedo  types 
in  the  model.  However,  the  water  surface  was  not  used 
because  it  violates  the  assumption  that  the  surface  is 
Lambertian.  The  differences  observed  should  represent  large 
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variations  due  to  surface  albedo  effects.  Figure  3.2.1. A 
shows  the  comparison  plots  of  the  modelled  daily  values  of 
the  global,  diffuse,  UV,  RG2,  and  tilted  surface  radiation. 

The  largest  variation  with  surface  albedo  occurs  during 
the  summer  months  and  the  smallest  during  the  winter  months. 
As  the  solar  zenith  angle  becomes  large,  the  surface  albedo 
increases.  Paltridge  and  Platt's  formulation  for  this  change 
is  nearly  exponential  as  the  solar  zenith  angle  increases 
and  the  differences  between  surface  types  will  begin  to 
disappear  as  the  solar  zenith  angle  becomes  large.  This 
means  the  differences  in  surface  albedo  will  be  most 
pronounced  at  low  solar  zenith  angles  and  least  pronounced 
at  high  solar  zenith  angles;  an  effect  demonstrated  in  every 
comparison  plot. 

The  maximum  amount  of  variation  due  to  changing  surface 
albedo  has  been  estimated  by  taking  the  largest  difference 
between  the  snow  albedo  case  and  the  wet  sand  case  and 
dividing  it  by  the  wet  sand  value.  The  variation  due  to  a 
change  in  surface  albedo  for  a  horizontal  surface  at 
mid-latitude  (Ann  Arbor)  is  approximately: 

Global  6.13% 

Diffuse  23.63% 

UV  19.14% 

RG2  2.17% 


For  a  surface  tilted  at  42.3  degrees,  the  variation  in  the 
global  total  radiation  is  23.21%.  While  the  percentage  for  a 
tilted  surface  is  nearly  the  same  as  for  a  diffuse 
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Figure  3. 2.1. A 


Comparison  Of  Daily  Irradiance  Under  Clear  Skies 
With  Snow  And  Wet  Sand  Surfaces 
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horizontal  surface,  the  amount  of  radiation  accumulated  by 
the  tilted  surface  is  much  greater.  The  direct  radiation 
arriving  at  the  south-facing  tilted  surface  is  greater  than 
that  arriving  at  a  horizontal  surface  due  to  the  plane 
orientation.  However,  the  differences  observed  are  due  to 
the  changes  in  the  sky  and  ground  diffuse  radiation.  The 
effects  of  low  solar  zenith  angle  can  be  seen  on  the  tilted 
surface  just  as  for  the  horizontal  surface.  The  differences 
are  greatest  in  Ann  Arbor's  summer  when  the  solar  zenith 
angle  is  smallest.  The  double  peak  in  the  tilted  surface 
plot  is  due  to  the  orientation  of  the  surface  relative  to 
the  solar  position  during  the  year.  The  shape  of  the  curve 
can  be  changed  by  adjusting  the  orientation  of  the  plane 
surface.  The  effects  of  different  surface  orientations  have 
been  presented  by  Kondratev  (1977)  and  Robinson  (1966). 

3.2.2  Spectral  Effects  Of  Surface  Albedo  With  A  Clear  Sky 

while  the  broadband  variations  of  surface  albedo  are 
interesting,  they  do  not  provide  the  detail  that  a  spectral 
example  can  show.  For  the  sake  of  simplicity,  a  day  near  the 
equinox  (day  74)  was  chosen  for  this  analysis;  a  day  when 
the  solar  zenith  angle  is  not  particularly  high  or  low.  To 
compare  the  cases,  the  following  parameters  were  held 


constant : 


SPECTRAL  IRRADIANCE  SPECTRAL  IRRADiANCI 

SNOW  SURFACE  «  WET  SAND  SURFACE 
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Comparison  Of  The  Global  And  Diffuse  Spectral  Irradiances 
Under  Clear  Skies  With  Snow  and  Wet  Sand  Surfaces 
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Turbidity  (500  nm) 
Water  Vapor 
Surface  pressure 
Hour  angle 
Ozone  amount 


0.100 
2.000  cm 
1013.000  mb 
0.000  (noon) 
0.364  cm 


The  cases  were  for  Ann  Arbor,  Michigan  (latitude  42.2833  N, 
longitude  83.7333  W) .  Figure  3. 2. 2. A  highlights  the 
differences  for  horizontal  snow  and  wet  sand  surfaces.  The 
variation  in  the  diffuse  radiation  has  caused  the 
differences  in  the  figures.  The  direct  radiation  reaching 
the  horizontal  surface  is  the  same  in  each  case.  The 
broadband  difference  (0.28  to  4.0  micrometers)  represented 
by  these  plots  is  47.094  watts  per  square  meter.  This 
represents  a  difference  of  6.37%  in  the  global  radiation  and 
a  change  of  45.77%  in  the  diffuse  radiation.  In  this  case 
(and  typically  in  most  other  cases),  the  largest  changes 
tend  to  occur  between  0.4  and  0.7  micrometers,  an 
atmospheric  window  where  absorption  by  the  atmosphere  is 
small.  Throughout  the  spectrum,  atmospheric  absorption 
causes  the  effects  of  multiple  scattering  to  be  enhanced 
«" ore  in  some  regions  than  others. 


3 . 3  Broadband  Effects  Of  Clouds  And  Surface  Albedo 
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The  modelling  of  clouds  and  surface  albedo  present  many 
problems.  One  of  which  is  trying  to  account  for  the  many 
parameters  which  influence  the  results.  To  analyze  the 
effects  of  one  or  two  parameters,  the  other  parameters  need 
to  be  held  constant  so  the  effect  of  varying  one  variable 
can  be  determined.  For  this  reason,  the  cases  in  this 
section  have  the  following  variables  fixed: 


Latitude 

42.2833  N 

Deg 

Longitude 

83.7333  W 

Deg 

Turbidity  (500nm) 

0.1 

Water  Vapor 

2.0 

cm 

Ozone 

0.3620 

cm 

Day 

174 

Hour  Angle 

0.0 

Noon 

Sky  Condition 

Overcast 

The  intent  of  the  following  sections  is  to  demonstrate  the 
effects  of  cloud  cover  on  the  radiation  reaching  a  plane  at 
the  Earth's  surface.  The  solar  zenith  angle  will  be  the  same 
in  each  case. 

3.3.1  Cloud  Thickness ,  Cloud  Type ,  and  Surface  Albedo 

The  purpose  of  this  section  is  to  demonstrate  the 
qualitative  effects  of  cloud  thickness,  cloud  type,  and 
surface  albedo  on  the  radiation  reaching  a  plane  at  the 
surface  of  the  Earth.  The  cloud  types  and  surface  albedos 
have  been  chosen  so  as  to  obtain  the  maximum  variation.  Of 
the  cloud  types  available  in  this  model,  stratus  is  the 
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least  absorptive  and  nimbostratus  is  the  most  absorptive.  Of 
the  surface  albedo  types,  snow  is  the  most  reflective  and 
wet  sand  the  least  reflective  .  Combinations  of  these  cloud 
and  surface  albedo  types  will  span  the  range  of  variation 
due  to  each  factor.  Figure  3.3.1. A  shows  the  same  situation 
with  a  snow  surface.  Figure  3.3.1.B  shows  the  effect  of 
varying  the  cloud  thickness  over  a  snow  surface  for  both  a 
stratus  and  nimbostratus  overcast. 

The  change  in  cloud  type  can  reduce  the  radiation  for  a 
given  thickness  by  up  to  50  percent.  For  global  radiation,  a 
change  in  cloud  thickness  can  reduce  the  radiation  by  up  to 
100  percent  for  thick  clouds.  The  diffuse  radiation  is  shown 
to  first  increase  for  thin  clouds  and  later  decrease  for 
thicker  clouds.  The  effect  is  caused  by  the  combination  of 
the  scattering  and  absorbing  properties  of  the  cloud.  As  the 
cloud  becomes  thicker,  the  absorption  becomes  large  enough 
to  compensate  for  the  increase  in  scattering  and  the  diffuse 
radiation  decreases.  When  the  cloud  is  thin,  the  absorption 
does  not  compensate  for  the  increase  in  scattering  by  the 
cloud  and  the  diffuse  radiation  increases.  The  maxima  in  the 
diffuse  radiation  occur  when  the  clouds  are  between  25  and 
50  meters  thick.  The  direct  radiation  plotted  as  a  function 
of  thickness  is  identical  in  both  figures.  The  global  and 
diffuse  radiation  are  affected  by  varying  the  cloud  type, 
cloud  thickness,  and  surface  albedo.  The  snow  or  high 
surface  albedo  causes  more  radiation  to  be  received  at  the 
surface  than  the  wet  sand  or  low  albedo  surface  regardless 
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of  cloud  type  or  cloud  thickness. 

Based  on  Figures  3.3.1. A  and  B,  the  following 
conclusions  can  be  drawn.  Of  the  three  parameters  (cloud 
thickness,  cloud  type,  and  surface  albedo),  the  cloud 
thickness  has  the  greatest  effect  on  the  radiation  reaching 
the  surface  of  the  Earth.  Surface  albedo  produces  the  second 
largest  effect  and  cloud  type  produces  the  smallest  effect 
of  the  three  parameters. 

3.3.2.  Cloud  Amount ,  Cloud  Type ,  And  Surface  Albedo 

For  the  comparisons  in  this  section,  the  same  parameter 
values  were  used  as  in  the  last  section.  To  quantify  the 
variation  in  cloud  amount,  the  cloud  thickness  was  constant 
at  100  meters.  The  clouds  were  assumed  to  be  in  a 
rectangular  array  similar  to  a  typical  field  of  cumulus 
clouds.  The  cloud  bases  were  chosen  to  be  100  meters  on  each 
side.  The  distances  between  the  clouds  a  were  adjusted  to 
produce  the  cloud  coverage.  The  radiation  reaching  the 
surface  is  the  average  amount  of  radiation  from  the  clear 
sky,  through  the  side  of  the  rectangular  cloud,  and  through 
the  top  of  the  rectangular  cloud.  In  other  words,  it 
represents  the  spatial  average  of  radiation  reaching  the 


surface.  Figures  3. 3. 2. A  and  B  represent  the  variation  in 
radiation  reaching  the  surface  of  the  Earth  due  to  changing 
cloud  amount,  cloud  type,  and  surface  albedo. 
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By  comparing  Figures  3. 3. 2. A  and  B,  one  can  determine 
that  cloud  amount  has  the  largest  influence  on  the  radiation 
reaching  the  surface  of  the  Earth.  In  comparison,  surface 
albedo  has  the  second  largest  effect;  while  cloud  type  has 
the  smallest  effect. 

3.3.3  Cloud  Thickness ,  Cloud  Amount ,  Surface  Albedo,  and 
Cloud  Type 

By  comparing  the  figures  in  sections  3.3.1  and  3.3.2, 
one  can  quantify  which  of  the  four  parameters  has  the 
largest  effect  on  the  radiation  reaching  the  surface  of  the 
Earth.  Qualitatively,  it  is  straight-forward  to  determine 
from  the  figures  that  cloud  amount  plays  the  dominant  role. 
The  other  parameters  in  descending  order  of  importance  are 
cloud  thickness,  surface  albedo,  and  cloud  type.  These 
conclusions  were  drawn  based  on  the  reduction  in  the 
radiation  reaching  the  surface  of  the  Earth.  It  is  possible 
to  generate  cases  where  the  cloud  thickness  will  be  the 
dominant  factor  instead  of  the  cloud  amount.  However,  in 
most  instances  the  order  of  influence  of  the  parameters  will 
be  as  listed. 

3 . 4  Spectral  Effects  Of  Clouds  And  Surface  Albedo 

The  spectral  effects  of  varying  cloud  thickness,  cloud 
amount,  surface  albedo,  and  cloud  type  will  be  discussed. 
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Many  parameters  have  been  held  constant  so  that  the  effects 
of  cloud  thickness,  cloud  amount,  surface  albedo,  and  cloud 
type  can  be  demostrated.  For  these  cases,  the  following 
variables  were  constant; 


Latitude 
Long  i  tude 

Turbidity  (500  nm) 
Water  Vapor 
Ozone 
Day 

Hour  Angle 


42.2833  N  Deg 
83.7333  W  Deg 
0.1 

2.0  cm 
0.3639  cm 
74 

0.0  Moon 


In  each  of  the  following  sections,  a  different  parameter 
will  be  evaluated  in  combination  with  the  surface  albedo. 
This  will  give  some  insight  into  the  interrelationship 
between  cloud  parameters  and  the  surface  albedo. 


3.4.1  Spectral  Effects  Of  Cloud  Thickness  And  Surface  Albedo 


Figures  3. 4.1. A  through  D  illustrate  the  effects  of 
varying  the  cloud  thickness  and  changing  the  surface  albedo. 
Figures  3.4.1. A  and  B  show  the  effects  of  varying  the 
thickness  of  overcast  nimbostratus  cloud  with  a  snow 
surface.  Figures  3.4. 1.C  and  D  demonstrate  the  same 
situation  with  a  wet  sand  surface.  Figures  3.4. 1.B  and  D 
indicate  that  once  the  cloud  thickness  becomes  large,  the 
differences  observed  are  due  to  changes  in  the  diffuse 
radiation  alone.  This  is  because  there  is  so  little  direct 
radiation  coming  through  a  'thick'  cloud.  Figures  3. 4.1. A 
and  C  illustrate  two  effects?  the  effect  of  cloud  thickness 
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The  global  radiation  and  dilfuse  radiation  decrease  with  the  change  in  the  cloud  thickness 
(100  m  thickness:  top  line;  500  m  thickness:  bottom  line).  The  wavelength  is  in  micro- 
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and  the  effect  of  changing  the  surface  albedo.  As  shown 
previously  with  the  broadband  cases,  the  effect  of  cloud 
thickness  is  to  continually  reduce  the  global  radiation  and 
to  first  increase  and  then  decrease  the  diffuse  radiation. 
These  same  trends  are  observed  in  the  spectral  plots.  The 
largest  changes  in  the  spectral  cases  occur  in  the  region  of 
the  atmospheric  window;  although  changes  occur  throughout 
the  spectrum. 

Comparing  the  global  and  diffuse  radiation  plots  for 
cloud  thicknesses  less  than  100  meters  (Figures  3. 4.1. A  and 
D) ,  one  can  determine  that  the  effect  of  high  surface  albedo 
is  to  produce  minimal  change  in  the  global  and  diffuse 
radiation  reaching  the  surface.  The  decline  in  the  global 
and  diffuse  radiation  with  thickness  of  the  cloud  occurs 
slowly  with  a  high  albedo  surface.  In  comparison,  a  low 
albedo  surface  shows  more  pronounced  variations  with 
thickness  changes  in  both  the  global  and  diffuse  radition. 
High  surface  albedos  tend  to  slow  the  reduction  in  both 
global  and  diffuse  radiation  due  to  thickness  changes. 

3.4,2  Spectral  Effects  Of  Cloud  Amount  And  Surface  Albedo 

Figures  3. 4. 2. A  and  B  demonstrate  the  spectral  effects 
of  varying  the  cloud  amount.  In  the  cases  presented,  the 
cloud  thickness  was  a  constant  100  meters,  the  cloud  base 
was  100  meters  square,  arJ  the  distance  between  the  clouds 
was  varied  to  obtain  the  desired  cloud  amount.  For  the  cases 
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bOX  stratus:  middle  line;  overcast:  bottom  line).  The  diffuse  radiation  becomes  larger 
with  each  increase  in  cloud  amount  (clear  sky:  bottom  line;  507„  stratus:  middle  line; 
overcast:  top  line).  The  wavelength  is  in  micrometers. 
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shown,  a  stratus  cloud  was  used.  The  degree  of  variation  is 
most  pronounced  for  a  low  albedo  (wet  sand)  surface.  Varying 
the  cloud  amount  over  a  snow  surface  produces  larger 
variations  in  the  diffuse  radiation  than  by  varying  the 
thickness.  The  reason  is  the  sky  is  not  completely  filled 
with  cloud.  The  diffuse  radiation  increases  in  proportion  to 
the  amount  of  cloud  cover.  By  the  same  token,  the  direct 
radiation  decreases  proportionately.  The  result  is  a  larger 
variation  in  the  diffuse  component  for  a  change  in  cloud 
amount  when  compared  to  a  change  in  cloud  thickness. 

The  surface  albedo  exhibits  similar  effects  with  the 
change  in  cloud  amount  as  it  did  with  the  change  in 
thickness.  The  global  radiation  tends  to  decline  more  slowly 
with  increasing  cloud  amount  over  a  high  albedo  surface  than 
over  a  low  albedo  surface.  The  diffuse  radiation  increases 
more  rapidly  over  a  high  albedo  surface  with  changes  in 
cloud  amount  than  when  over  a  low  albedo  surface.  The  low 
albedo  surface  tends  to  cause  significant  differences  in  the 
global  radiation  for  a  change  in  cloud  amount  while  causing 
relatively  small  changes  in  the  diffuse  component.  A  high 
albedo  surface  tends  to  cause  relatively  small  changes  in 
the  global  radiation  and  relatively  large  changes  in  the 
diffuse  radiation  for  a  change  in  cloud  amount. 


3.4.3  Spectral  Effects  Of  Cloud  Types 
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The  effects  of  cloud  type  and  surface  albedo  can  be 
seen  in  some  of  the  previous  figures.  For  comparison 
purposes,  the  global  radiation  for  a  stratus  overcast  and  a 
nimbostratus  overcast  for  a  wet  sand  surface  are  shown  in 
Figure  3. 4. 3. A.  A  wet  sand  surface  was  chosen  for  this 
section  because  the  lower  albedo  surface  causes  the  most 
differentiation  between  the  cases.  The  cloud  thicknesses  are 
25  meters  and  100  meters  respectively.  The  changes  due  to 
the  cloud  type  become  larger  as  the  cloud  thickness 
increases.  The  difference  in  the  global  radiation  at  the 
peak  in  each  curve  is  approximately  20  percent  of  the 
stratus  value  for  a  cloud  thickness  of  25  meters.  For  a 
cloud  thickness  of  100  meters,  the  difference  is 
approximately  25  percent.  Depending  on  the  spectral 
interval,  the  percentage  may  be  higher  or  lower.  However, 
the  most  significant  effects  take  place  in  the  region 
between  0.4  and  0.7  micrometers. 

Figure  3.4. 3. B  highlights  the  significance  of  cloud 
type  as  it  relates  to  cloud  amount.  The  plots  demonstrate 
about  a  15  percent  difference  in  the  global  radiation 
between  the  peaks  of  each  curve  for  50  percent  cloud 
coverage.  For  an  overcast,  the  difference  at  the  peak  when 
compared  to  the  value  of  the  stratus  curve  at  the  peak  is 
approximately  33  percent.  The  degree  of  influence  of  cloud 
type  is  obviously  interrelated  with  the  other  parameters. 
However,  a  reasonable  estimate  of  the  change  caused  by  cloud 
type  through  either  cloud  thickness  or  cloud  amount  may 
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range  from  15  to  50  percent  depending  on  the  precise 
circumstances. 

3 . 5  Comparison  With  An  Empirical  Model 

The  Air  Resources  Laboratory  (ARL)  developed  an 
empirical  model  for  hourly  solar  radiation  based  on  solar 
zenith  angle,  cloud  amount,  and  preciptation.  This  model 
does  not  include  many  of  the  parameters  in  the  solar 
radiation  model  developed  for  this  dissertation.  Factors  of 
importance  which  are  left  out  of  the  ARL  model  are  cloud 
thickness,  cloud  type,  water  vapor  amount,  and  turbidity. 

The  ARL  model  calculates  the  clear  sky  and  cloudy  sky  hourly 
irradiances  in  kilojoules  per  square  meter.  The  empirical 
equation  for  the  clear  sky  model  is: 

SRC  =  Aq  +  cos(ZA)  +  A2cos2(ZA)  +  A3cos3(ZA) 

The  empirical  equation  for  cloudy  skies  is: 

SR  =  SRC  (  Cq  +  C20PQ  +  C^OPQ2  +  C4OPQ3  *  C.RN) 

where 

SRC  -  Solar  radiation  for  clear  skies  (kilo-joules  per  square 
meter) 

SR  =  Solar  radiation  for  cloudy  skies  (kilo-joules  per  square 
meter) 

ZA  =  Solar  zenith  angle 

OPQ  =  Number  of  tenths  of  opaque  cloudiness  divided  by  10 

RN  =  Precipitation,  0  when  no  precipitation  is  reported  and  1 
when  some  form  of  precipitation  is  reported 
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The  equation  for  cloudy  skies  is  for  opaque  clouds  and  is 
the  form  of  the  equation  used  when  the  number  of  minutes  of 
sunshine  is  not  available.  When  the  number  of  minutes  of 
sunshine  is  available,  another  empircal  equation  is  used 
which  contains  the  constant  C,.  The  number  of  minutes  of 
sunshine  was  not  included  in  this  model,  so  the  above 
equation  was  used  for  comparison.  The  empirical  constants  in 
the  equation  were  determined  for  specific  sites  around  the 
United  States.  There  are  two  sets  of  coefficients;  one  for 
the  morning  and  one  for  the  afternoon.  Ann  Arbor,  Michigan 
was  not  one  of  the  sites  for  which  coefficients  were 
determined.  This  comparison  will  be  completed  for  Madison, 
Wisconsin.  The  empirical  constants  for  the  equations  can  be 
obtained  for  a  number  of  locations  for  each  month  of  the 
year  from  the  National  Climatic  Center,  Air  Resources 
Laboratory,  or  the  Solar  Energy  Measurements  and 
Instrumentation  Training  Course  from  the  Department  of 
Atmospheric  and  Oceanic  Science  at  The  University  of 
Michigan. 

There  are  a  number  of  parameters  which  must  be  applied 
to  the  solar  model  that  cannot  be  used  in  the  empirical 
model.  A  list  of  the  parameters  used  for  a  horizontal 
surface  follow: 


Latitude 

Longitude 

Day 

Solar  Zenith  Angle 
Surface  Type 
Cloud  width 
Cloud  separation 


43.083  degrees  N 
89.416  degrees  W 
166 

31.677  degrees 
Bare  soil 
500.0  meters 

Varies  depending  on  cloud  amount 


i 


CLOUD 

AMOUNT 

(TENTHS) 

REGRESSION 

MODEL 

(MJ/M2) 

SOLAR 

MODEL 

(MJ/M2) 

CLOUD 

THICKNESS 

(M) 

WATER 

VAPOR 

(CM) 

TURBIDITY 
( 500  nm) 

CLEAR 

3.038 

3.243 

0.5 

0.1 

3.196 

- 

1 .0 

0.1 

3.122 

- 

2.0 

0.1 

2 

3.010 

2.838 

50.0 

2.0 

0.1 

2.642 

150.0 

2.0 

0.1 

2.473 

2.0 

0.1 

5 

2.977 

2.413 

50.0 

2.0 

0.1 

1.924 

150.0 

2.0 

0.1 

1 .499 

300.0 

2.0 

0.1 

6 

2.967 

2.270 

50.0 

2.0 

0.1 

1.683 

150.0 

2.0 

0.1 

1.252 

300.0 

2.0 

0.1 

7  2.957 

2.127 

50.0 

2.0 

0.1 

1.442 

150.0 

2.0 

0.1 

1.084 

300.0 

2.0 

0.1 

10 

2.930 

1.742 

50.0 

2.0 

0.1 

0.947 

150.0 

2.0 

0.1 

0.536 

300.0 

2.0 

0.1 

Table  3. 5. A 


Regression  Model  Versus  Solar  Spectral  Model  Calculations 
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Table  3. 5. A  shows  the  comparison  of  the  regression  or 
empirical  model  versus  the  solar  spectral  model.  The 
comparison  demonstrates  differences  due  to  the  water  vapor 
amount  for  a  clear  sky  and  differences  due  to  cloud 
thickness  for  cloudy  skies.  The  results  in  Table  3. 5. A  are 
in  mega joules  per  square  meter.  The  empirical  model  gives 
reasonable  results  based  on  the  few  parameters  used  in 
calculating  the  solar  radiation.  The  solar  spectral  model 
demonstrates  similar  results  and  shows  the  variation  due  to 
cloud  thickness  and  cloud  amount  as  well.  While  an  empirical 
model  may  be  used  for  very  rough  estimates  of  the  amount  of 
solar  radiation,  they  do  not  represent  the  variation  in  the 
solar  radiation  which  would  be  of  interest  to  someone 
building  a  solar  collector  or  other  solar  oriented  device. 
Solar  models  which  have  more  physical  parameters  allow  one 
to  better  evaluate  the  changes  due  to  a  variety  of  factors. 

3.6  Radiance  At  The  Top  Of  The  Atmosphere 

This  model  calculates  the  solar  radiance  at  the  top  of 
the  atmosphere  as  might  be  received  by  a  satellite  sensor. 
There  are  no  known  sets  of  observed  spectral  data  of  the 
type  needed  to  compare  with  these  calculations.  However, 
there  are  calculations  of  the  radiance  at  the  top  of  the 
atmosphere.  Bowker,  Davis,  Myrick,  Stacy,  and  Jones  (1985) 
presented  a  plot  of  such  calculations  (See  Figure  3. 6. A).  It 
is  impossible  to  reproduce  the  plot  with  the  solar  spectral 
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model  due  to  many  unknown  factors  and  differences.  The  solar 
zenith  angle  directly  influences  the  amount  of  illumination 
of  a  target  at  the  surface.  While  the  solar  zenith  angle  is 
known,  the  starting  extraterrestrial  solar  irradiance  is  not 
and  there  are  several  possible  data  sets  in  use.  The  amount 
of  surface  illumination  is  critical  in  determining  the 
radiance  at  the  top  of  the  atmosphere.  Other  factors  of 
importance  are  the  aerosol  distribution,  turbidity,  water 
vapor  amount,  etc.  To  illustate  that  the  solar  spectral 
model  computes  radiances  in  the  correct  range  of  values,  two 
cases  were  calculated.  One  had  a  constant  spectral  albedo  of 
0.7  and  the  other  case  had  a  value  of  0.1.  The  location  was 
at  20  degrees  north  latitude  and  assumed  to  be  at  an  hour 
angle  of  0.0  degrees  (local  solar  noon).  The  two  cases  are 
plotted  in  Figure  3.6.B.  When  one  compares  the  solar 
spectral  model  results  with  Figure  3. 6. A,  one  can  see  that 
the  model  radiances  are  higher  in  one  portion  of  the 
spectrum  and  lower  in  another.  Figure  3. 6. A  had  the 
absorption  features  'removed  for  clarity’.  This  accounts  for 
some  of  the  differences  between  the  plots.  In  Figure  3. 6. A, 
beam  radiance  is  the  radiation  which  comes  directly  from  the 
target  and  total  radiance  is  all  the  radiation  which  reaches 
the  top  of  the  atmosphere.  In  Figure  3.6.B,  global  radiance 
is  equivalent  to  total  radiance  and  direct  radiance  is 
equivalent  to  beam  radiance.  The  diffuse  radiance  is  due  to 
radiation  scattered  by  the  atmosphere. 
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Lg  -  Beam  radiance 

L.j,  -  Total  radiance 

p  -  Surface  albedo 

-  Observation  zenith 
angle 


Radiance  At 


d--  Solar  zenith  angle 
a/o-  Single  scattering  albedo 
V  -  Visual  range 
TOA  -  Top  of  atmosphere 

Figure  3. 6. A 

The  Top  Of  The  Atmosphere 


The  beam  radiance  is  the  radiation  received  directly  from  the  target. 
The  total  radiance  is  the  radiation  received  directly  from  the  target 
plus  any  radiation  scattered  to  the  sensor  by  the  atmosphere. 

(Bowker,  Davis,  Myrick,  Stacy,  and  Jones,  1985) 


MILLIWATTS  /  SQ  CM  /  S*  /  MICROMCTCR 
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RADIANCE  -  TOP  Of  ATVOSFFCRE 
CONSTANT  REFLECTANCE  0.7 


RADIANCE  -  TOP  OF  ATMOSPHERE 


Legend 

A  GLOBAL 
X  DIRECT 
□  DIFFUSE. 


Figure  3.6.B 

Radiance  At  The  Top  Of  The  Atmosphere  For  Surfaces 
With  Constant  Spectral  Reflectance 


The  wavelength  is  in  micrometers. 
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In  general,  it  appears  the  solar  spectral  model  does  a 
reasonable  job  of  estimating  the  radiance  at  the  top  of  the 
atmosphere  for  a  clear  sky.  It  is  impossible  to  determine 
the  accuracy  of  the  model  in  computing  radiances  for  the  top 
of  the  atmosphere  until  observations  have  been  obtained. 
However,  it  appears  the  model  does  generate  answers  within 
an  order  of  magnitude  based  on  the  very  loose  comparison 
with  the  Bowker,  Davis,  Myrick,  Stacy,  and  Jones  plot. 

3.7  Proqation  Of  Error  In  The  Model 

The  amount  of  error  propagated  by  the  model  based  on 
the  uncertainties  of  the  independent  variables  is  a  legimate 
concern.  An  analysis  was  performed  to  quantify  the  degree  of 
error  one  might  expect  if  one  of  the  independent  variables 
is  not  known  accurately.  This  analysis  will  demonstrate 
which  variables  need  to  be  measured  most  accurately  to 
minimize  model  error.  Because  portions  of  this  model  are 
based  on  numerical  techniques,  it  is  not  practical  to  obtain 
an  analytic  expression  for  the  errors  due  to  each 
independent  variable.  However,  an  estimate  of  error  can  be 
obtained  by  choosing  a  particular  set  of  parameters  and 
calculating  the  value  of  global  radiation  at  the  surface. 

The  global  radiation  values  compared  with  the  original  case 
will  determine  the  error  caused  by  uncertainties  in  the 
independent  variables. 
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The  particular  case  used  as  the  basis  for  this  analysis 
used  the  following  parameter  values: 


Latitude 

42.2833 

deg 

Longitude 

-83.7333 

deg 

Surface  Type 

Bare  Soil 

Hour  Angle 

0.0 

deg  (noon) 

Ozone 

0.3639 

cm 

Precipitable  Water 

2.0 

cm 

Turbidity 

0.15 

Cloud  Thickness 

300.00 

m 

CLoud  Width 

500.00 

m 

Cloud  Amount 

0.30 

percent 

Cloud  Height 

1000.00 

m 

Table  3. 7. A  lists  the  independent  quantities  which  were 
varied.  With  the  exception  of  the  surface  pressure,  all  the 
variables  were  deviated  by  10  percent  from  their  original 
value.  The  surface  pressure  was  assumed  to  be  known  within 
15  millibars.  Table  3. 7. A  gives  the  original,  high,  and  low 
values  for  each  variable. 

Table  3.7.B  shows  the  results  of  this  analysis  in  terms 
of  global  radiation.  The  percentage  error  is  the  difference 
between  the  new  global  value  minus  the  original  global  value 
divided  by  the  original  global  value.  By  scanning  the 
percentage  error  column,  it  is  quite  easy  to  determine  which 
variables  it  is  important  to  know  accurately.  Cloud  amount 
and  cloud  thickness  should  be  measured  as  precisely  as 
possible  if  one  expects  to  obtain  accurate  answers  from  the 
model.  The  variables  which  most  influence  the  model 
calculations  by  uncertainty  in  observation  are  (most 
influence  to  least  influence):  cloud  amount,  cloud 
thickness,  cloud  width,  precipitable  water,  turbidity, 
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VARIABLE 

ORIGINAL 

CASE  VALUES 

LOW 

VALUE 

HIGH 

VALUE 

PERCENTAGE 

OF  ORIGINAL  VALUE 

Ozone  (cm) 

0.3639 

0.3275 

0.4003 

10 

Free  ip  itab le 
Water  (cm) 

2.0000 

1.8000 

2.2000 

+ 

10 

Turbidity 

0.1500 

0.1350 

0.1650 

_ , 

10 

Surface 

Pressure  (mb) 

1013.0 

998.0 

1028.0 

+ 

1.5 

Cloud 

Thickness  (m) 

300.0 

270.0 

330.0 

10 

Cloud 

Width  (m) 

500.0 

450.0 

550.0 

- 

10 

Cloud 

Amount 

0.50 

0.40 

0.60 

- 

10 

Cloud 

Height  (m) 

1000.0 

900.0 

1100.0 

- 

10 

Figure  3. 7. A 


Assumed  Uncertainties  In  Model  Variables 
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VARIABLE 

ORIGINAL 

GLOBAL 

VALUE 

LOW 

GLOBAL 

VALUE 

PERCENTAGE 
ERROR  OF 
ORIGINAL 

HIGH 

GLOBAL 

VALUE 

PERCENTAGE 
ERROR  OF 
ORIGINAL 

Ozone 

316.799 

316.219 

0.1831 

317.393 

0.1875 

Precipitabl 

Water 

e  316.799 

315.536 

0.3986 

318.117 

0.4160 

Turbidity 

316.799 

315.555 

0.3927 

318.057 

0.3971 

Surface 

Pressure 

316.799 

316.455 

0.1086 

317.146 

0.1095 

Cloud 

Thickness 

316.799 

309.944 

2.1638 

325.507 

2.7487 

Cloud 

Width 

316.799 

315.936 

0.2724 

318.522 

0.5438 

Cloud 

Amount 

316.799 

273.892 

13.5439 

368.295 

16.2551 

C  loud 

Height 

316.799 

316.792 

0.0022 

316.817 

0.0057 

Table  3.7.B 

Propagation  Of  Error  In  The  Model 
The  global  radiation  values  are  watts  per  square  meter. 
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surface  pressure,  ozone  amount,  and  cloud  height.  The  cloud 
amount  and  cloud  thickness  are  by  far  the  two  most  important 
parameters  to  know  accurately. 


I 

f 
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CHAPTER  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 

A  spectral  insolation  model  was  designed  which  includes 
clear  skies,  cloudy  skies,  and  inclined  surfaces.  Several 
typical  spectral  surface  albedos  are  available  in  the  model 
as  well  as  several  different  cloud  types.  Only  one  aerosol 
distribution  (rural)  is  included  in  the  model.  Interesting 
features  of  the  model  include  the  ability  to  sum  the 
radiation  over  any  desired  wavelength  band  between  0.28  and 
4.0  micrometers,  and  the  ability  to  average  the  radiation 
over  time.  The  flexibility  designed  into  this  interactive 
model  is  one  of  its  strongest  attributes. 

The  solar  spectral  model  was  compared  against  actual 
observations  for  clear  skies.  The  observations  were  both 
spectral  and  broadband.  The  modelled  values  were  within  the 
accuracy  goals  of  ±  5%  for  the  global  radiation  and  +  15% 
for  the  diffuse  radiation.  The  aerosol  distribution  was  the 
most  important  factor  in  the  clear  sky.  If  the  aerosol 
distribution  of  the  atmosphere  is  not  compatible  with  the 
rural  aerosol  distribution,  then  the  model  error  will 
increase  substantially  depending  on  the  amount  of  disparity 
between  the  observed  and  modelled  aerosol  distributions.  A 
recommendation  for  further  study  is  to  include  more  than  one 
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aerosol  distribution  in  the  model.  A  second  recommendation 
is  to  use  more  than  one  turbidity  value  to  determine  the 
aerosol  distribution.  Turbidity  measurements  at  several 
different  wavelengths  should  be  used  to  more  accurately 
determine  the  observed  aerosol  distribution  and  to  determine 
which  model  aerosol  distribution  to  use. 

The  solar  spectral  model  was  compared  against  cloudy 
sky  observations.  The  disagreement  between  modelled  and 
observed  values  was  large.  Most  of  the  differences  are 
thought  to  be  due  to  the  inability  of  conventional  surface 
observations  to  adequately  describe  the  cloud  situation. 
Specialized  observations  are  required  in  order  to  compare 
the  calculations  from  a  model  of  this  detail  with 
observations.  The  observations  should  include  the  dimensions 
of  the  clouds,  the  geometry  of  the  clouds  relative  to  the 
instruments  taking  the  measurements,  and  the  time  the 
instruments  are  in  the  cloud  shadow  versus  the  direct 
sunlight.  The  model  should  be  modified  to  compute  a  time  as 
well  as  spatial  average  and  should  include  the  use  of  the 
number  of  minutes  of  direct  sunshine. 

The  cloudy  sky  comparison  and  the  cloud  albedo 
comparison  indicated  the  solution  for  the  cloud  radiative 
transfer  is  not  adequate.  This  was  exhibited  by  the  sharply 
reduced  radiation  transmitted  through  the  cloud  with 
increasing  cloud  thickness  and  the  distinctly  different 
shape  of  the  cloud  albedo  curve  when  compared  with 
Chel'tsov's  measurements.  While  Paltridge  and  Platt 
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recommended  the  solution  because  of  its  simplicity,  it  has 
shortcomings.  This  solution  is  primarily  useful  when  the 
absorption  is  small  and  when  the  optical  depth  of  the  layer 
is  large  (greater  than  5).  While  the  absorption  is  small  in 
the  visible  portion  of  the  spectrum  and  most  water  droplet 
clouds  in  the  low  atmosphere  have  large  optical  depths,  the 
simple,  two  stream  solution  appears  inadequate.  The  solution 
results  in  too  rapid  attenuation  through  the  cloud  and  too 
rapid  a  change  in  cloud  albedo  with  a  change  in  cloud 
thickness.  For  future  work,  I  would  recommend  another 
solution  to  the  radiative  transfer  equation  be  used.  It  was 
not  possible  to  establish  the  degree  of  inaccuracy  due  to 
the  lack  of  suitable  observations. 

Even  though  the  cloud  modelling  is  in  error,  the 
effects  of  varying  the  cloud  parameters  and  surface  types 
could  still  be  evaluated.  The  solar  zenith  angle  has  the 
greatest  effect  of  all  the  parameters.  The  solar  zenith 
angle  affects  the  surface  albedo,  the  path  distance  through 
the  cloud,  and  the  amount  of  direct  radiation  received  on  an 
inclined  plane.  The  effect  of  solar  zenith  angle  is 
demonstrated  in  the  modelled  albedos  of  overcast  clouds  at 
Archangel,  Soviet  Union  where  the  day  of  the  year  determined 
the  solar  zenith  angle  which  was  largely  responsible  for  the 
changes  in  the  calculated  albedos.  The  effect  of  solar 
zenith  angle  is  particularly  strong  as  the  zenith  angle 
becomes  large,  especially  in  cloudy  skies. 
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While  it  is  possible  to  choose  a  situation  in  which  any 
particular  parameter  can  be  demonstrated  to  have  a  large 
effect,  cloud  amount  has  the  largest  influence  of  any  of  the 
cloud  or  surface  parameters.  The  other  factors  in  decreasing 
order  of  influence  are  cloud  thickness,  surface  albedo,  and 
cloud  type.  In  typical  situations,  the  cloud  amount  will 
dominate  the  scenario  in  terms  of  the  radiation  reaching  the 
surface  of  the  Earth.  There  are  a  number  of  situations  where 
the  cloud  thickness  can  be  the  dominant  parameter.  These 
situations  are  usually  when  the  cloud  amount  is  large  and 
most  of  the  radiation  is  reaching  the  surface  through  the 
clouds.  Cloud  type  is  the  least  influential  of  the  cloud 
parameters  and  was  used  to  determine  the  cloud  droplet 
distribution.  As  long  as  the  cloud  drop  size  distribution  is 
broad  the  effects  of  clouds  are  generally  characterized 
fairly  well.  It  is  the  tail  of  the  drop  size  distribution 
for  the  large  drop  sizes  which  is  important.  Relatively 
small  changes  in  the  number  of  large  drops  produce  the 
greatest  changes.  As  a  consequence,  future  cloud  modelling 
should  concentrate  on  understanding  how  the  distribution  of 
the  largest  drops  change  with  cloud  type  and  cloud 
lifecycle.  These  same  thoughts  concerning  drop  size  are  also 
stated  by  Welch,  Cox,  and  Davis  (1980;  and  Ulaby,  Fung  and 
Moore  (1981). 

The  effects  of  surface  albedo  have  been  particularly 
interesting.  The  surface  albedo  is  influenced  by  solar 
zenith  angle.  For  large  solar  zenith  angles,  changing  the 
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surface  type  results  in  small  changes  of  the  solar  radiation 
reaching  the  surface.  High  surface  albedos  insure  the 
decrease  in  radiation  which  reaches  the  surface  is 
minimized.  Low  surface  albedos  exhibit  larger  changes  in  the 
radiation  reaching  the  surface. 

While  this  model  has  shortcomings,  it  is  a  high 
resolution  spectral  model  which  includes  clouds.  Although 
this  model  is  deficient  in  several  areas  due  to  the 
solutions  used  in  the  modelling,  it  provides  a  means  of 
assessing  the  influence  of  a  variety  of  parameters.  When 
compared  to  empirical  models  the  detail  and  variation  one 
can  obtain  in  the  answers  provides  more  information  for 
users.  Since  this  model  is  based  on  physical  principles,  it 
provides  insight  into  problems  where  empirical  models  fail 
short . 

In  summary,  this  model  provides  a  high  degree  of 
spectral  resolution  for  both  clear  and  cloudy  skies,  while 
the  two  stream  solution  for  clouds  is  inadequate,  the 
general  results  of  the  calculations  can  be  used  to  infer 
amounts  of  radiation  as  well  as  general  trends.  Analysis  via 
changing  parameters  can  be  of  particular  help  to  designers 
of  solar  and  remote  sensing  equipment.  Future  research 
should  attempt  to  better  model  the  cloud  radiative  transfer 
properties.  This  model  did  not  include  the  radiation 
reflected  from  the  sidewalls  of  the  clouds.  The  reflection 
from  the  cloud  sidewalls  should  be  included  in  future  work. 
The  technique  of  using  simple  geometry  clouds  and  weighting 


the  respective  areas  of  the  clouds  by  the  radiation  entering 
those  areas  seems  to  have  merit.  A  more  definitive 
assessment  of  the  technique  cannot  be  made  until  more 
detailed  observations  are  taken. 
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QUANTITIES: 

A 

E 

H 

o 

I 

L 
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*1 


e 


2 


9 


3 


9 


4 


P 

1 

<t> 


CALCULATION  OF  THE  IRRADIANCE  ARRIVING  AT 
THE  SATELLITE  OBSERVATION  POINT 


-  Area  of  (object)  (m  ) 

2 

-  Irradiance  on  (object)  (W/m  ) 

2 

-  Irradiance  on  top  of  atmosphere  (W/m  ) 

-  Intensity  (W/Sr) 

2 

-  Radiance  of  (object)  (W/m  -Sr) 

2 

-  Reflected  irradiance  from  (object)  (W/m  ) 

-  Distance  from  satellite  to  target  (m) 

-  Solar  zenith  angle  relative  to  a 
flat  surface 


-  Solar  zenith  angle  relative  to  the 
normal  of  a  tilted  plane 

-  Satellite  zenith  angle  relative  to  the 
normal  of  a  tilted  plane 

-  Satellite  zenith  angle  relative  to  a 
flat  surface 

-  Albedo  or  reflectance  of  (object) 

-  Pi  (3.14159) 

-Flux  (Watts) 


Tatml 


Total  atmospheric  transmittance  along  the 
path  from  the  sun  to  the  target  location 


Tatm2 


total  atmospheric  transmittance  along  the 
path  from  the  target  location  to  the 
sate l 1 ite 


13  -  Solid  angle  of  (object) 
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DERIVATION: 

(1)  Calculate  flux  caught  by  plane: 

^(caught  by  plane)  =  Hq  A(plane)  cosCfl^)  Tatm2 

(2)  Irradiance  on  plane: 

E(on  plane)  =  d>  (caught  by  plane)  /  A(plane) 

(3)  Irradiance  reflected  from  plane: 

M(from  plane)  =  p(plane)  E(on  plane) 
ASSUMING  Lambertian  type  scatter: 

L  (plane)  =  M(from  plane)  /  n 

(4)  Intensity  in  the  direction  of  view: 

Kin  the  direction  of  view)  =  L(plane)  A(plane)  cos(d^) 


(5)  Flux  at  the  satellite  observation  point: 


( sate  1 1 ite )  -  I(in  the  diri  :tion  of  view)  Tatm2  .  A  (sensor)  /  R 
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(6)  Irradiance  at  the  sateLLite  observation  point: 

E(satellite)  =  <t>  (satellite)  /  A(sensor) 


(7)  Solve  for  E(satellite)  by  backward  substitution  from  step  (6): 


,  1  A(sensor)  A( plane )  .cos  (0,). 

E( satellite)  =  A(sensor)  2  A(plane)  ATM 2 

Rs 

P (plane) .H^ .A(plane).cos(02^ '  ath1 

71 


(8)  Reduce  Equation: 


E( sate l 1 ite )  =  A( Plane) ■  .  cos(  0., ) .  cosi  0., ) .  /0(plane).H 

2  TATM1  tATM2  J  t  ° 

K 

S 


2 

Note:  Solid  angle  of  plane  =Q(plane)  =  A(plane)  /Rg 

L(of  plane)  =p(plane)  /t 


E(satellite)  =  fj(plane).  Tatmi’  TaCm2  C05^  &  2^  cos  ^  ^3^  L(plane) 
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(9)  Determining  the  irradiance  from  the  atmosphere  : 

The  irradiance  arriving  at  the  satellite  consists  of  two 
components 

1.  Radiation  from  the  target  plane 

2.  Radiation  from  the  atmosphere 

The  viewing  geometry  will  determine  the  radiance  of  the 
atmosphere.  To  simplify  calculations,  two  assumptions  have 
been  made. 

Assumption  1:  If  the  angle  formed  by  connecting  the  sun, 

target  plane  location,  and  the  satellite  is  less 
than  or  equal  to  90  degrees,  then  the 
atmospheric  radiance  is  calculated  from  the 
radiation  backscattered  by  the  atmosphere. 

Otherwise,  if  the  angle  formed  is  greater  than 
90  degrees,  then  the  atmospheric  radiance  is 
derived  from  the  radiation  forward  scattered  by 
the  atmosphere. 

Assumption  2:  Once  the  amount  of  atmospheric  radiation 

backscattered  (forward  scattered)  is  computed, 
the  atmospheric  radiance  is  calculated  by 
assuming  the  radiation  uniformly  scattered  over 
the  backscattering  (forward  scattered) 
hemisphere  (»r  steradians). 

Summary:  If  the  ( Sun-targe t-sate 1 1 i te )  angle  is  less  than 

or  equal  to  90  degrees  then  the  atmospheric 
radiance  is: 

L( atmosphere)  =  E( backscattered)  /  t 
Otherwise : 

L( atmosphere )  -  E(forward  scattered)  /t 

(10)  Total  irradiance  arriving  at  the  satellite  observation  point: 

E(  sate  1 1  i  te )  =fl(Plane).  TisiTM2'C0S  'C0S^2^ '  (PIane)*H0 


where 

2 

fl(sensor)  =  A(sensor)  /R 

s 


L(atmosphere) fl( sensor) 
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Figure  A.l 

Solar  Zenith  Angles  Relative  To  Flat  And  Tilted  Surfaces 

H  =  Extraterrestrial  solar  irradiance 
°A 

0 ^  =  Solar  zenith  angle  relative  to  flat  surface 

6 2  =  Solar  zenith  angle  relative  to  tilted  surface 

=  Satellite  zenith  angle  relative  to  tilted  surface 

0 ^  =  Satellite  zenith  angle  relative  to  flat  surface 

tATM^  =  Transmission  for  the  atmospheric  path  between  the  sun  and 
target  plane 

TATMj  =  Transmission  for  the  atmospheric  path  between  the  satellite 
and  the  target  plane 
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Figure  A. 2 


Sources  of  radiance  reaching  the  satellite 
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L(Atmosphere)  =  L( Backsca t te red )  =  E( Backscattered)  /t 

Figure  A. 3 

When  the  radiance  of  the  atmosphere  is  derived  from 
backscattered  radiation 
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L(Acmosphere)  =  L( Forward  Scattered)  =  E( Forward  Scattered)  !i 


Figure  A.^ 

When  the  radiance  of  the  atmosphere  is  derived  from  forward 
scattered  radiation 
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APPENDIX  B 

THE  ADDING  TECHNIQUE  FOR 
COMBINING  LAYERS 


TERMS  AND  SYMBOLS: 


TERMS: 

D  =  Downward  transmission  incident  on  top  of  bottom  layer 

R^  =  Reflectivity  of  top  layer 

R^  =  Reflectivity  of  bottom  layer 

Rj2  =  Combined  reflectivity  of  the  two  layers 

T^  =  Transmissivity  of  top  layer 

T2  =  Transmissivity  of  bottom  layer 

T 1 2  -  Combined  transmission  of  the  two  layers 

U  =  Upward  transmission  incident  on  bottom  of  top  layer 
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ADDING  TECHNIQUE  WITH  SOURCE  ABOVE  THE  LAYERS 


Diagram  of  the  components  of  radiation  using  the  Adding 
technique  for  combining  two  layers  with  known  transmissivities 
and  reflectivities  with  the  source  above  the  layers. 


EQUATIONS: 

R 


12 


R1  *  T1  R2  *  T1  R2RiR2  +  T1  R2RiR2RiR2  + 
=  Rj  *  Tt2R2  (1  *  R^  +  (R1R2)2*  ...) 


=  R. 


Tt  R  /(l-R,R  ) 


U  =  V2  .  TlR2R2R2  *  t1r2r1r2r1r2 

=  TrR2  (1  *  R2Rl  *  (R2Rl)2  *  ... 
=  T1R2/(1-R1R2) 


d  -  Tt  .  TlR2Rl  *  t1h2R1R2R1  * 

«  Tj  (1+R2R1  *  (R2R1)2  -  ...) 


t1/(i-r1r2) 


T12  '  T1T2  +  T1T2R2R1  +  TlT2R2RlR2Rl  * 
*  TtT2  (1  *  R2Rt  *  (R2R1)2  +  ...) 


TtT2/( l-RlR2) 
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ADDING  TECHNIQUE  WITH  THE  SOURCE  BELOW  THE  LAYERS 


Figure  B.3 

Diagram  of  the  components  of  radiation  using  the  Adding 
technique  for  combining  two  layers  with  known  transmissivities 
and  reflectivities  with  the  source  below  the  layers. 


=  a2  +  T2  R1  +  T2  R1R2R1  +  T2  RiR2R2  + 

=  r2  +  t22rl  (i  *  r1r2  +  (r1r2)2  +  ...) 

=  r2  ♦  t22r1/(i-r1r2) 


u 


=  T2  "  T2R1R2  +  T2R1R2R1R2 
=  t2  (i  ♦  r1r2  +  (r  r  )2  * 


T2/(1'R1R2) 
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ADDING  TECHNIQUE  WITH  THE  SOURCE  BETWEEN  THE  LAYERS 


Diagram  of  the  components  of  radiation  using  the  Adding 
technique  for  combining  two  layers  with  known  transmissions  and 
reflectivities  with  source  between  the  layers. 


EQUATIONS: 

R 


12 


-  T2Rt  4.  T^R^  +  T2RiR2R1R2R1 


=  t2r1  (i  +  r2r1  4.  (r2rl) 
=  t2r1/(i  -  rlr2) 


U  -  Rj R2  +  R!R2R1R2  +  R1R2R1R2R1R2  + 
=  R1R2  (1  +  R:R2  4.  (RtR2)2  4-  .  .  .  ) 

=  R  R2/(l  -  R1R2) 


D  -  Rt  +  RiR2R1  +  R1R2R1R2R1  + 

=  Rj  (1  4-  R  R2  *  (R1R2)2  4-  .  .  .  ) 

=  R^d  -  R  R2) 


12 


=  T!  4-  TiRiR2  +  t1r1r2RiR2  +  t1r1r2r1R2RiR2  * 


=  Tj  (1  4  R1R2  4.  (RtR2)  +  ...) 


=  Tj/d  -  R  R2) 
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APPENDIX  C 

GEOMETRIES  ASSOCIATED  WITH  ILLUMINATING  AN  ARBITRARILY 
INCLINED  PLANE  ON  THE  EARTH 

Adapted  and  condensed  directly  from:  Solar  Energy  Engineering 

Edited  by  A.A.M.  Sayigh 
Chapter  2  by  Enrico  Coffari 
Academic  Press,  New  York,  1977 


TERMS  AND  SYMBOLS: 

A  -  Azimuth  measured  from  true  (geographical)  north, 
positive  in  the  eastward  (clockwise)  direction. 

Range:  0  to  360  degrees 

0  to  2  t  radians 

h  -  Hour  angle 

(  +  )  westward  from  the  noon  meridian 
(-)  eastward  from  the  noon  meridian 
Range:  +180  to  -180  degrees 

+  n  to  n  radians 

h  -  Hour  angle  of  sunrise 

sr  ° 

h  -  Hour  angle  of  sunset 

ss 

i  -  Elevation  angle  of  the  incident  radiation  on  the 

inclined  plane  measured  from  the  surface  of  the  plane  to 
the  position  of  the  sun. 

(+)  angle  in  front  of  plane's  surface  (sun 
i 1 luminated) 

(-)  angle  behind  the  plane's  surface 
Range:  -90  to  -90  degrees 

+  T/2  to  -T/2  radians 

a  -  Azimuth  angle  the  inclined  plane  is  facing  in  the 
coordinate  system  used  by  E.  Coffari 

(+)  eastward  azimuth  from  true  north 
(-)  westward  azimuth  from  true  north 
Range:  +90  to  -90  degrees 

+  T/2  to  -nH  radians 

Note  the  range  does  not  need  to  be  larger  due  to  the 
definition  of  the  zenith  angle. 

/I  -  Solar  elevation  angle 

Range:  +90  to  -90  degrees 

-T/2  to  T/2  radians 
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a  -  Zenith  angle  of  the  inclined  plane  measured  from  the 
local  vertical  to  the  normal  vector  of  the  plane's 
surface . 

(  +  )  southward  facing  plane 
(-)  northward  facing  plane 
Range:  +90  to  -90  degrees 

+  */2  to  ~nt 2  radians 

S  -  Declination  of  the  sun 
(+)  north 
(-)  south 

Range:  +90  to  -90  degrees 

+  T/2  to  -T/2  radians 

©  -  Latitude 

(+)  north 
(-)  south 

Range:  +90  to  -90  degrees 

+  T/2  to  -  71  / 2  radians 

SOLAR  ALTITUDE  (ELEVATION): 

sin/}=  sin5sin©  +  cos  8  cos  h  cos©  (C.l) 

At  solar  noon  the  hour  angle  is  zero  and  Eqn  (C.l)  reduces  to: 

sin/J=  sin«5sin?>  +  cos  S  cos  ©  =  cos  (5-©)  =  cos  ( q>-8 )  (C.2) 

(J  =  arcsin  (  cos(<5-©)  )  =  arcsin  (  cos(®-<5)  )  (C.3) 


Beta  (  fi  )  ranges  from  0  to  90  degrees  with  negative  values  of 
meaning  that  Che  sun  is  below  the  horizon  (polar  night). 

SOLAR  AZIMUTH 


The  solar  azimuth  angle  A  measured  from  true  (geographical)  north 
positively  in  the  eastward  direction: 

A  =  arctan  ( _ s^-n  h _ )  (C.9) 

-  cos  ©  tan  A  +  sinocos  h 

with  -90  -  A  il  +90  degrees 

To  obtain  the  correct  azimuth: 

If  h>0  (afternoon  azimuth)  and  A<180  then  A  =  A  +  180 

If  h<0  (morning  azimuth)  and  A>160  then  A  =  A  -  180 

SUNRISE  AND  SUNSET  TIMES 

h  =  arccos  (  -  tan  i  tan  »)  (C.5) 

h  =  -h 
sr 


ss 


=  h 
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But  if: 


-  tan  &  can  <3  +1 


=  +1 


=  -1 


<r 


-1 


SUNRISE  AND  SUNSET  AZIMUTHS 


The  sun  will  neither  rise  nor  set  for 
the  day,  polar  night. 

The  sun  will  be  on  the  horizon  for  an 
instant  only  (at  solar  noon).  This 
occurs  the  last  day  before  a  polar 
night  and  the  first  day  after. 

The  sun  will  be  on  the  horizon  for  an 
instant  only  (at  midnight,  h=  _^180). 
This  occurs  che  day  before  and  the  day 
after  a  polar  day. 

The  sun  will  neither  rise  nor  sec  for 
the  day;  polar  day. 


A  =  arcsin  (  -  sin  h  cos  8  ) 


( C .  6 ) 


If  8  >  0  then  the  sun  rises  and  sets  in  che  northern  quadrants. 

If  S  =  0  then  the  sun  rises  due  east  and  sets  due  west. 

If  S  <.  0  then  the  sun  rises  and  sets  in  the  southern  quadrants. 


CALCULATION  OF  THE  ANGLE  OF  INCIDENCE  OF  DIRECT  RADIATION  ON  AN  INCLINED 
PLANE 


cos  i 


cos  i 


where 


(cos  ?  sin  ®  -  s  in  r  cos  cos  ) 
+•  (cos  :  cos  ®  +  s in  s  s in  0  cos  a  ) 
+  sins  sinacos  5  sin  h 


s  i  n  £ 

cos  S  cos  h 


(C.7) 


=  sinS  +  cos  5  cos  h 


C cos  5  s  in  h 


Cj  =  cos-  sin<3  -  s  in  s;  cos  <3  cos  ct 
C2  =  cos-gos<3  +  sin  ;  sin  iocosa 


C^  =  s in  •  s in  a 

ZERO  IRRADIATION  SOLAR  TIMES,  ALTITUDES,  AND  AZIMUTHS  OF  AN  INCLINED 
PLANE 


From  Eqn  (C.7),  E.  Coffari  developed  a  solution  by  transforming  the 
sines  and  cosines  into  tangents,  and  because  of  the  one-to-one 
correspondence  of  tangent  functions  wihtin  their  period,  the 
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mathematical  solutions  coincide  with  the  physical  positions  of  the  hour 
angle  of  the  sun  at  zero-irradiation  times. 

The  two  zero-irradiation  hour  angles  of  a  given  day  and  for  a  given 
tilted  plane  divide  the  circle  representing  the  24  hours  of  the  dav  into 
two  sectors:  in  one  the  sun  is  above  the  plane,  in  the  other  it  lies 
under  the  plane.  The  lack  of  real  solutions  indicate  that  the  sun  is 
above  or  below  the  tilted  plane  for  the  whole  day.  To  determine  which 
sector  represents  the  time  of  day  when  the  sun  is  above  (or  under!  the 
tilted  plane,  or,  in  the  latter  case,  to  determine  whether  the  plane  is 
in  the  sun  or  in  the  shade  for  the  day,  it  is  sufficient  to  calculate 
the  angle  of  incidence  of  the  solar  radiation  for  any  hour  angle  except 
for  the  hour  angles  of  the  zero-irradiation  times.  This  is  most 
conveniently  done  for  solar  noon  as  Eqn  (C.7)  reduces  to: 

cos  i  =  (cos  t  sir©  -  sin  i  cos  0  cos  a  )  sin  <  C . 3  ’ 

♦  (cos  tcostf  *  s  in  t:  s  in  id  cos  a  )  cos  £ 

and  if  cos  i  >0,  the  sector  in  which  the  selected  hour  angle  lies  will 

be  the  sector  in  which  the  sun  is  above  the  tilced  plane. 

If  cos  i<  0,  the  will  be  under  the  tiltec  plane  in  the  other  sec'. 

Sunrise  and  sunset  also  divide  the  circle  representing  the  I-  hours 
of  a  day  into  two  sectors.  The  sector!  3  )  in  which  the  sur.  is  above 
both  the  tilted  plane  and  the  horizon  is  (are'  the  one1  s1  in  which  the 
surface  is  exposed  to  the  sun. 

Eqn  (C.7)  may  be  reduced  to  a  form  simpler  to  handle  bv  letting: 
x  =  sin  ;  sin  it  cos  i 

v  =  (cos  ;cosu  -  sin  b  sin  ocos  'i  '  co?  6 

2  -  icos  :  sin  u  -  sin  :  cos  u  cos  a  1  sini 

cos  i  •  x  sin  h  •  v  cosh  -  2  -  d  I .  - 
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The  condition  of  parallelism  is  reached  when  the  solar  beam  lies  on  the 
plane  of  the  slope,  i.e.  the  angle  of  incidence  is  equal  to  90  degrees, 
or  cos  i  =  0: 

x  sin  h  +  y  cos  h  +  z  =  0  (C.10) 
and  the  hour  angles  h  of  the  zero-irradiation  times  may  be  calculated  as 
follows,  letting  h  =  2ur- 


sin  h  =  sin  2af~  2  sinorcos ur=  2 


tan  or 


Vi  2 

V  1  +  tan  a r 


'v7 


tan  ur 


V 


tan  or 


1  -f  tan 


2  2 

cos  h  =  cos  2 ur=  cos  ur-  sinur= 


2 

tan  co 


tan  ur 


„  2 

tan  ur 


1 


tan 


ur 


1  +  tan  ur 
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Substituting  the  above  equations  in  Eqn  (C.10),  we  obtain 


and  rearranging,  we  have 


(z  -y)  tan +■  2x  tan ur  +  (y  +  z)  =  0 


hj  -  2  arctan(  tanor^)  =  2  arctan(  tanar^)  (C.ll) 
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The  hour  angles  h^  and  h^ ,  are  the  angles  at  which  the  sun's  rays  are 
parallel  to  the  inclined  plane;  complex  and  conjugate  solutions  indicate 
that  the  sun's  rays  are  never  parallel  to  the  given  plane  for  the  day. 

The  next  step  involves  calculating  the  cosine  of  the  angle  of 
incidence  for  an  hour  angle  that  is  not  a  zero-irradiation  hour  angle, 
by  using  Eqn  (C.9).  For  simplicity,  we  shall  calculate  the  angle  of 
incidence  of  the  solar  radiation  for  an  hour  angle  of  0  degrees, 

i 

corresponding  to  solar  noon. 

If  cos  i >0,  the  plane  will  be  exposed  to  solar  radiation  at  solar 
noon  and  hence,  the  sun  will  shine  on  it  from  h^  or  from  sunrise, 
whichever  occurs  later,  and  until  h^  or  sunset  whichever  occurs  earlier 
(from  sunrise  to  sunset  if  h^  and  1^  are  complex  numbers). 

If  cos  i-c  0,  the  tilted  plane  is  in  the  shade  at  solar  noon  (all  day 
if  h^  and  h^  are  complex  numbers),  and  will  be  exposed  to  the  sun  from 
to  h^  (clockwise  direction)  except,  of  course,  for  the  interval  of 
time  between  sunrise  and  sunset. 

The  solar  altitude  and  azimuth  at  the  hour  angle  corresponding  to 
the  zero-irradiation  times  may  be  calculated  from  Eqns  (C.l)  and  (C.4). 
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APPENDIX  D 

TWO  STREAM  SOLUTION  OF  THE  EQUATION  OF  RADIATIVE  TRANSFER 
(Adapted  From  Radiation  Commission,  1975) 

TERMS  AND  SYMBOLS: 

A  =  Albedo 

b  =  Backscattered  radiation 

E+  =  Downward  radiation  flux 

E  =  Upward  radiation  flux 

f  =  Forward  scattered  radiation 

p(f*)  =  Scattering  phase  function 

T  =  Transmission 

t  =  Optical  thickness  of  parallel  slab 
uT0  =  Single  scattering  albedo 
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TWO-STREAM  APPROXIMATION  (OR  SHUSTER-SCHWARZSCHILD  METHOD) 

The  two-scream  approximation  has  been  used  by  many  authors  to 
achieve  a  quick  approximate  solution  to  the  equation  of  transfer  by 
decomposing  the  radiation  field  into  two  opposing  streams.  The 
treatments  vary  slightly,  and  a  comparison  of  the  various  methods 
follow. 

Chu  and  Churchill  (1955)  consider  a  beam  passing  through  a  parallel 
slab  of  optical  thickness  ,  incident  on  the  surface  at  an  angle  0Q. 
The  radiation  is  broken  into  two  fluxes:  the  first  is  in  the  direction 
of  the  incident  beam,  while  the  second  is  in  the  backwards  direction. 

For  these  fluxes  the  equation  of  transfer  integrated  over  azimuth 
yields  the  two  equations 


u  =|cos  e  | 
o  o 

Here,  is  assumed  positive  for  the  downwards  directions  and  the  factors 
b  and  f  represent  the  fraction  of  radiation  scattered  backward  and 
forward,  respectively.  Clearly,  we  have 
b  +  f  =  1 

and  Chu  and  Churchill  make  the  logical  choice 


f 


(u) 


c  u 


i 
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A 


=  E  (o)  = 
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The  transmission  (including  the  incident  beam)  is  given  by 


T  * 


Mo 


(l-C2)e”ltTj 


1-G 


2 

ForST0  =  It  the  secular  equation  give  k  =  0.  The  analysis  must  be 

repeated  for  this  case.  We  easily  find 

b(r1  -r  V0 
E“  -  - 

*♦  brl'M0 


U  b(T1  -T)/^ 

!+bV*o 


This  yields 


bTl^o 

A  »  — — — 

1+bTiV0 


T  »  1 


bv*o 

UbV"o 
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APPENDIX  E 

CALCULATION  OF  SOLID  ANGLE  AS  RELATED  TO  A  LAMBERTIAN  SURFACE 


Figure  E.l  shows  the  illustration  of  the  solid  angle  and  its 
representation  in  polar  coordinates.  Under  this  formulation  the  solid 
angle  is: 


fl  = 


where 

o—  -  Surface  area  on  the  spherical  surface 
R  -  Radius  of  the  sphere 
From  the  diagram,  it  an  be  seen  that: 

do-  =  (R  d®  )  (R  sin#  d</>  ) 


n 


d  a~ 


R 


2 


sin  9  d  9  d  o 


the  irradiance  is  defined  by  the  normal  component  of  the  radiance 
integrated  over  the  entire  solid  angle  and  may  be  written: 


E 


L  (0,0) 


cos  9  d  12 


E 


V  2 

L  (9,o) 


cos  9 


s  in  9 


d  9 


d  v 


o  o 

For  isotropic  radiation,  L  is  not  a  function  of  angle  and  the  value  of 
the  irradiance  becomes: 

E  =  t  L 

When  a  Lambertian  surface  is  assumed,  the  surface  is  presumed  to  reflect 
isotropically  in  all  directions.  Hence,  the  solid  angle  is  weighted  by 


219 


the  cosine  factor  as  above.  To  determine  the  radiance  of  a  Lambertian 
surface,  multiply  the  incident  irradiance  by  the  surface  reflectivity 
and  divide  by  t  . 


L 


pE 


Figure  E.l 

Illustration  of  a  solid  angle  and  its  representation  in  polar 
coordinates.  Also  shown  is  a  pencil  of  radiation  through  an 
element  of  area  d A  in  directions  confined  to  an  element  of 
solid  angle  dQ  . 


GLOSSARY  OF  TERMS 


Terms  with  a  firsc  or  second  subscript,  \  ,  are  spectral  quantities  or 
functions  of  wavelength. 

2 

A  -  Area  of  (object)  in  m 

a  ,  -  Absorption  coefficient  for  ozone 

oA 

a  ,  -  Absorption  coefficient  for  the  uniformly  mixed  gases 

uA 

a  .  -  Absorption  coefficent  for  water  vapor 

wA 

-  Correction  factor  for  Bird  model  formulation  of  the  diffuse 
radiation  on  a  horizontal  surface 

-  Day  of  the  year  (Jan  1  =  0;  Dec  31  =  364) 

2 

E  -  Irradiance  on  (object)  in  (W/m  ) 

F  -  Forward  to  total  scattering  ratio  for  the  aerosol 

cl 

2 

H  ,  -  SoLar  irradiance  at  the  top  of  the  atmosphere  (W  /  m  ) 

oA 

I  -  Intensity  in  (W/Sr) 

1^  ^  -  Aerosol  scattered  component  on  a  horizontal  surface 

I,  -  Direct  normal  component  of  solar  radiation 
dA 

I  -  Rayleigh  scattered  component  on  a  horizontal  surface 
rA 

I ^  -  Scattered  or  diffuse  component  on  a  horizontal  surface 

2 

L  -  Radiance  of  (object)  in  (W/m  -Sr) 

M  -  Relative  air  mass  in  atmosphere 

or  Reflected  irradiance  from  (object)  in  (W/m^) 

M  -  Relative  air  mass  for  ozone  in  atmosphere 
o 

M'  -  Pressure  corrected  relative  air  mass  in  atmosphere 
0^  -  Ozone  amount  in  a  vertical  path  (cm) 

P  -  Surface  pressure  in  millibars  (mb) 

P  -  Reference  surface  pressure,  1013.0  mb 


Actual  Earth-Sun  distance 
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R 

m 

R 

s 


T 

WA 

w 

o 

w 


a 


n 

fin 

0 


or 


T 

P  - 


rg 

?S 


Mean  Earth-Sun  distance 

The  distance  between  the  target  location  and  the  satellite 
location 

Transmittance  function  for  aerosol  extinction 
Transmittance  function  for  ozone  absorption 
Transmittance  function  for  the  uniformly  mixed  gases 
(Oxygen  and  carbon  dioxide) 

Transmittance  function  for  water  vapor 
Single  scattering  albedo 

Precipitable  water  in  a  vertical  path  (cm) 

Turbidity  exponent  (Angstrom  formalism) 

Turbidity  coefficient  (Angstrom  formalism) 

Expression  defined  in  terms  of  the  day  of  the  year  for  use  in 

computing  the  astronomical  factors 

Solar  zenith  angle  relative  to  a  flat  surface 

Solar  zenith  angle  relative  to  the  normal  of  a  tilted  plane 

Satellite  zenith  angle  relative  to  the  normal  of  a  tilted 

plane 

Satellite  zenith  angle  relative  to  a  flat  surface 
Wavelength  in  micrometers 
Longitude 
Pi  (3.14159) 

Albedo  or  reflectance  of  (object) 

Albedo  of  the  ground 
Albedo  of  the  air 


<D 


Flux  in  (Watts) 
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